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SUMMARY 


NUCLEAR  WINTER:  WHERE  DO  WE  STAND? 

R.  P.  Turco 

Four  years  of  research  on  nuclear  winter  yhas  greatly  improved 
our  understanding  of  this  complex  phenomenon^  Studies  have 
confirmed  the  possibility  of  significant  temperature  decreases  and 
other  severe  environmental  perturbations  following  a  nuclear  war, 
with  potentially  critical  implications  for  human  survival. 
Nevertheless,  important  uncertainties  remain  to  be  resolved. ,  The 
major  issues  and  their  present  research  status  are  reviewed. 

b^Fuel  Inventories^  Fuel  burdens  in  rural  and  urban  settings 
appear  to  be  known_to_jw4rth in  a  factor  “2,  and  are  lower  than 
earlier  estimates  by  a  factor  of  2  or  so; 

6-»Fuel  Impactiont*  The  quantities  of  fuels  affected  by 
nuclear  explosions  are  sensitive  to  the  scenario  adopted;  although 
definitive  targeting  analyses  are  forthcoming,  issues  of  targeting 
policy  may  never  be  fully  resolved  and  must  necessarily  remain 
uncertain;*' ignition  and  burning  of  rubblized  fuels,  while  likely 
after  a  nuclear  attack,  may  be  of  lower  intensity;  it  is  presently 
understood/ that  a  major  nuclear  attack  against  identified 
mil itaryy industrial  targets  would  cause  extensive  collateral 
damage/ in  urban  areas; 

OOmoke  Emission  Factor:  Burning  petroleum,  plastics  and 
related  materials  can  emit  5%  or  more  of  their  mass  as  soot; 
recent  experiments  reveal  that  the  combustion  of  wood  under 
restricted  ventilation  can  convert  up  to  2%  to  soot;  such  emission 
factors  imply  blacker  smoke  than  has  been  presumed  in  most 
existing  studiesj^soot  emission  factors  for  large  scale 
nuclear-induciecL'-o'i  1  refinery  and  urban  fires  are  still  uncertain, 
but  have  more  likely  been  underestimated  than  overestimated; 

o^piume  Heights:  Simulations  and  observations  indicate 
initial  smoke  injection  as  high  as  15  kilometers;*  soot  generated 
under  intense  flaming  conditions  is  likely  to  be  injected  into  the 
middle  and  upper  troposphere,  while  smoke  produced  by  lower 
intensity  combustion  would  be  deposited  in  the  middle  and  lower 
troposphere; 

o^Prompt  Scavenging:  The  immediate  rainout  of  the  sooty 
component  of  smoke  emissions  is  probably  less  than  50%,  because  of 
its  poor  nucleation  properties  relative  to  other  materials,  and 
because  of  the  likely  overseeding  and  reduced  precipitation 
efficiencies  of  smoke  cloudsf'actual  scavenging  processes  remain 

be  tested  in  full-scale  field  tests; _ _ _  _ 

o^Mesoscale  Dispersion^  Observations  of  large  smoke  plumes  } 
caused  by  fire  complexes  suggest  rapid  regional  dispersion;  the 
role  of  water  vapor  in  modifying  soot  properties  and  initiating 
precipitation  requires  investigation;  smoke  aging  by  coagulation 
and  chemical  reaction  must  also  be  carefully  defined  in  this 
regime; 


o'Acute  Climate  Change:  The  extent  of  land  surface  cooling 
for  a  specific  smoke^  injection  can  presently  be  estimated  to 
within  a  factor  of  r2;  the  primary  unknowns  concern  the  surface 
boundary  layer  response  to  a  sudden  loss  of  solar  energy, 
including  the  formation  of  fogs  and  inversion  layers;  critical 
changes  in  precipitation  rates  may  also  occur,  but  are  currently 
less  well  understood; 

o^Long-term  Climate  Change:  Preliminary  investigations  show 
important  couplings  between  smoke  injections,  oceanic  responses, 
ice  formation  and  other  processes  that  drive  persistent  climate 
fluctuations  over  periods  of  years;  residual  stabilized  smoke 
layers  could  greatly  enhance  these  effects;  uncertainty  about  the 
long-term  chemical  interaction  between  soot  and  ozone  must  be 

resolved|j  . 

brother  Environmental  Effects*  Substantial  ozone  layer 
depletions  are  expected  following  4  nuclear  war  due  to 
smoke-induced  perturbations  in  atmospheric  temperatures, 
composition  and  circulation,  but  the  effect  has  not  been 
quantified;  radiological  impacts  would  be  severe  over  large  areas 
because  of  local  fallout,  although  global  health  effects  need  to 
be  more  clearly  defined;  release  of  chemical  toxins  leading  to 
hazardous  local  pollution  of  air,  soil  and  water  must  be 
quantified;  ■ 

o'Biological  Impacts:  The  human  consequences  of  nuclear  war 
as  described  in  the  SCOPE  Report  must  be  extended  and  refined; 
an  organized  research  program  should  be  developed  in  the  U.S.  to 
place  the  biological  conclusions  on  a  firmer  analytical 
foundation . 


PREFACE 


The  Defense  Nuclear  Agency  has  collected  and  printed  the  attached  papers 
from  the  April  7-9  1987  Global  Effects  review  as  a  service  to  the  community. 

The  Defense  Nuclear  Agency  takes  this  opportunity  to  express  its  gratitude  to 
the  numerous  participants  in  the  Global  Effects  review. 

The  technical  papers  enclosed  include  all  those  which  were  received  by  DNA 
prior  to  the  closing  date  of  15  May  1987.  Where  papers  are  missing  their 
place  is  occupied  by  the  abstract  received  prior  to  the  meeting. 

The  inclusion  of  a  paper  in  this  proceedings  does  not  necessarily  imply 
endorsement  of  the  results  of  the  research  reported  or  conclusions  which  might 
be  drawn  from  that  research.  It  is  the  opinion  of  the  Defense  Nuclear  Agency 
that,  while  good  progress  is  being  made  in  improving  our  understanding  of 
Global  Effects,  the  results  to  date  are  tentative  and  preliminary  and  should 
not  be  used  for  planning  beyond  the  planning  of  future  research. 
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ANALYTIC  EQUILIBRIUM  PREDICTIONS 
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Soot 


Mass  density 
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CURRENT  CALCULATIONS 


1-0  RADIATIVE  CONVECTIVE 
HEAT  TRANSFER  SIMULATION 


RADIATIVE  TRANSPORT 
2  STREAMS 

1  GROUP  VISIBLE,  4  GROUPS  INFRARED 


CONVECTION 

-  6.5*K/KIL0METER  LAPSE  RATE  IS  MAINTAINED 
BY  EXCHANGING  MASS  8ETWEEN  LAYERS 


HATER  LATENT  HEATS  FUSION,  EVAPORATION 

BLACK  BODY  (ZERO  HEAT  CAPACITY)  OR 
CONSTANT  TEMPERATURE  EARTH  SURFACE 


SOOT  CLOUD  1-D  SIMULATION 
INITIAL  CONDITIONS 


MASS  DENSITIES  (GM/M**3) 


SOOT  CLOUD  AT  3,5  DAYS 


SOOTY  WATER  CLOUD 
1  -D  SIMULATION 
INITIAL  CONDITIONS 


SOOTY  WATER  CLOUD  AT  3.5  DAYS 


CONCLUSIONS 


EFFECTS  OF  CONDENSED  WATER  ON  ATMOSPHERIC 
HEATING  BY  SOOT  CLOUDS  WHICH  HAVE  BEEN 
SUPPORTED  BY  NUMERICAL  CALCULATIONS 


WARMER  EARTH  SURFACE  TEMPERATURES 

COOLER  CLOUD  TEMPERATURES 

REDUCTION  IN  CLOUD  RISE 

(THIS  MAY  OCCUR  TO  A  LESSER  EXTENT 

FROM  WATER  VAPOR) 


ON-GOING  WORK  WILL  EXAMINE  THESE  EFFECTS  OVER 
THE  RANGE  OF  CLOUD  DENSITIES  OF  INTEREST 
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Simulation  of  Coastal  Flow  Fields 
when  the  Incident  Solar  Radiation 

Is  Obscured 

by 


3^.  'M.o-lmkcbm ft 

Lawrence  Livermore  National  Laboratory 


SIMULATION  OF  COASTAL  FLOW  FIELDS 


WHEN  THE  INCIDENT  SOLAR  RADIATION  IS  OBSCURED* 


Charles  R.  Molenkamp 

Lawrence  Livermore  National  Laboratory 
Livermore,  CA  94550 


ABSTRACT 


In  the  aftermath  of  a  large  scale  nuclear  exchange,  smoke  from  many 
intense  fires  would  be  injected  into  the  upper  troposphere  and  stratosphere 
in  such  amounts  that  very  little  incident  solar  radiation  would  reach  the 
ground.  It  has  been  suggested  that,  in  coastal  regions,  the  air  over  land 
would  cool  much  more  rapidly  than  that  over  the  ocean.  In  response  to  this 
thermal  gradient,  a  super  land  breese  would  form,  lifting  the  moist  air 
over  the  ocean  and  producing  a  zone  of  anomalous  persistent  precipitation. 

Using  an  enhanced  version  of  the  Colorado  State  University  Mesoscale 
Model,  the  development  of  coastal  flow  fields  are  being  simulated  during 
periods  of  extended  obscuration  of  solar  energy  by  completing  one  normal 
diurnal  cycle  starting  at  dawn  and  then  extending  the  simulation  for  an 
additional  24-48  hours  with  no  incident  solar  radiation.  For  these  simula¬ 
tions  two  changes  have  been  made  to  the  model.  First,  a  cloud  module  has 
been  added  that  allows  for  the  formation  of  a  cloud  or  fog  whenever  the 
atmosphere  is  saturated.  Second,  the  longwave  cooling  module  has  been 
rewritten  to  improve  its  accuracy  and  allow  for  the  effects  of  clouds. 

For  a  typical  west  coast  situation,  with  a  moderate  westerly  mean 
flow,  the  surface  layer  over  land  cools  a  few  degrees  during  the  night  to  a 
temperature  where  saturation  occurs  and  a  thin  ground  fog  forms.  This  fog 
protects  the  ground  surface  from  further  rapid  cooling  and  grows  in  height 
as  a  reasonably  well-mixed  layer  forms  over  land  with  a  cloud  layer  at  the 
top.  Over  the  ocean,  moisture  is  vertically  mixed  upward  into  layers  that 
are  cooling  radiatively  so  a  similar  cloud  forms.  With  clouds  over  both 
land  and  sea  there  are  only  small  thermal  gradients  so  the  hypothesized 
precipitation  zones  do  not  occur. 

Simulations  for  a  typical  east  coast  situation  are  currently  underway. 


‘This  work  was  performed  under  the  auspices  of  the  U.S.  Department  of 
Energy  by  the  Lawrence  Liveraore  National  Laboratory  under  Contract  W- 7405- 
Eng-48. 
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prevents  the  development  of  large  thermal  gradients 
and  significant  induced  flows. 


-  Surface  heat  budget  (land  and  sea  surfaces) 
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-  Sea  surface  temperature  is  289K 


MODIFICATIONS  TO  THE  CSU  MESOSCALE  MODEL 
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Clouds  act  as  boundaries  for  the  clear  air  flux 
calculations 


POTENTIAL  TEMPERATURE 
AT  A  POINT  27KM  INLAND 


Mesoscale  Aspects  of  the  Climatic  Effects  of  a  Large  Nuclear  War 


\ 

I  Filippo  Giorgi 


National  Center  for  Atmospheric  Research 
P.0.  Box  S000,  Boulder,  Colorado  80S07 


)  Most  assessments  of  the  climatic  consequences  of  a  major  nuclear  conflict  have  been 

|  made  using  1-dimensional  (ID)  radiative-convective  models  or  3-dimensional  general  cir- 

f  culation  models  (GCM).  In  these  studies,  initial  distributions  of  nuclear  smoke  and  dust 

*  uniform  over  large  scale  domains  were  employed.  The  assumption  was  that  the  thousands 

of  smoke  and  dust  plumes  generated  in  the  aftermath  of  the  conflict  would  rapidly  merge, 
!  due  to  the  atmospheric  winds,  to  form  a  widespread  cloud.  A  few  preliminary  studies, 

j  however,  have  indicated  that  atmospheric  perturbations  can  be  induced  by  the  nuclear 

(  aerosols  on  the  mesoscale  (10-1000  km),  which  can  strongly  affect  the  aerosol  distribu- 

j  tion  and  properties.  This,  in  turn,  will  influence  both  local  or  regional  “acute”  effects 

(by  perturbing  the  local  meteorology)  and  long  term  “chronic”  effects  (by  modifying  the 
|  amount  and  properties  of  the  injected  material).  In  this  study  it  is  proposed  to  employ  a 

|  mesoscale  model  to  investigate  behaviour,  interactions,  and  meteorological  effects  of  the 

aerosol  plumes  in  the  first  several  days  (5-10)  after  the  conflict.  Results  will  be  compared 
with  GCM  predictions  of  acute  climatic  effects,  and  the  implications  on  the  initial  smoke 
distribution  and  properties  to  be  used  in  GCM  studies  will  be  analyzed.  An  augmented, 
;  2D  version  of  the  Penn  State/NCAR  mesoscale  model  will  be  initially  used.  A  number  of 

;  physics  parameterizations  will  be  added  to  the  current  version  of  the  model.  These  include 

>  a  detailed  description  of  surface  physics  and  radiative  transfer,  an  improved  description  of 

\  the  atmospheric  water  cycle,  and  aerosol  transport,  microphysics,  and  interactive  radia¬ 

tive  effects.  The  sensitivity  of  the  results  to  model  resolution  and  physical  domain,  physics 
parameterizations,  and  aerosol  injection  scenarios  will  be  examined.  Depending  on  the 
indications  of  these  2D  simulations  the  possibility  of  performing  3D  simulations  will  be 
considered. 
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MESOSCALE  ASPECTS  OF  THE  CLIMATIC  EFFECTS  OF  A 

LARGE  NUCLEAR  WAR 


MAIN  OBJECTIVES 

•  Study  of  regional,  "acute” .  climatic  effects  of  massive 
aerosol  injections  from  a  large  nuclear  war  (5-10  days  after 
the  conflict) 

-  Prediction  of  atmospheric  perturbations  on  the 
mesoscale  (10-1000  km),  and  related  perturbations  of 
surface  climatic  variables  (e.g.,  surface  temperature) 

-  Comparison  with  general  circulation  model  (GCM) 
predictions 

•  Study  of  the  effect  of  aerosol-induced  mesoscale  atmo¬ 
spheric  perturbations  on  the  aerosol  properties  and  distri¬ 
bution;  implications  for  injection  scenarios  used  in  global 
scale  studies. 
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MODELING  APPROACH 


•  Penn  State/NCAR  mesoscale  model,  version  MM4 
(2D  version  for  now) 

-  Primitive  equations 

-  a  vertical  coordinate 

-  Hydrostatic 

-  Compressible 

•  Physics  included  in  the  current  version  of  the  model 

-  Moisture 

*  Cumulus  parameterization  (Anthes,  1977) 

*  Explicit  liquid  water  scheme  (Hsie  et  al.,  1984) 

-  Boundary  layer 

*  Bulk  PBL  (Deardorff,  1972) 

*  High  resolution  PBL  (Zhang  and  Anthes.  1982) 

-  Radiative  transfer  (simple  parameterization) 

-  Ground  temperature  calculated  via  an  energy  balance 
equation  (Blackadar.  1979) 

-  Horizontal  and  vertical  sub-grid  scale  eddy  diffusion 
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Model  enhancements  planned  for  this  study 

-  Improve  surface  physics  description 

-  Improve  radiative  transfer  calculations 

-  Include  simple  ice  physics  parameterization 

-  Include  aerosol  transport,  microphysics,  and  interac¬ 
tive  radiative  effects 


SURFACE  PHYSICS  MODEL  ENHANCEMENTS 


•  Apply  to  MM4  the  surface  physics  package  developed  by 
R.  Dickinson  and  collaborators 


SENSIBLE  RADIATIVE  EXCHANGE  MOISTURE 


RADIATIVE  TRANSFER  MODEL  ENHANCEMENTS 


Apply  to  MM4  the  radiative  transfer  package  of  the  Commu¬ 
nity  Climate  Model  (version  CCM1).  This  package  includes: 
Clear  sky 

*  Longwave  — *  H20,  CO2.  O3 

*  Shortwave  — ►  Rayleigh  scattering,  absorption  by  O3. 

H2O.  CO2.  O2 

Cloud-radiation  interactions 

*  Longwave  — *  Clouds  are  assumed  to  be  black  for 
G  >0.5,  for  G  <0.5  clouds  are  assumed  to  be  black 
only  if  the  cloud  water  content  exceeds  10  g/m2 

*  Solar  — ►  Cloud  albedo  depends  on  zenith  angle  and 
cloud  type.  Accounts  for  multiple  reflection  between 
different  cloud  layers  and  between  clouds  and  surface. 
Gaseous  absorption  within  clouds  included. 


MOISTURE  PHYSICS  MODEL  ENHANCEMENTS 

•  A  prognostic  equation  for  cloud  ice  is  added  to  those  for 
cloud  water  and  rain  water  present  in  the  current  explicit 
liquid  water  scheme  (after  Rutledge  and  Hobbs.  1983.  1984) 


AEROSOL  PHYSICS 


•  Prognostic  equations  are  solved  for  the  hygroscopic  (in¬ 
cludes  also  large  wettable  particles  and  IN)  and  hydropho¬ 
bic  fractions  of  a  given  nuclear  aerosol  (smoke,  dust).  The 
equations  include  transport  and  several  microphysical  pro¬ 
cesses.  For  both  aerosol  fractions  the  size  distribution  is 


prescribed. 
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EXPERIMENTS  PLANNED 


•  Test  the  new  physics  parameterizations;  use  as  a  "control 
case"  the  sea-breeze  simulations  of  Yan  and  Anthes  (1987) 

•  Sensitivity  to  injection  scenarios;  use  results  from  plume 
dynamics  simulations 

-  Plume  distribution 

-  Injection  height 

-  Properties  of  injected  aerosol 

-  Aerosol  loading 

•  Sensitivity  to  model  physics  options 

-  Surface  physics  options 

-  Radiation  transfer  options 

-  Moisture  physics  options 

•  Sensitivity  to  model  domain  and  resolution 

-  Surface  characterstics 

-  Seasonal  effects 

-  Initial  conditions 

-  Boundary  conditions 
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A  CASE  STUDY  OF  A 
FOREST  FIRE  SMOKE  PLUME 


Douglas  L.  Westphal 

Pennsylvania  State  University 
Owen  B.  Toon 
NASA  Ames  Research  Center 


PROJECT  SUMMARY 


•  Goals 

—  Understand  the  regional-scale  interactions  between 
large  soot  clouds  and  atmospheric  circulations 

—  Determine  the  impact  of  soot  clouds  on  local  cli¬ 
mate 


•  Method 

—  Study  natural  analogs  to  nuclear  winter 

—  Develop  numerical  models;  verify  their  accuracy 
by  modeling  the  natural  analogs 

—  Use  the  models  to  study  possible  nuclear  winter 

l 

;  scenarios 

i 

* 

> 

> 
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FORT  NELSON  FIRE  CHRONOLOGY 

June  19  -  Fire  started  by  lightning 
400  ha  (100  ha  =  1  km2). 

June  24  -  20,000  ha 

July  15  -  78,000  ha 

July  26  -  87,000  ha 

July  28  -  88,000  ha 

July  29  -  160,000  ha  (40  km  x  40  km) 

Fire  doubled  in  size  in  24  hrs. 

16  km  run  in  4  hours. 

Convective  column  was  observed. 

(~  0.07  to  0.4  Tg  soot  produced) 


FORT  NELSON  FIRE  DATA  SUMMARY 


•  Fire  Data 

—  Fire  spread  maps,  fuel  loading  (Lawson,  Canadian  Forest 
Service) 

—  Lodi  Canyon  fire  and  Fresno  area  fire  (Peuschel  et  al.) 

—  Landsat  and  AVHRR  satellite  imagery 

•  Meteorological  Data 

—  Surface  observations 

—  Upper  air  (NMC  analyses  and  Navy  data  archives) 

—  MOS  analyses  (Robock,  U.  of  Mary.) 

•  Aerosol  Data 

—  Aircraft  pilot  reports 

—  Surface  radiometric  measurements  (Ryznar,  U.  of  Mich.) 
—  Satellite  data  (Fraser  and  Ferrare,  Goddard) 

—  Nephelometer  measurements  (Eloranta,  U.  of  Wise.) 


CONCLUSIONS 

Significant  dataset  is  available 

—  Source  term  information 

—  Soot  cloud  optical  properties 

—  Observations  of  perturbations  due  to  smoke  cloud 

Transport  modified  by  small-scale  wave 
Model  shows  skill  at  forecasting  small-scale  waves 
Ready  to  carry  out  interactive  calculations 


FIGURE  CAPTIONS 

Fig.  1  Outline  of  the  area  where  surface  meteorological  stations  re¬ 
ported  “visibility  reduced  by  smoke”  during  the  period  4  am 
(PST)  29  July  to  4  am  30  July,  1982. 

Fig.  2  Contours  of  the  analyzed  (~  observed)  height  (gpm)  of  the 
500  mb  surface  at  4  pm  (PST)  31  July,  1982. 

Fig.  2  Contours  of  the  analyzed  (~  observed)  height  (gpm)  of  the 
500  mb  surface  at  4  pm  (PST)  2  August,  1982. 

Fig.  4  Contours  of  the  predicted  height  (gpm)  of  the  500  mb  surface 
at  4  pm  (PST)  2  August,  1982.  Compare  with  Fig.  3. 

Fig.  5  Contours  of  the  analyzed  (~  observed)  surface  temperature 
(°C)  at  4  pm  (PST)  2  August,  1982. 

Fig.  6  Contours  of  the  predicted  surface  temperature  (°C)  at  4  pm 
(PST)  2  August,  1982.  Compare  with  Fig.  5. 
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Microphysical  and  Radiative  Simulations  of  Optically  Thick  Smoke  Clouds 

on  the  Continental  Scale 

Owen  B.  Toon1 ,  Douglas  Westphal 1  *2,  Krishnamurthy  Santhanam3 

1  NASA  Ames  Research  Center  Moffett  Field  Ca  94035 
^Pennsylvania  State  University ,  University  Park,  Pa.  16802 
■^Sterling  Software,  Palo  Alto  Ca.  94303 

We  are  using  the  Penn  State/NCAR  regional  scale  dynamical  model  in 
conjunction  with  radiative  transfer  and  aerosol  microphysical  models 
developed  at  Ames  Research  Center  in  order  to  simulate  optically  thick 
smoke  clouds  as  they  evolve  over  several  day  time  periods.  Our  eventual 
goal  is  to  study  the  interactions  of  the  plumes  from  many  large,  separate 
fires  representing  the  results  of  a  nuclear  conflict.  However,  we  first 
wish  to  establish  the  ability  of  the  model  to  simulate  the  plume  from  a 
single  large  fire.  We  have  therefore  chosen  to  model  the  July- August  1982 
smoke  plume  which  originated  over  British  Columbia  and  spread  over  the 
East  Coast  of  the  United  States.  The  particle  size,  optical  depth,  single 
scattering  albedo,  horizontal  extent  of  the  cloud  and  effect  of  the  cloud 
on  the  surface  radiation  budget  and  surface  temperature  have  all  been 
determined  by  a  variety  of  researchers  from  ground  based,  satellite  and 
aircraft  observations.  We  are  presently  assemblying  a  comprehensive 
set  of  data  on  this  smoke  cloud.  Here  we  present  two  dimensional 
simulations  of  the  cloud.  We  use  a  wind  profile  and  water  vapor 
concentration  observed  at  one  station  as  a  representative.  We  then  advect 
smoke  over  a  several  day  time  period.  Our  goals  are:  to  determine  the 
relative  importance  of  vertical  wind  shear  and  duration  of  the  fire  for  the 
geographic  area  covered  by  smoke  at  any  one  time;  to  determine  the 
Importance  of  coagulation  to  smoke  particle  size;  and  to  estimate  the 
radiative  heating  rates  in  the  cloud  and  the  effects  of  the  cloud  on  the 
surface  energy  budget.  These  preliminary  calculations  will  be  used  to 
guide  more  sophisticated  three  dimensional  interactive  calculations. 


REGIONAL  SCALE  SIMULATIONS  OF  OPTICALLY  THICK 

SMOKE  CLOUDS 


GOAL:  Perform  interactive  simulations  of  dense  smoke 
clouds  including  dynamics,  radiative  transfer,  aerosol 
physics  and  water  vapor  condensation. 

Steps: 

1.  Identify  natural  test  case  and  gather  data. 

% . 

2.  Perform  3-d  non-interactive  simulations  of 
meteorology. 

3.  Perform  2-d  simulations  of  aerosol  physics  and 
radiative  heating  to  determine  which  processes  are 
important. 

4.  Perform  3-d  non-interactive  simulations  of  aerosol 
physics  and  radiative  transfer. 

5.  Perform  3-d  interactive  simulations. 

6.  Perform  simulations  of  alternative  fire 
possibilities  (multiple  fires,  shorter  duration  fires, 
more  intense  fires). 


Questions  addressed  by  2-d  aerosol  and  radiative 

simulations 


1.  Does  coagulation  affect  the  optical  properties  of  the 
cloud? 


2.  Are  the  heating  rates  within  the  cloud  significant? 

3.  What  are  the  properties  of  the  initial  cloud  that  are 
needed  to  explain  the  observed  properties 
downwind  (particle  size,  refractive  indices)? 


4.  Is  dry  deposition  significant  in  the  few  day  time 
period? 


Initial  conditions 

1.  Cloud  originally  located  over  12°  of  longitude  at 
altitudes  of  1  to  3  km. 

2.  Initial  size  distribution  chosen  to  reproduce  Lodi 
spectral  optical  depth  data.  (Log  normal  with  number 
mean  size  of  0.06pm,  sigma  =  1 .6). 

3.  Wind  speed  increases  in  proportion  to  the  inverse 
pressure  with  a  value  of  lOm/sec  in  the  middle  of  the 
cloud. 

4.  Vertical  diffusion  is  constant  at  104  cm2  s~l. 

5.  Coagulation  and  sedimentation  assume  unit  density 
spheres. 


Sensitivity  tests 

1.  No  coagulation 

2.  No  vertical  wind  shear 

3.  High  diffusion  (106  cm2  s”1 )  in  lowest  km 


Figure  Captions 


Fig.  1  Panels  a-d  illustrate  the  evolution  of  the  optical 
depth,  mass  concentration,  number  concentration,  and 
extinction  coefficient  for  three  days  of  cloud  dispersal. 
The  particle  concentration  decreases  by  a  factor  of 
about  20  due  to  coagulation  and  the  optical  depth  first 
increases,  due  to  the  particle  size  increases,  and  then 
decreases  due  to  spreading  of  the  cloud.  Optical  depths 
of  about  2  were  observed  in  the  July  82  smoke  plume 
over  the  East  coast  of  the  United  States. 

Fig.2  Illustrates  the  size  distribution  initially  and  the 
size  distribution  after  72  hours.  Coagulation  leads  to  an 
increase  of  size  into  the  range  observed. 

Fig.  3  Illustrates  the  properties  of  the  cloud  after  72 
hours  if  coagulation  is  not  included.  Comparison  with 
Fig.  Id  shows  that  coagulation  leads  to  an  increase  in 
optical  depth. 

Fig.  4  Illustrates  the  cloud  properties  after  72  hours 
without  vertical  wind  shear.  In  this  case  the  optical 
depth  is  about  the  same  as  after  24  hours  (Fig.  lb) 
indicating  that  shear  is  reducing  the  optical  depth  by 
dispersing  the  cloud.  However  further  coagulation 
beyond  24  hours  is  not  affecting  the  optical  depth. 

Fig.  5  Illustrates  the  cloud  properties  after  72  hours 
with  a  boundary  layer  having  high  diffusion  coefficients. 
The  optical  depth  is  somewhat  reduced,  indicating  that 
dry  deposition  may  be  important. 

Fig.  6  Illustrates  the  additional  heating  rate  when  the 
cloud  is  present.  At  noon  (72  hours)  significant  heating 
rates  occur  in  the  smoke  cloud.  At  night  the  cloud  has 
almost  no  effect  on  the  heating  rate. 
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Fig.  7  Part  a  Illustrates  the  surface  net  radiation  budget 
with  no  cloud  while  Part  b  illustrates  the  radiation  with  a 
cloud.  At  noon  the  cloud  reduces  the  solar  energy 
reaching  the  ground  by  almost  50%  and  the  total  energy 
by  about  25%. 


RADIUS  (cm) 


RADIUS  (cm) 


CONCLUSIONS 


1 .  Coagulation  leads  to  an  increase  in  cloud  optical  depth 
over  the  first  few  days. 

2.  Dry  deposition  is  an  important  loss  process. 

3. An  imaginary  refractive  index  in  the  range  of  0.01  to 

0.02  and  an  initial  particle  size  distribution  similar  to 
those  observed  at  Lodi  lead  to  sizes  and  single 
scattering  albedos  after  a  few  days  which  are  close 
to  thos£  observed. 

4.  Significant  heating  rates  occur  in  the  cloud  and  a 
substantial  light  loss  occurs  at  the  surface. 
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SURFACE  TEMPERATURE  EEFECTS  OF  FOREST  FIRE  SMOKE  PLUMES 


Alan  Robock 

Department  of  Meteorology- 
University  of  Maryland 
College  Park,  Maryland  20742 

Paper  presented  at  DNA  Global  Effects  Program 
Technical  Meeting.  Santa  Barbara.  CA.  7-9  April  1987 

ABSTRACT 

The  location  of  smoke  clouds  from  extensive  Canadian  forest 
fires  is  determined  from  satellite  imagery  for  cases  in  the 
summers  of  1981  and  1982.  As  these  smoke  clouds  pass  over  the 
Midwestern  United  States,  surface  air  temperature  effects  are 
determined  by  comparing  actual  surface  air  temperatures  with 
those  forecast  using  Model  Output  Statistics  (MOS)  by  the  United 
States  National  Weather  Service.  A  cooling  of  1.5  to  4°C  is 
found  in  the  daytime  under  the  smoke  plumes.  No  effect  is  found 
at  night.  This  corresponds  to  theoretical  estimates  of  the 
effects  of  an  elevated  smoke  plume.  It  serves  as  observational 
confirmation  of  portions  of  the  nuclear  winter  theory. 

1 .  Introduction 

Crutzen  and  Birks  (1982)  suggested  that  smoke  from  forest 
fires  ignited  by  nuclear  weapons  would  be  extensive  enough  to 
block  out  significant  amounts  of  sunlight.  Subsequent  work, 

such  as  Crutzen  et  al .  (1984) ,  Turco  _ al .  (1984) ,  National 

Academy  of  Sciences  (1985) ,  and  SCOPE  (1986) ,  has  pointed  out 
that  the  smoke  from  urban  and  industrial  fires  would  be  much 


more  effective  at  preventing  solar  radiation  from  reaching  the 
earth's  surface  than  forest  fire  smoke  after  a  large-scale 
nuclear  war.  With  both  urban  and  rural  targets ,  not  only  would 
more  urban  smoke  be  generated,  but  its  optical  properties  would 
make  it  more  effective  at  blocking  sunlight. 

The  effects  of  forest  fire  smoke  are  still  interesting, 
however,  for  two  reasons.  First  of  all,  a  lot  of  forest  fire 
smoke  would  still  be  generated  in  many  nuclear  war  scenarios, 
especially  those  that  include  only  non-urban  military  targets. 
The  optical  properties  and  surface  temperature  effects  of  this 
smoke  is  an  important  part  of  the  study  of  nuclear  winter. 

Secondly,  in  order  to  understand  nuclear  winter,  it  would 
be  nice  to  have  some  actual  observations  of  the  effects  of  smoke 
to  compare  to  the  theoretical  models.  Extensive  contemporary 
urban  and  industrial  smoke  plumes  are  not  available  for  study. 
Each  year,  however,  forest  fires  are  generated  by  lightning,  and 
in  some  cases  can  produce  extensive  smoke  plumes .  Ginsburg  and 
Golitsyn  (1986)  presented  observations  of  daytime  surface 
cooling  from  extensive  Siberian  forest  fires  in  1915.  Wexler 
(1950)  presented  anecdotal  evidence  of  daytime  cooling  at  the 
surface  of  2-5°C  from  "The  Great  Smoke  Pall"  caused  by  extensive 
Canadian  forest  fires.  In  this  paper,  two  such  cases,  in  which 
smoke  from  western  Canadian  forest  fires  drifted  across  the 
midwestern  United  States  are  investigated  to  attempt  to  measure 
their  surface  temperature  effects.  It  is  believed  that  "The 
Great  Smoke  Pall"  and  the  Siberian  fires  of  1915  were  more 
optically  thick  cloud  than  the  ones  investigated  here. 


The  location  of  smoke  clouds  from  extensive  Canadian  forest 


fires  was  determined  from  satellite  imagery  for  cases  in  the 

summers  of  1981  and  1982.  The  smoke  clouds  are  easy  to  detect 

with  visible  wavelength  imagery  (Figs.  2-7,  9),  but  are 

invisible  in  infrared  imagery  (Fig.  8,  which  can  be  compared 

with  the  visible  image  in  Fig.  3),  implying  a  much  greater 

optical  depth  in  the  visible  as  compared  to  the  infrared. 

Within  2-3  days  after  being  produced,  the  smoke  plumes  in  the 

£  2 

1982  case  grew  into  an  extensive  shield  covering  about  10  km  . 
There  was  some  patchiness  evident  (Figs.  2-7),  but  the  smoke 
veil  maintained  a  coherent  structure  for  several  days.  In  the 
1981  case,  the  smoke  appeared  to  be  in  streaks,  with  patches  of 
less  dense  smoke  in  between  (Fig.  9) .  The  total  area  covered 
was  about  the  same  as  the  1982  case. 

As  these  smoke  clouds  passed  over  the  midwestern  United 
States,  surface  air  temperature  effects  were  determined  by 
comparing  actual  surface  air  temperatures  with  those  forecast 
using  Model  Output  Statistics  (MOS)  by  the  United  States 
National  Weather  Service  (Glahn  and  Lowry,  1972) .  This  MOS 
error  technique  has  been  used  successfully  before  to  determine 
surface  air  temperature  effects  of  the  Mount  St.  Helens  volcanic 
eruption  of  1980  (Robock  and  Mass,  1982) .  The  analysis  was  done 
in  regions  of  high  pressure  where  synoptic  disturbances  were  not 
affecting  the  temperature.  The  errors  are  attributed  to  the 
presence  of  aerosols  in  the  atmosphere,  since  the  aerosol 
content  is  not  a  MOS  predictor.  The  locations  of  all  the  MOS 
stations  used  in  this  analysis  are  shown  in  Fig.  1. 


Intense  forest  fires  burned  in  British  Columbia,  Canada, 
on  29  and  30  July  1982  (Eerrare,  personal  communication) . 
Smoke  from  these  fires  was  transported  by  the  prevailing  winds, 
and  crossed  the  U.S.  border  in  North  Dakota  on  31  July.  It 
proceeded  across  the  Midwest  and  then  east  over  the  Mid-Atlantic 
states.  Smoke  from  these  fires  was  reported  over  Western  Europe 
on  5  and  6  August.  Such  was  also  the  case  for  The  Great  Smoke 
Pall  (Wexler,  1950)  and  demonstrates  the  ability  of  the 
atmosphere  to  transport  smoke  over  long  distances  before  it  is 
removed . 

Figures  2-7  show  both  visible  wavelength  GOES  satellite 
images  and  MOS  surface  temperature  forecast  errors  for  31  July 
and  1  August  1982  at  1500,  1800  and  2100  GMT  (10  am,  1  pm  and 
4  pm  CDT  -  local  time  in  the  Midwest)  each  day.  Smoke  plumes 
appear  gray  in  these  images,  and  can  easily  be  detected  over  the 
Midwest.  On  the  original  images,  the  smoke  plume  which  was 
headed  for  Europe  can  easily  be  seen  on  1  August  over  the 
Atlantic  Ocean. 

Forecast  errors  of  -2  to  -4°C  are  evident  under  the  smoke 


plumes.  Forecast  errors  are  not  evident  under  the  smoke  plumes 
at  other  times  of  the  day.  Since  the  errors  are  only  evident  at 
times  of  maximum  daily  insolation,  their  predominant  effect  on 
shortwave  radiation  is  demonstrated. 

Forecast  errors  of  the  same  or  greater  amplitude  are  also 
evident  in  other  locations  in  the  Figures.  Upon  close 
inspection,  each  can  be  attributed  to  the  presence  or  absence  of 
water  clouds  or  mesoscale  features  which  were  not  well  forecast 
by  the  LFM  numerical  forecasting  model,  which  provides  the 
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predictors  for  the  MOS  forecasts.  In  clear  areas,  however,  the 
only  large  errors  are  under  the  smoke  cloud.  Eor  example,  in 
Fig.  2,  for  1500  GMT  on  31  July,  the  negative  area  in  eastern 
Michigan  can  be  attributed  to  an  unforecast  mesoscale  region  of 
cloudiness.  The  positive  error  region  in  western  Tennessee  and 
Kentucky  is  associated  with  a  mesoscale  clear  area  in  an 
otherwise  overcast  region. 

In  Fig.  3,  for  1800  GMT,  31  July,  the  same  error  regions 
are  present  as  in  Fig.  2,  except  that  an  additional  negative 
area  appears  in  Virginia  under  a  heavy  water  cloud  bank,  and 
small  negative  areas  appear  east  and  west  of  Lake  Ontario. 
These  latter  areas  become  more  extensive  and  of  higher  amplitude 
three  hours  later  in  Fig.  4.  They  are  associated  with  strong 
cold  advection  behind  a  rapidly  developing  mesoscale  Low,  the 
center  of  which  is  located  on  the  Vermont -Canadian  border.  This 
is  easily  seen  on  the  synoptic  weather  maps  (not  presented  here) 
and  presumably  was  not  well  forecast  by  the  LFM,  which  does  not 
handle  weather  systems  of  this  small  scale  well.  In  Figs.  3  and 
4  the  moderating  effects  of  Lakes  Erie  and  Ontario  are  evident, 
as  the  maximum  cooling  associated  with  this  cold  advection  is 
over  the  land  areas  between  the  lakes.  The  only  other  negative 
error  area,  and  by  far  the  largest  one,  is  under  the  forest  fire 
smoke  veil.  Areas  of  fair  weather  cumulus  over  Illinois, 
Indiana  and  Ohio  do  not  produce  large  MOS  errors. 

On  1  August  (Figs.  5-7)  the  only  large  error  regions  (both 
negative)  are  associated  with  a  small  region  of  water  clouds 
over  Wisconsin  and  the  forest  fire  smoke  veil.  The  region 
around  the  Great  Lakes  under  the  smoke  veil  has  errors  between 


-1°C  and  -2°C,  evidence  again  of  a  moderating  effect  of  the 
lakes.  The  thickest  patch  of  smoke  appears  over  Baltimore, 
Maryland,  in  Fig.  7  at  2100  GMT,  and  has  an  error  of  -4.1°C 
associated  with  it. 

The  results  presented  here  are  from  forecasts  made  from  LFM 
runs  12  to  24  hours  before  the  forecast  time.  Comparisons  were 
made  with  forecasts  made  with  earlier  and  later  model  runs,  and 
the  results  were  virtually  the  same.  The  errors  (cooling)  found 
were  much  larger  than  the  run  to  run  differences. 

Attempts  were  made  to  estimate  the  height  of  the  smoke  as 
it  passed  over  the  Midwest.  Winds  were  examined  at  different 
levels  during  this  time  and  compared  to  the  motion  of  the 
smoke.  The  winds  corresponding  most  closely  to  the  motion  were 
at  700  mb  (3.2  km) .  The  jet  of  smoke  headed  to  Europe  on  1 
August  corresponded  to  winds  between  the  700  mb  and  500  mb 
levels  (3  -  5.5  km) . 

Eloranta  (personal  communication)  happened  to  be  conducting 
a  vertical  profile  study  of  lower  level  haze  on  the  next  day 
over  the  eastern  shore  of  Chesapeake  Bay,  due  east  of  Baltimore 
at  1700  GMT.  The  portion  of  the  smoke  veil  that  was  over 
Kentucky  on  1  August  appears  from  satellite  images  to  have  been 
over  Maryland  at  that  time.  He  flew  up  into  the  base  of  the 
forest  fire  smoke  plume  at  an  altitude  of  approximately  3.5  km 
and  at  4.5  km  (the  highest  altitude  reached)  he  was  still  in  the 
smoke . 

A  preliminary  analysis  of  the  12  August  1981  case  (Fig-  9) 
also  shows  daytime  cooling  under  the  smoke.  In  this  case,  the 
smoke  is  not  as  uniform,  but  appears  in  streaks.  The  cooling  in 


$;( 

*:*?*:■ 


m 


9 
I 

V 


\7 


this  case  is  -1.5  to  -3C.  A  preliminary  examination  of  the 


satellite  images  from  this  case  shows  shorter  shadows  at  sunrise 


and  cumuli  form  clouds  casting  shadows  on  the  smoke  later  in  the 


day.  These  both  suggest  that  the  smoke  was  at  a  lower  elevation 


on  these  days . 


A  cooling  of  1.5  to  4  C  is  found  in  the  daytime  under 


forest  fire  smoke  plumes  in  two  cases  in  1981  and  1982.  No 


effect  is  found  at  night.  This  corresponds  to  theoretical 


estimates  of  the  effects  of  an  elevated  smoke  plume  (Ginsburg 


and  Golitsyn.  1986)  with  optical  depth  of  approximately  2. 


Ferrare  (personal  communication)  reports  theoretical  estimates 


of  the  optical  depth  of  this  smoke  cloud.  based  on 


multi -spectral  satellite  measurements,  of  2  to  3. 


This  result  serves  as  observational  confirmation  of 


portions  of  the  nuclear  winter  theory.  Smoke  from  burning 


cities,  which  could  result  from  a  large  scale  nuclear  war,  would 


be  a  better  absorber  of  shortwave  radiation  than  forest  fire 


smoke,  and  would  therefore  produce  even  larger  surface 


temperature  drops.  The  relative  amounts  of  forest  fire  and  city 


fire  smoke  resulting  from  a  nuclear  war  would  be 


scenario -dependent . 


Several  more  aspects  of  this  research  are  planned.  These 


include : 


-  complete  study  of  11-13  August  1981  case.  Satellite 


images  and  MOS  forecasts  have  already  been  obtained. 


search  for  other  cases  to  study  that  can  use  the  same 
research  technique.  They  would  have  to  have  extensive 
smoke  veils  over  clear  areas  of  the  U.S.  One 

possibility  is  10  September  1981  (Fig-  6,  Chung  and 
Lee ,  1984) . 

use  a  radiative-convective  model  to  calculate  the 
theoretical  surface  cooling  that  would  result  from  the 
cases  studied  here.  This  would  require  knowledge  of 

the  height  of  the  smoke  cloud  and  the  radiative 

parameters  of  the  smoke  particles.  Ferrare  (personal 
communication) ,  in  collaboration  with  Kaufman  and 

Fraser  at  NASA/GLA  will  complete  within  the  next  month 
a  study  of  the  1982  case  presented  here,  in  which  they 
use  a  multi-spectral  satellite  technique  to  obtain  the 
optical  depth,  single-scatter  albedo  and  mean  mass 
radius  of  the  smoke  particles  at  a  number  of  locations 
on  several  days.  I  intend  to  obtain  the  data  from 
them  and  do  these  calculations. 

I  have  obtained  data  from  the  U.S.  network  of  surface 
radiation  stations  for  the  1982  case,  and  will  search 
them  for  information  which  may  be  useful  to  compare  to 
the  Ferrare  data. 

work  with  Toon  and  Westphal  (personal  communication)  who 
are  using  a  mesoscale  dynamic  and  radiative  model  to 
study  the  1982  case.  I  will  compare  the  location  of 
the  smoke'  plume  and  the  surface  temperatures  that  I 
have  obtained  to  those  that  they  calculate.  This  will 
help  to  verify  their  model  and  suggest  other 
observational  studies  of  this  case. 
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MODELS 


GLODEP2:  Estimates  gamma  ray  dose 
from  global  fallout. 
Authors:  L.  Edwards,  T.  Harvey. 

K.  Peterson 

Reference:  UCID-20033  (1984) 
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GRANTOUR:  Calculates  the  3-D  transport, 

6MdfUU-\  diffusion,  and  wet  and  dry 

deposition  and  dose  from  trope- 
(1  ^  spheric  global  fallout. 

^U^LM  Authors:  J.  Walton  and  M.  MacCracken 
^  Reference:  UCID- 19985  (1984) 

KDFOC2:  Calculates  local  fallout. 

Authors:  T.  Harvey  and  F.  Serduke 
Reference:  UCRL-5 2858  (1979) 
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BECAUSE  INJECTIONS  OCCUR  IN  POLAR  OR  EQUATORIAL  ATMOSPHERES , 

THERE  ARE  EIGHT  COMPARTMENTS 


Global  Fallout  Dose  Assessments  (Rads)  for  an 
Unperturbed  Atmosphere  With  No  Smoke 

Knox  (1983)  5300  Mt  baseline  nuclear  war 


Latitude  band 
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Global  population  dose 
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(xlO^O)  person-rads 
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A  a  Summer  injection  using  GL0DEP2 

B  a  Summer  injection  using  GRANTOUR  (run  for  30  days)  with  stratospheric 
contributions  from  GL0DEP2 
C  a  Winter  injection  using  CL0DEP2 

D  a  Winter  injection  using  GRANTOUR  (run  for  30  days)  with  stratospheric 
contributions  from  GL0DEP2 
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Global  Fallout  Oosa  Assessments  (Rads)  using  GRANTOUR; 
Comparing  Perturbed  Atmosphere  (Smoke)  and 
Unperturbed  Atmosphere  (No  Smoke) ; 

Land  Areas  Only 
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The  doses  here  include  the  stratospheric  contributions  ss  calculated  by 
GLOOEP2  and  GRANTOUR  (run  for  30  days) 


) 


\ 


101 


Lo  CAL  FfiLldu  T 


'CWS  ON  lOWPL&X 

rmrmt  co>u>/r*V& 


>  l&.llod' 
far*. 


6-pFecriu<r  cutfb  »& 

*tTwtT'r-%a«r  onTRiatfitoM 

-»lkrwirY-  H«n*rtr  &iir*i»«noMi 


50  km 


& 

M 

$3 

m 

m 

& 


3*%t»whpi  an* 

«o*?*r  fA<TM.  «r  2, 

utTHF&io  excGArw/i 


^fitHCr^e  $ 


L'  -rn-fZ  fecceA  JT 


'«?! 


?OH-A)/oy  boY 


JV\»GU=  u 


'Teapot  £s$ 


33 


94u/uy  -  0oY 


cpQPjo^a.  r 


ulttf 


£'/Ac.UAr*o*4  4  SttGLfV&fAjG.  /.  *-<r 

ft  <r*/t  FtC4tiy— 


6A mma  ooie.  tectn/eo  ttot>s) 


UWJ  CUJ>  W  \AAtL$&r  hj 

'dM. 


flu? 


1  hour  to  1  month  20-Knot  Wind 
4C0-*  Contour  1-MtonYiaid 


OifUnca  Downwind  (Milat) 


Tho  400- ft  iaodoaa  conaoun  for  on«  hour  to  on* 
Md  for  Wtaut  from  ■  lOOOMWiot  reactor, 
1-Mton  Mon  weapon,  wtd  combination. 


riim  1  1,000  MW(t)  LMTBR  and  1-Mton 


and  Cheater  (Raf.  IS) 


(dark  talar).  liman 
Taw  altar  Ik*  attack 


10  rm  par  year  la 
•  nactar  la  SSjOOO 


i  ntakabttad  for  a 


115 


Trrrrrr 


o  O  -  PI  N 

8  5  o  o  o  -• 

6  dodo  d 


V, 

w 

vl 

O 

:■$ 


I 


$ 

$ 

5: 

n 

& 

»»:M 

a 

w»a 

•X- 

$ 

*,yJ 


Vo 


C  M  u 

.2  3,- 

«*§  = 
c.S  < 

■2  Si 
*|! 
3  -S 

<*  «  M 

-awe 

-«•= 
3  V 

e  c  -r 

«  O  -S 
« —  . 
®  «  «tf 

J  .s 
-  o  t 

—  g» 

A  §  3 
6  *H  •» 
|  -  * 
<b  T3  >» 

o»  e  c 

- 

*  «  - 
w  e  c 

>»  5  3 

6  o-  3 
•*  3  5 

e  *  t 

i:  v  «* 

^-S.S 


-  1  2 
*  -S  O 

td  w  -3 

-as* 

o,2  3 

*3 

H  -5  > 


to 

o  o  -<  f  PI  N 

^  dodo  — 


«0  «  O  <D 

o  6  -<  o 


o 

« 

o 

o 

o‘ 

*4 

« 

<o 

« 

d 

o 

pi 

o» 

q 

« 

d 

ri 

■4 

—4 

o» 

N 

pi 

p> 

ai 

pi 

p» 

p> 

n 

oo 

N 

fi 

N 

%» 

P4 

•P> 

V 

p- 

-r 

d 

S  *  p-  'ft  c, 

5  £  2  £  as 

•/}  >  >J  '/)  u. 


116 


‘.'■vtaT 

-  , »  ...» .  • 


r  '»f.v  y  .v  ■ 

■■t  ,'j  . 


**  *’• 


Forage-cow-milk, 

Pasture-  meat  (including  steer,  lamb,  and  hog 
meat). 

Leafy  vegetable. 

Grains  and  other  non- leafy  crops, 

Fresh  and  salt  water  fish,  and  other  edible 
aquatic  animals. 


Surface  air  concentration  of  the  radionuclide. 
Deposition  of  Cs-137  on  the  grass  or  dry  feed. 
Fraction  of  the  nuclide  that  remains  on  the 
forage, 

Residence  time  of  the  nuclide  on  the  forage, 
Intake  rates  of  both  wet  (grass)  and  dry  (hay) 
feed  consumed  by  the  cow, 

Concentration  in  milk. 

Dose  per  unit  activity  ingested, 

Infant's  milk  consumption  rate. 


(a)  little  or  no  effect, 

(b)  slight  to  moderate  effect, 

(c)  moderate  to  large  eftacd. 

(d)  large  effect. 
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INTRODUCTION 


[Slide  11  Because  of  continuing  interest  in  acute 
radiobiological  effects  of  ionizing  radiation,  we  reviewed 
the  available  literature  to  develop  a  symptomatology  basis 
for  assessing  early  functional  impairment  of  individuals  who 
may  be  involved  in  civil  defense,  medical  care,  and  various 
military  activities  in  the  event  of  nuclear  attack. 

This  paper  will  summarize  our  information  sources  and 

describe  how  we  distilled  from  them  descriptions  of  human 

* 

signs  and  symptoms  for  prompt  ionizing  dose  in  the  range  50 
to  3000  rads  (cGy)  raidline  body  tissue. 

Before  I  get  started,  let  me  first  mention  that  the 
authors  would  like  to  acknowledge  the  Radiation  Policy 
Directorate  of  the  Defense  Nuclear  Agency,  in  particular 
Drs.  David  Auton  and  Robert  Young,  for  supporting  our  ef¬ 
fort. 


This  review  was  undertaken  as  a  group  effort  where 
varying  individual  interpretation  of  frequently  imprecise  or 


* 

"Symptoms"  is  used  to  mean  both  subjective  evidence 
and  objective  signs  of  radiation  sickness. 


lacking  information  was  resolved  by  arriving  at  a  consensus 
in  constructing  our  representations  of  acute  radiation 
sequelae.  Rather  than  restating  the  well  known  variability 
in  acute  human  response  to  radiation,  specifically,  our 
objective  was  to  analyze  the  available  data  and  correlate 
radiation  sickness  symptoms  with  dose  levels  and  time  in 
order  to  develop  a  consensus  about  typical  symptomatology 
for  incidence,  onset  time,  severity,  and  duration. 


SYMPTOMATOLOGY  OF  ACUTE  RADIATION  EFFECTS 
IN  HUMANS  AFTER  DOSE  OF  50  TO  3000  RADS  (cGy) 
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SOURCES  OF  INFORMATION 


[Slide  2]  Data  on  human  radiation  sickness  signs  and 
symptoms  (for  convenience  we  simply  refer  to  symptoms)  are 
both  diverse  and  sparse.  Our  information,  representing  a 
review  of  well  over  100  references,  fall  into  four  general 
categories:  (1)  case  studies  of  nuclear  accident  victims, 

(2)  records  of  total-body  radiation  therapy  patients,  (3) 
"composite"  analyses  based  on  data  from  a  variety  of 
sources,  and  (4)  expert  opinion,  sometimes  elicited  in 
private  communication.  Additional  information  comes  from 
survivors  of  the  atomic  bombings  in  Japan  and  those  acciden¬ 
tally  irradiated  in  nuclear  tests  in  the  South  Pacific. 

No  one  source  category  provides  a  comprehensive  picture 
of  the  incidence,  severity,  and  duration  of  radiation  sick¬ 
ness  in  humans.  For  example,  the  data  on  atomic  bomb  sur¬ 
vivors  are  usable  for  delayed  hematological  effects  (  1  week 
postexposure)  but  are  more  uncertain  for  early  acute  ef¬ 
fects;  records  of  symptoms  during  the  first  few  postexposure 
days  were  constructed  some  time  after  the  fact.  To  refine  a 
plausible  description  of  acute  symptoms,  we  evaluated  the 
data  from  a  variety  of  sources.  The  handouts  contain  a  list 
of  our  references. 
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With  the  next  few  slides,  I  will  point  out  what  aspects 


we  focused  on  to  construct  our  typical  sequelae  repre¬ 
sentations. 
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ANALYTICAL  CONSIDERATIONS 


[Slide  31  The  period  of  acute  illness  is  convention¬ 
ally  divided  into  the  prodromal  or  initial  phase  (1  to  3 
days  after  exposure)  and  manifest  illness  phase  (1  to  6 
weeks  after  exposure);  these  times,  of  course,  depend  upon 
dose  level. 

The  main  prodromal  symptoms  relevant  to  functional 
impairment  are  nausea,  vomiting,  anorexia,  diarrhea,  fluid 
loss,  and  electrolyte  imbalance.  Concomitant  effects,  either 
a  direct  result  of  radiation  or  secondary  to  fluid  loss,  are 
headaches,  fainting,  and  prostration.  Other  early  effects 
with  a  different  pathophysiological  base  are  fatigue  (more 
specifically  fatigability)  and  weakness. 

The  manifest-illness  phase  is  dominated  by  bleeding, 
fever,  infection,  and  ulceration  due  to  injury  of  the 
hemopoietic  system.  Higher  doses  can  produce  hypotension, 
dizziness,  and  disorientation.  Fluid  loss,  electrolyte 
imbalance,  and  delayed  diarrhea  recur  after  relatively  high 
doses  because  of  damage  to  the  intestines. 

Symptoms  were  also  grouped  according  to  physiological 
relationship  designated  by  the  double  capital  letters  in 


parentheses,  i.e.,  upper  gastrointestinal  distress  (UG) , 
lower  gastrointestinal  distress  (LG) ,  fatigability  and 
weakness  (FW) ,  etc. 

[Slide  4]  To  express  symptom  severity,  we  use  the 
terms  "mild,"  "moderate,"  and  "severe" — common  clinical 
distinctions — but  we  are  unable  to  quantify  the  terms  for 
most  of  the  symptoms.  For  example,  by  what  measure  could  one 
distinguish  moderate  from  severe  nausea  or  headaches?  Some 
quantitative  links  are  possible  but  we  mainly  rely  on 
qualitative  understanding  of  the  distinctions. 

Mild  vomiting  or  diarrhea  may  mean  a  single  to  a  few 
episodes;  moderate,  several  episodes;  and  severe,  many  and 
profuse  episodes.  Fatigue  and  weakness  are  potentially 
quantifiable  since  exertion,  necessary  to  reveal  those 
symptoms,  is  measurable;  however,  few  data  have  been  col¬ 
lected.  With  regard  to  hypotension,  "mild"  may  refer  to  a  10 
percent  drop  in  blood  pressure;  and  "severe,"  a  drop  of  30 
percent  or  more.  We  quantify  hemopoietic  injury  by  drops  in 
lymphocyte,  granulocyte,  and  platelet  counts.  Net  continued 
fluid  losses  of  up  to  2  liters  are  considered  mild  to 
moderate;  more  than  2  liters  can  be  severe. 
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ACUTE  RADIATION  SYMPTOMS 


Nausea  i  upper  Gastrointestinal  Distress  (UG) 
Vomiting  ) 

•  Diarrhea  -  Lower  Gastrointestinal  Distress  (LG) 


Fatigability  t  (FW) 
Weakness  ’ 


•  Hypotension  ) 
Dizziness  >  (HY) 
Disorientation  ' 


•  Infection  ) 

Bleeding  ( (IB) 

Fever  ( 

Ulceration  ) 

•  Fluid  loss 

Electrolyte  imbalance  ) 
Headache  V  (PL) 

Fainting  I 

Prostration  ' 
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Also,  our  correlations  of  severity  level  with  time  are 


gross  representations;  even  though  severity  may  vary  con¬ 
siderably  over  time,  we  did  not  attempt  to  develop  detailed 
time- intensity  profiles  because  of  the  lack  of  appropriate 
quantitative  data. 


[Slide  5]  For  symptom  incidence  we  largely  utilized 
results  from  probit  analysis  of  data  from  the  accounts  of 
accident  victims  and  over  2000  therapy  patients.  However, 
the  probit  data  are  not  specifically  correlated  with  postex¬ 
posure  time,  but  we  were  able  to  infer  the  time  dependence 
of  symptom  incidence  from  various  portions  of  the  body  of 
literature  reviewed. 

Results  of  probit  analysis  relating  symptom  incidence 
to  radiation  dose  are  shown  on  lognormal  probability  paper; 
the  straight-line  curves  assume  a  lognormal  distribution  and 
the  concave  curves,  a  normal  distribution.  The  lognormal 
form  generally  fits  the  data  better.  Exceptions  are  the 
incidence  of  diarrhea  at  doses  above  250  rads  and  the  in¬ 
cidence  of  fatigue  at  doses  above  300  rads,  for  which  we  use 
the  normal  curve  estimates. 
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SYMPTOM  SEVERITY 
QUALITATIVE  REPRESENTATION 


Mild 

Mild  to  moderate 
Moderate 


Moderate  to  severe 


Severe 


The  probit  curves  for  diarrhea  reflect  data  gathered 
over  a  six-week  period  following  exposure.  Since  the  data 
are  not  time-resolved,  the  curves  include  later  as  well  as 
early  occurrences.  However,  clinical  experience  indicates 
that  delayed  diarrhea  occurs  a  few  days  to  a  week  after 
exposure.  Early  diarrhea  most  likely  doesn't  occur  unless 
the  dose  is  at  least  300  rads  where  about  10  percent  would 
experience  one  or  two  episodes  3  to  6  hours  after  exposure; 
at  2000  to  3000  rads,  30  percent  could  have  early  diarrhea. 

[Slide  6]  Because  of  the  sparseness  and  quality  of 
human  data,  quantification  of  incidence  of  death  due  to 
radiation  is  uncertain.  The  three  curves  shown  here  il¬ 
lustrate  the  uncertainty:  the  LD5Q  ranges  from  256  to  325 
rads.  We  use  the  right-hand  curve  estimated  for  healthy 
adult  humans  based  on  an  LD5q  of  450  rads  (cGy)  air  dose; 
the  slope  is  from  experiments  with  large  animals,  since 
lethality  dose  curves  tend  to  be  parallel  for  dogs,  swine, 
sheep,  goats,  etc. 


[Slide  7]  For  mean  time  of  death,  we  use  the  ranges 
given  by  the  heavily  outlined  boxes  to  indicate  the  period 
over  which  fatalities  are  likely  to  occur  for  a  given  dose 
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DOSE,  RADS  (cGy) 
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level.  Data  are  insufficient  to  specify  a  time  distribution 
of  mortality.  These  ranges  are  consistent  with  a  consensus 
of  fairly  recent  opinion  and  some  more  recent  data  points 
for  human  fatalities.  For  comparison  a  curve  previously 
drawn  through  three  earlier  human  data  points  is  shown;  the 
curve's  shape  is  inferred  from  animal  data.  Relative  to  the 
curve,  a  trend  toward  earlier  death  with  doses  above  around 
1300  rads  is  suggested. 

[Slide  81  Estimates  of  the  onset  of  prodromal  symptoms 
are  based  on  the  data  points  logarithmically  plotted  here 
against  dose  indicated  by  symbols  representing  the  four 
categories.  Accident,  therapy,  and  expert  opinion  symbols 
indicate  firm  information;  lines  connecting  the  composite 
symbols  indicate  fairly  wide  ranges  of  uncertainty;  arrows 
indicate  open-ended  values  and  thus  reflect  an  even  larger 
measure  of  uncertainty. 

Because  the  data  vary  greatly  in  density  and  precision, 
numerical  curve  fitting  techniques  (such  as  regression 
analysis)  were  not  applied.  Instead,  we  drew  a  curve  through 
the  data  points  to  represent  "typical”  individuals  favoring 
accident  and  therapy  data  over  the  other  two  categories, 
which  shows  onset  time  inversely  proportional  to  dose. 
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(The  exact  trend  at  the  high  end  of  the  dose  range  is 
uncertain  because  of  the  lack  of  empirical  data.  However, 
the  curve  is  supported  by  the  accident  data  point  at  8800 
rads,  and  expert  opinion  suggesting  that  persons  exposed  to 
several  thousand  rads  will  probably  show  the  eucire  range  of 
prodromal  symptoms  within  5  to  15  minutes.) 

Prodromal  symptom  duration,  manifest  illness  phase 
onset  and  duration  were  related  to  dose  plotting  appropriate 
data  in  a  similar  manner.  I  should  like  to  point  out  that 
our  estimates  of  duration  account  for  symptom  occurrence  in 
three  ways:  continuously  at  one  level  of  severity  (e.g., 
anorexia);  continuously  with  varying  severity  (e.g., 
fatigue);  and  intermittently  (e.g.,  vomiting  and  diarrhea). 

Let  me  briefly  comment  on  the  variety  of  dose  rate 
contained  in  the  data.  Accident  victims  were  exposed  to  from 
hundreds  to  thousands  of  rads  in  a  fraction  of  a  second  as 
compared  to  therapy  patients  undergoing  total-body  irradia¬ 
tion  exposed  to  from  a  few  to  30  rads/min  given  over  periods 
of  minutes  to  hours. 

Given  the  great  differences  in  dose  rate  between  acci¬ 
dent  victims  and  therapy  patients,  we  might  expect  their 
prodromal  symptoms  to  begin  at  different  times.  However, 
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correlation  of  onset  time  with  dose  level  shows  no  marked 


difference  between  the  two  groups  suggesting  that  onset 
depends  more  on  total  dose  rather  than  dose  rate  in  the 
therapy  range.  We,  therefore,  believe  data  from  therapy 
patients  are  applicable  for  deriving  typical  acute  response 
patterns  for  symptoms  such  as  nausea,  vomiting,  and 
fatigability. 


SYMPTOMATOLOGY  DESCRIPTIONS 


[Slide  9)  We  subdivided  the  dose  range  of  50  to  3000 
rads  into  eight  subranges  reflecting  important  pathophsyical 
events.  Doses  in  the  lowest  range,  50  to  100  rads,  cause 
minor  acute  damage  to  the  hemopoietic  system  and  mild 
prodromal  effects  (nausea,  vomiting,  anorexia)  in  a  small 
number  of  irradiated  persons. 

In  the  dose  range  100  to  200  rads,  prodromal  effects 
and  injury  to  the  hemopoietic  system  (primarily  the  bone 
marrow  stem  and  precursor  cells)  increase  significantly. 
However,  victims  will  probably  survive,  except  for  the  2  to 
5  percent  who  will  die  after  doses  of  about  200  rads.  The 
probability  of  death  increases  if  victims  are  already 
weakened  by  other  conditions,  such  as  infection. 

Although  survival  is  possible  from  200  to  350  rads, 
prodromal  effects  become  pronounced.  Victims  also  suffer 
moderate  to  severe  bone  marrow  damage.  As  the  dose  reaches 
about  325  rads,  50  percent  who  do  not  receive  appropriate 
medical  care  may  die  within  60  days. 


Prom  350  to  550  rads,  symptoms  are  more  severe  affect¬ 
ing  nearly  all  exposed.  If  untreated,  50  to  99  percent  may 


K 


die,  primarily  because  of  extensive  injury  to  the 
hemopoietic  system,  manifested  by  overwhelming  infections 
and  bleeding. 

Between  550  and  750  rads,  responses  begin  to  reflect 
the  combined  effects  of  gastrointestinal  and  hemopoietic 
damage.  Survival  is  almost  impossible  short  of  a  compatible 
bone  marrow  transplant  and  extensive  medical  care.  Nearly 
everyone  irradiated  at  this  level  suffers  severe  prodromal 
effects  during  the  first  day  after  exposure. 

Between  750  and  1000  rads,  injuries  are  much  more 
severe  due  to  greater  depletion  of  bone  marrow  stem  cells, 
increased  gastrointestinal  damage,  and  systemic  complica¬ 
tions  from  bacterial  endotoxins  entering  the  blood  stream. 

At  1000  to  2000  rads,  death  results  in  less  than  two 
weeks  from  septecemia  due  to  severe  gastrointestinal  injury, 
complicated  by  complete  bone  marrow  damage  and  the  cessation 
of  granulocyte  production.  Above  about  1300  rads,  death  may 
occur  sooner  from  electrolyte  imbalance  and  dehydration  due 
to  vomiting  and  diarrhea,  especially  in  hot  or  humid  condi¬ 
tions.  Extremely  severe  gastrointestinal  and  cardiovascular 
damage  causes  death  within  2  to  5  days  after  doses  of  2000 
to  3000  rads. 
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DOSE  RANGES  AND 

ASSOCIATED  PATHOPHYSIOLOGICAL  EVENTS 


Pathophysiological  Events 

Dose  Range 
rads(cGy) 

Prodomal 

Effects 

Manifest-Illness  Effects 

Survival 

50*100 

Mild 

Slight  decrease  in  blood 
cell  count 

Virtually  certain 

100-200 

Mild  to 
moderate 

Beginning  symptoms  of 
bone  marrow  damage 

Probable  (>  90  percent) 

200-350 

Moderate 

Moderate  to  severe  bone 
marrow  damage 

Possible- 

Bottom  third  of  range: 
LD5/60 

Middle  third:  LD10/60 
Top  third:  LD50/6Q 

350-550 

Severe 

Severe  bone  marrow 
damage 

Death  within  3%-6  weeks- 
Bottom  half:  LD90/6O 
Top  half:  LDgg/60 

550-750 

Severe 

Bone  marrow  pancytopenia 
and  moderate  intestinal 
damage 

Death  within  2-3  weeks 

750-1000 

Severe 

Combined  gastrointestinal 
and  bone  marrow  damage; 
hypotension 

Death  within  1  -2%  weeks 

1000-2000 

Severe  gastrointestinal  damage 

Upper  half  of  range;  early  transient 
incapacitation;  gastrointestinal  death 

Death  within  5-12  days 

2000-3000 

Gastrointestinal  and  cardiovascular  damage 

Death  within  2-5  days 
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[Slide  101  Here  is  an  example  of  the  symptom  charts  we 
developed  for  the  350  to  500  rads  dose  range.  Symptom  charts 
for  all  eight  dose  ranges  are  included  in  the  handouts. 

The  charts  portray  symptom  onset  severity,  duration, 
and  incidence.  Symptom  incidence  is  expressed  as  percentage 
affected.  Both  symptom  incidence  and  severity  are  expressed 
broadly  with  time  because  of  the  lack  of  specific  time- 
resolved  data. 

(Essentially,  this  chart  shows  nearly  all  persons 
exposed  experience  moderate  to  severe  symptoms  the  first  day 
postexposure,  moderating  within  a  day  or  so.  Fluid  loss 
without  replacement  from  severe  or  prolonged  vomiting  can 
affect  electrolyte  balance.  With  the  exception  of 
fatigability  and  weakness,  there  is  a  period  of  symptom 
remission  from  about  2  days  to  a  week  following  exposure.) 

[Slide  11]  Next,  I  will  present  some  summary  slides 
which  illustrate  the  overall  picture  we  formed  for  acute 
radiation  sickness  symptomatology;  these  illustrations  are 
in  the  handouts  so  I  won't  spend  much  time  discussing  them. 


146 


I 


n 

13$ 

Mi 

V. 

>» 


1 


,V’< 


Y  *>  J 

iV1 


■»> 


SJ 


O’l 

A'* 

«,iy 

AV 

Yi‘i 

’A'f 


SYMPTOMS  FOR  DOSE  RANGE  350  TO  550  RAOS  (cGy) 


Po»c«xpo»urt  Tlao 


12  3  4  5  6  7 


0  4  3  12  16  20  24 


12  14  5  6 


*S«v«r«  drop  la  placolocs:  Proa  1  <  105/«aJ  to  0.1  «  l05-0/aaJ. 
^Svrart  drop  la  sranulocycoo:  Proa  6  «  l0}/na^  to  0.5  *  l03-0/«aJ. 


“Sovoro  drop  la  lyaphocycoa: 


1  «  lOJ/ni 


CO  0.4-0. 1  *  to-’/* 
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This  is  a  summary  chart  of  the  six  symptom  groups  for 
all  eight  dose  ranges  for  times  up  to  six  weeks  following 
exposure.  Each  of  the  six  groups  include  symptoms  indicated 
in  Slide  3;  symptom  incidence  is  also  indicated  along  with 
severity,  denoted  by  shading. 


[Slide  121  This  chart  shows  the  combined  symptom 
course  of  acute  radiation  sickness  for  the  eight  dose 
ranges.  Notice  that  beginning  with  the  350  to  550  rads  dose 
range,  fatigability  and  weakness  (FW)  continues  with  no 
remission  period. 


[Slide  131  Here  are  some  surface  plots  which  il¬ 
lustrate  symptom  severity  related  to  dose  level  and  postex¬ 
posure  time.  Because  dose  and  time  are  expressed  logarithmi¬ 
cally  (to  include  the  full  range  of  dose  and  time)  visual 
distortion  is  introduced  where  the  slope  of  the  surface  for 
small  doses  and/or  early  times  appears  much  more  pronounced 
than  that  of  large  doses  and/or  times. 

(Slide  14]  These  are  additional  surface  plots  of 
symptom  severity  plus  a  composite  representation  of  all  the 
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ACUTE  RADIATION  SYMPTOMS 
PostMpoiur*  Tim* 


Do** 

R*cta  Symptom 


Days 

1  2  3  4  5  6  7 


-•Gy)  Group  0  4  8  12  16  20  24 


1  2  3  4  5  8 


symptom  groups  referred  to  as  "Acute  Radiation  Sickness." 


In  conclusion,  we  point  out  that  general  application  of 
these  descriptions  carries  some  reservations  because  of 
differences  in  population  characteristics,  environmental 
conditions,  and  medical  attention.  Accident  victims  and 
therapy  patients  had  the  benefit  of  medical  care,  which  may 
not  be  available.  On  the  other  hand,  factors  such  as  youth, 
vigor,  and  motivation  may  to  some  extent  counter  the  effects 
of  radiation  injury  on  functional  capability.  The  data  do 
not  permit  a  quantitative  assessment  of  the  tradeoffs  be¬ 
tween  postexposure  medical  care  and  preexposure  robustness. 
Although  both  are  important,  we  believe  that  preexposure 
health  condition  may  be  relatively  less  important  in  resist¬ 
ing  the  debilitating  effects  of  radiation  than  postexposure 
medical  care,  barring  prior  infection  (bacterial  or  viral) 
or  serious  disease.  The  medical  care  (environmental  control, 
broad  spectrum  antibiotics,  antifungal,  antiviral,  fresh 
platelets,  total  parenteral  feeding,  electrolyte  replace¬ 
ment)  given  to  the  Chernobyl  accident  victims  seems  to  bear 
this  out  based  on  larger  LD5Q  estimates  (400  to  600  rads) 
made  by  some.  However,  complications  from  secondary  injury 
(thermal  and  beta  radiation  skin  burns)  as  well  as  dose 
protraction  (few  to  several  hours)  compounds  such  an  assess- 


ACUTE  RADIATION  SICKNESS 
DOSE-TIME  SYMPTOM  SEVERITY 
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[Slide  15]  Our  study  of  acute  radiation  sickness 


provides  the  basis  for  estimating  the  expected  functional 


capability  of  individuals  after  a  brief  exposure  to  ionizing 


radiation.  Here  we  can  get  an  idea  of  the  gross  functional 


impairment  of  individuals  who  may  be  involved  in  emergency 


civil  defense  activities.  Physically  undemanding  tasks  are 


those  that  require  cognitive  capability  such  as  communica¬ 


tion,  navigation  and  plotting,  counting,  fine  motor  coor¬ 


dination,  etc.  "Effective"  may  correspond  to  from  normal  to 


75%  capability;  "degraded",  from  25  to  75%  capability;  and 


'ineffective",  less  than  25%  capability.  The  low  portion  of 


the  curves  from  10-15  hrs  postexposure  corresponds  to  the 


maximum  expression  of  initial  symptoms  of  acute  radiation 


sickness.  The  high  portion  of  the  curves  from  2  to  4  days 


represents  a  transient  period  of  remission  of  symptoms 


followed  by  a  second  onset  due  to  the  latent  expression  of 


symptoms  caused  by  radiation  damage  to  the  blood-forming 


organs. 


[Slide  16]  Physically  demanding  tasks,  for  example, 


lifting  and  moving  heavy  objects,  climbing  flights  of 


stairs,  running  or  walking  appreciable  distances,  etc.,  are 


more  adversely  affected  by  radiation  than  unotmanding  ones. 


For  example,  functional  capability  for  performing  physically 
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TIME  AFTER  EXPOSURE  (min) 


demanding  tasks  would  continue  to  be  affected  after  several 
hours  postexposure  for  doses  as  low  as  150  rads  (cGy) . 


ACUTE  RADIATION  SYMPTOM  CHARTS 


DISCUSSION 


Doses  of  50  to  100  Rads  (cGy) 

Acute  radiation  effects  at  this  level  are  mild  and 
occur  only  during  the  first  day  after  exposure.  Blood  cell 
counts  may  drop  slightly,  but  typical  victims  will  surely 
survive. 

Doses  of  100  to  200  Rads  (cGy) 

Between  20  and  70  percent  of  exposed  persons  develop 
nausea  and  vomiting.  About  30  to  60  percent  complain  of 
fatigue  and  weakness.  Significant  destruction  of  bone  marrow 
stem  cells  may  lead  to  a  25  to  35  percent  drop  in  blood  cell 
production.  As  a  result,  mild  bleeding,  fever,  and  infection 
may  occur  during  the  fourth  and  fifth  postexposure  weeks.  Up 
to  5  percent  may  die  5  to  6  weeks  after  exposure  to  200  rads 
(cGy)  . 

Doses  of  200  to  350  Rads  (cGy) 

Commencing  in  this  dose  range,  about  10  percent  may 
experience  one  or  two  episodes  of  moderate  diarrhea  4  to  6 
hours  postexposure.  Most  victims  tire  easily  and  experience 
mild  to  moderate  weakness  intermittently  over  the  6  weeks. 
Under  normal  conditions,  vomiting  and  diarrhea  are  not 
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enough  to  cause  serious  fluid  loss  and  electrolyte  im¬ 
balance.  In  hot  or  humid  conditions,  however,  combined  fluid 
loss  and  electrolyte  imbalance  could  become  serious. 

Injury  to  the  hemopoietic  system  is  indicated  by 
moderate  bleeding,  fever,  infection,  and  ulceration  3  to  5 
weeks  postexposure;  more  than  50  percent  of  those  exposed 
are  affected.  During  the  fourth  and  fifth  weeks,  moderate 
diarrhea  may  complicate  their  condition.  Five  to  50  percent 
of  nontreated  persons  may  die  during  the  fifth  week;  death 
comes  earlier  to  those  with  preexisting  infections,  for 
example,  of  the  upper  respiratory  tract. 

Doses  of  350  to  550  Rads  (cGy) 

In  this  next  higher  dose  range,  symptoms  are  more 
severe  affecting  nearly  all  exposed  persons.  Severe  and 
prolonged  vomiting  can  take  a  toll  on  electrolyte  balance 
which  may  be  accelerated  by  perspiration  through  heat, 
humidity,  or  activity. 

Nearly  all  show  moderate  to  severe  fatigue  and  weakness 
for  many  weeks.  If  untreated,  50  to  99  percent  may  die, 
primarily  because  of  extensive  injury  to  the  hemopoietic 
system,  manifested  in  overwhelming  infections  and  bleeding 
during  the  third  to  sixth  weeks.  Nausea,  vomiting,  and 
anorexia  may  recur  at  that  time.  Diarrhea,  electrolyte 
imbalance,  and  headaches  affect  at  least  half.  The  condition 
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of  the  lethally  irradiated  during  their  last  days  may  be 
complicated  by  dizziness,  disorientation,  fainting,  prostra¬ 
tion,  and  symptoms  of  infection  and  bleeding. 

Doses  of  550  to  750  Rads  (cGy) 

Virtually  all  exposed  persons  experience  severe  nausea 
and  vomiting  the  first  postexposure  day,  moderating  over  the 
next  day  or  two.  During  that  time  they  also  experience 
dizziness  and  disorientation. 

With  the  near-maximum  destruction  of  bone  marrow  stem 
cells  and  absence  of  granulocytes,  untreated  persons  lose 
their  defense  against  infection.  By  the  end  of  the  first 
postexposure  week,  infection  is  rampant  from  endogenous 
bacteria  that  have  escaped  from  the  injured  gastrointestinal 
tract. 

The  combination  of  hemopoietic  damage  and  gastrointes¬ 
tinal  lesions  reduces  the  survival  of  all  untreated  persons 
to  2  to  3  weeks.  During  the  entire  time  they  suffer  from 
severe  fatigue  and  weakness.  Toward  the  end  of  the  first 
week,  nausea,  vomiting,  and  anorexia  recur.  Moderate  to 
severe  diarrhea  may  begin  as  early  as  the  fourth  day.  Severe 
bleeding,  headaches,  hypotension,  dehydration,  electrolyte 
imbalance,  and  fainting  complicate  the  condition  of  all 
prior  to  death. 
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The  survival  time  for  untreated  persons  diminishes  to  2 
to  2  1/2  weeks.  Symptoms  resemble  those  experienced  at  the 
preceding  dose  range,  with  the  following  notable  dif¬ 
ferences: 

o  Severe  nausea  and  vomiting  may  continue  into  the 
third  day  before  moderating. 

o  During  the  first  day,  hypotension  affects  about  80 
percent;  moderate  fever,  30  to  45  percent. 

o  Electrolyte  imbalance  is  a  persistent  problem  from 
the  sixth  hour  on. 

o  All  have  moderate  to  severe  headaches  during  the 
first  day. 

o  Nearly  three-quarters  are  prostrate  before  the  end 
of  the  first  week. 

Much  of  our  description  of  symptoms  at  this  dose  range 
derives  from  postexposure  observations  of  patients  treated 
with  total-body  irradiation  for  leukemia.  There  are  undeni¬ 
able  difficulties  in  extrapolating  from  sick  people  under 
close  medical  attention  to  otherwise  healthy  individuals. 
Nevertheless,  therapy  patients  constitute  the  only  substan¬ 
tial  number  of  irradiated  persons  whose  reactions  have  been 
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Severe  nausea  and  vomiting  affect  all  within  30  minutes 
of  exposure  and  continue  intermittently,  along  with 
anorexia,  until  death  the  second  week.  Severe  headaches 
begin  after  about  4  hours  and  continue  for  2  to  3  days. 
Symptom  severity  may  diminish  somewhat  during  days  3  to  5. 
Gastrointestinal  injury  predominates,  manifested  4  to  6  days 
after  exposure  by  the  abrupt  return  of  severe  nausea,  vomit¬ 
ing,  anorexia,  and  diarrhea,  along  with  high  fever,  ab¬ 
dominal  distension,  and  undetectable  peristalsis  (ileus). 
During  the  second  week,  severe  dehydration,  hemoconcent ra¬ 
tion,  and  circulatory  collapse,  compounded  by  septicemia, 
lead  to  coma  and  death. 

Doses  of  2000  to  3000  Rads  (cGy) 

Symptoms  are  more  severe  versions  of  those  described 
for  the  preceding  dose  range.  Gastrointestinal  injury 
predominates,  complicated  by  cardiovascular  lesions. 
Prodromal  effects,  including  severe  headache  and  drowsiness, 
appear  almost  immediately  after  exposure  and  may  persist  as 
the  gastrointestinal  syndrome  develops.  Severe  dehydration 
and  electrolyte  imbalance  are  manifested  several  hours  after 
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Symptoms  for  dose  range  50  to  100  rads  (cGy) 


Postexposure  Time 


12  3  4  5  6  7 


Symptom 


Vomiting 

(retching) 

Anorexia 

Diarrhea 

(cramps) 


Fatigue 

Weakness 


Hypotension 
Dlzz i ness 
Disorientation 


Bleeding 


Infection 

Ulceration 


Fluid  loss/electro- 
lyte  imbalance 

Headache 
Fainting 
Prostrat ion 


0  4  8  12  16  20  24 


1  2  3  4  5  6 


References  for  this  group  of  symptoms:  l,  7,  15,  22,  26-30,  )3,  37,  42,  50,  62, 
66,  71.  72,  76,  78,  80,  85,  92,  96,  97,  105-107,  HO.  111.  These  symptoms  not  ob¬ 
served  in  American  servicemen  exposed  co  approximately  78  rads  (>C.v)  of  lallout 
radiation,  according  to  Refs.  26-29. 

^  n  c - ~  r»  rn..n  ..  f  ~  •  [  7,  14,  15,  26“)0,  33,  50,  64,  72, 


References  for  this  group  of  symptoms: 
101,  105,  106,  108. 


^Slight  drop  in  lymphocyte,  platelet,  and  granulocyte  counts;  no  overt  symptoms. 
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Symptoms  for  dose  range  100  to  200  rads  (cGy) 
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Bleed lngJ 


Infection 

Ulceration 


Fluid  loss/electro 
lyte  Imbalance 

Headache 

Fainting 

Prostration 


References  for  this  group  of  sympcoms:  4,  6,  7,  15,  17,  21,  22,  26-29,  33,  37, 
42,  44,  45,  50.  51.  56,  62,  65,  66,  74,  76,  78-82,  85-87,  90-92,  95,  96,  104,  105, 
107,  111. 

Ren  percent  of  the  Marshallese  victims  exposed  to  175  rads  (cCy)  experienced 
diarrhea  during  the  first  postexposure  day,  according  co  Ref.  4. 

References  foe  this  group  of  symptoms:  7,  50,  60,  65,  81,  85,  86,  90,  101,  102. 

References  for  this  group  of  symptoms:  6,  7,  14,  18,  21,  25,  26-29,  33,  35,  60, 
62,  64,  65.  75,  76,  78,  79.  81,  82,  85,  89,  101,  104,  105,  107,  110. 

Slight  to  moderate  drop  In  platelets:  from  3  *  lOVmml  to  1.8-0. 8  x  105/mral. 
^Slight  to  moderate  drop  in  granulocytes:  from  6  *  loVmm^  co  4. 5-2.0  *  10^/mm^. 
^Slight  to  moderate  drop  in  lymphocytes:  from  3  *  loRmm^  co  2.0-1.0  *  10^/mm^. 
References  for  this  event:  4,  14,  61,  81. 
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Symptoms  for  dose  range  200  to  350  rads  (cGy) 


References  for  this  group  of  symptoms:  4,  6,  7,  12,  15,  18,  21,  22,  32,  33,  35, 
37.  41-44,  47,  50,  51,  53,  56,  60,  62,  65,  75,  76,  80-82,  85-88,  90.  92.  95-97, 
105-111. 

References  for  this  group  of  symptoms:  1,  7,  14,  43,  65,  71,  76,  81,  101,  104, 
109,  110. 

References  for  this  group  of  symptoms:  2,  6,  7,  14,  15,  33,  35,  39,  54,  56,  60, 
62,  64,  65,  71-,  75,  76,  79,  81,  82,  85,  88,  95,  97,  101,  105.  106,  110. 

^Moderate  drop  In  platelets:  from  3  x  loVmml  to  0.8-0. 1  *  lO^/mm^. 

"Moderate  drop  in  granulocytes:  from  6  *  10^/mm^  to  2. 0-0. 5  *  10^/mra-*. 

Moderate  to  severe  drop  in  lymphocytes:  from  3  *  10^/mm^  to  1.0-0. 4  *  101/  mra^. 
^Epilation. 

References  for  this  event:  4,  7,  14,  61,  81. 
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Symptoms  for  dose  range  350  to  550  rads  (cGy) 

I 


Postexposure  Tin* 

Hours  Days  Weeks 

12  3  4  5  6  7 

Symptom  0  4  8  12  16  20  24  1  2  3  4  5  6 


^References  for  this  group  of  symptoms:  l,  4,  6,  7,  14,  15,  21.  22.  32,  33,  37, 
41,  42.  50.  51,  53.  56.  58.  62,  76,  77,  79-82.  85,  87,  89.  90,  95-97.  101,  105.  106, 
U0,  111. 


^References  for  this  group  of  sympcoms:  1,  6,  7,  14,  47,  51,  53,  65,  78.  81.  85, 

90,  101. 

C 

References  for  this  group  of  symptoms:  77,  89. 

^References  for  this  group  of  symptoms:  l,  6,  7,  14,  15,  18.  31.  33,  35,  54,  58, 

62,  64,  65,  67,  71,  75,  76.  77,  81,  85.  90,  95,  101,  105-107.  110,  111. 

^Severe  drop  in  platelets:  from  3  *  10^/ora^  to  0.1  *  lO^-0/mra^. 
f Severe  drop  In  granulocytes:  from  6  *  loVma^  to  0.5  *  10^-0/mra^. 

^Severe  drop  in  lymphocytes:  from  3  *  10^/mm^  to  0.4-0. 1  *  10^/mm^. 
hEpilation. 

^References  for  this  group  of  symptoms:  7,  14,  71,  81,  89,  101. 

^Mlld  intestinal  damage. 

Le 

References  for  this  event:  4,  14,  61,  81. 
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Symptoms  for  dose  range  750  to  1000  rads  (cGy) 


Poscexposure  Time 


Hours  Days  Weeks 

1  2  3  4  5  6  7 

Symptom  0  4  8  12  16  20  24  1  2  3  4  5  6 


References 

for 

this 

group  of 

symptoms: 

6. 

7. 

14,  15,  21.  22.  33,  41, 

,  44,  47 

50, 

►  51.  53,  57, 

.  62, 

,  65. 

71.  76, 

79-82.  85. 

87. 

90, 

,  93,  95,  96.  101,  106. 

85, 

References 
,  101. 

for 

this 

group  of 

symptomr • 

1.  8. 

7,  14,  32.  47,  51,  53. 

65.  67, 

References 

for 

this 

group  of 

symptoms: 

6. 

7, 

14,  15,  33,  50,  61,  62, 

,  65,  71 

76,  79,  80,  82,  85,  87,  89,  93,  101,  105.  106. 

Rlood  pressure  drops  25  percent;  temperature  increases  to  102*F,  according  to 


Ref.  33. 


References  for  this  group  of  symptoms:  6,  7,  14,  15,  21,  30,  31,  33,  51,  53, 
61.  62.  66,  71.  72,  79-*2,  85,  87,  89.  90,  93,  95,  101.  106. 

Platelet  count  drops  to  zero. 

*Crsnulocyte  count  dropa  to  zero. 

^Lymphocyte  count  drops  to  zero. 

1. 


Ep list  Ion. 

References  for  this  group  of  symptoms: 
72,  76,  79-82,  85,  89.  101,  106. 

Sloderace  to  severe  intest inei  damage. 
1. 


6,  7,  14,  15.  33.  50.  61.  62.  65,  71. 


References  for  this  event:  4,  14,  61.  77,  81. 
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Symptoms  for  dose  range  1000  to  2000  rads  (cGy) 


Poscexposure  Time 


Symptom 


Hours  Days 

1  2  i  4  5 

8  12  16  20  ZU 

n  !  I  l" 


Weeks 


T 


Nausea 

Vomiting 

(retching) 

Anorexia 

Diarrhea 
(cramps ) 


Katigue^ 

Weakness 


Hypotens ionc 

Dizziness 

Disorientation 


Bleeding^ 

Fever 

Infection 

Ulceration 


Fluid  loss/electro- 
iyte  imbalance*1 

Headache 

Fainting 

Piostracion 


DeathJ 


»l00t  severe* 
■100Z  severe* 


-4- 


100Z  severe 


-100Z- 


20Z 

severe 


•  100Z  severe  • 

'  I00Z  severe • 


^100Z 

— 100Z* 
severe 


-100Z  moderate' 

to  severe 

•  I00Z  severe  — 

>L00Z  severe - 


__45-80Z  — 
moderate  to 
severe 


(e 

Kf)— II 


-(«)- 


100Z  severe 
100Z  severe 


100Z  severe 


»100Z  moderate 

to  severe 

■— 100Z  severe 


-J _ L. 


(D 


I"  ■"  ,  I 


100Z  severe 

100Z  severe 
80Z  severe 
80Z  severe 


-<k) - IP—  100Z- 


Referencee  for  this  group  of  symptoms:  6, 

62,  65,  67,  71,  76,  79-82,  85,  87,  90,  95,  101. 

References  for  this  group  of  symptoms:  6, 

101. 

References  for  this  group  of  symptoms:  6, 

76,  79-82.  85,  87,  89,  101,  105.  106. 

References  for  this  group  of  symptoms:  6, 

65.  67,  71.  76.  79-82,  85,  87,  90,  95.  101.  105 

*Platelet  count  drops  to  zero. 

^Granulocyte  count  drops  to  zero. 

'Lymphocyte  count  drops  to  zero. 

Referehees  •for  this  group  of  symptoms:  7, 
79-82,  85,  95.  101,  105,  106,  108. 

Revere  Intestinal  damage. 

References  for  this  svent:  4,  14,  61.  81.  98. 
Renal  failure,  according  to  Ref.  98. 


14. 

15, 

21, 

22. 

31, 

33. 

47, 

50. 

105. 

106. 

\  14, 

47, 

53. 

65. 

71. 

78, 

85, 

90, 

14. 

15. 

33, 

50. 

61. 

62, 

65. 

71. 

\  14, 

15, 

21, 

22, 

33. 

47. 

50, 

62, 

106. 

4,  15 

.  33. 

.  50, 

,  61 

.  65, 

,  68. 

.  75. 

.  76 

171 


*%* 

'V 

$ 

ft 


'4 

'& 

( 

k 

& 

G 

I 

*4 

H1 

$ 

i 

•1*5 

0 

♦:* 


Symptoms  for  dose  range  2000  to  3000  rads  (cGy) 


Poscexposure  Time 

Hours  Days 

1  2  3  4  5  6  7 

8  12  16  20  24  1 


3  4  5  6 


Symptom 


Vomiting 

(retching) 


Diarrhea 

(cramps) 


Fatigueb 

Weakness 


Hypotension** 

Dizziness 

Disorientation 


Blee<lingd 

Fever 

Infection 

Ulceration 


Fluid  loss/electro¬ 
lyte  imbalance* 

Headache 

Fainting 

Prostration 


“References  for  this  group  of  symptoms:  6,  7,  14,  15,  21,  22,  42.  44,  50,  62, 
65,  71,  76,  79-82,  85,  87.  90,  95,  101.  106,  108. 

References  for  this  group  of  symptoms:  6,  7,  14,  15,  21,  22,  50.  62,  65,  71, 
76,  79-82,  85,  87,  90,  94,  95,  101,  105,  106,  108. 

References  for  this  group  of  symptoms:  14,  50,  61,  71,  76,  79,  80.  82,  85,  87, 
94,  95,  101,  105,  108. 

References  for  this  group  of  symptoms:  7,  14,  50,  61,  71,  75,  76,  79-82,  85, 
87,  94,  95,  101,  105,  106,  108. 

Crsnulocyte  counc  drops  to  zero. 

^Lymphocyte  count  drops  to  zero. 

References  for  this  group  of  symptoms:  14,  15,  33,  50,  61,  62,  71,  75,  76,  79, 
80,  82,  85,  87,  95,  101,  105,  108. 

References  for  this  event:  4,  14,  61,  77,  81. 


172 


REFERENCES 


Adams,  R.,  and  S.  Cullen,  The  Final  Epidemic;  Physicians 
and  Scientists  on  Nuclear  War,  Educational  Foundation  for 
Nuclear  Science,  Chicago,  Illinois,  1981. 

Adelstein,  S.  J.,  and  J.  B.  Dealy,  Jr.,  "Hematologic 
Responses  to  Human  Whole-Body  Irradiation,"  Am.  J. 
Roentgenol.,  Radium  Ther.  and  Nucl.  Med.,  Vol.  93, 

1965,  pp.  927-934. 

Ainsworth,  E.  J.,  G.  F.  Leong ,  K.  Kendall,  and  E.  L.  Alpen, 
"Cooperative  Lethality  Response  of  Neutron-  and  X-Irradiated 
Dogs,"  Radiat.  Res.,  Vol.  26,  1965,  pp.  33-34. 

Alpen,  E.  L. ,  Radiological  Hazard  Evaluation:  A  Critical 
Review  of  Present  Concepts  and  a  New  Approach  Thereto,  U.S. 
Naval  Radiological  Defense  Laboratory,  San  Francisco,  Calif¬ 
ornia,  Report  USNRDL-TR-186,  1957. 

Anderson,  A.  C.,  and  L.  S.  Rosenblatt,  "The  Effect  of  Whole 
Body  X-Irradiation  on  the  Median  Life  Span  of  Female  Dogs 
(Beagles) ,"  Radiat.  Res.,  Vol.  39,  1969,  pp.  177-200. 

Andrews,  G.  A.,  "Total-Body  Irradiation  in  the  Human  Being," 
Excerpta  Medica,  International  Congress  Series  105,  1965, 
pp.  1583-1589. 

Andrews,  G.  A.,  and  R.  J.  Cloutier,  "Accidental  Acute  Radiation 
Injury,"  Arch.  Environ.  Health,  Vol.  10,  1965,  pp.  498-507. 

Anno ,  G .  H . ,  Acute  Radiation  Response  in  Humans:  Informal 
Comments  by  Physicians  and  Radiobiologists,  Pacific-Sierra 
Research  Corporation,  Note  492  (rev.),  June  1983. 

Anno,  G.  H.,  H.  L.  Brode,  and  R.  Washton-Brown ,  Initial  Human 
Response  to  Nuclear  Radiation,  Pacific-Sierra  Research 
Corporation,  Note  477,  April  1982. 

Barret,  A.,  "Total  Body  Irradiation  (TBI)  before  Bone  Marrow 
Transplantation  in  Leukemia:  A  Cooperative  Study  from  the 
European  Group  for  Bone  Marrow  Transplantation,"  Br .  J. 
Radiol. ,  Vol.  55,  1982,  pp.  562-567. 

Baum,  S.  J.,  "Erythrocyte  Stem  Cell  Kinetics  in  the  Post¬ 
irradiation  Rat,"  Radiat.  Res.,  Vol.  30,  1967,  pp.  316-324. 


Blakely,  J.,  The  Care  of  Radiation  Casualties,  Charles  C. 
Thomas  Co.,  Springfield,  Illinois,  1968. 

Bond,  V.  P.,  "Radiation  Mortality  in  Different  Mammalian 

Species,"  in  V.  P.  Bond  adn  T.  Sugahara  (eds.).  Comparative 
Cellular  and  Species  Radiosensitivity,  Igaku  Shoin,  Tokyo, 
1969. 

Bond,  V.  P.,  T.  M.  Fliedner,  and  J.  0.  Archambeau,  Mammalian 
Radiation  Lethality:  A  Disturbance  in  Cellular  Kinetics, 
Academic  Press,  New  York,  1965 . 

Bond,  V.  P.,  T.  M.  Fliedner,  and  E.  P.  Cronkite,  "Evaluation 
and  Management  of  Heavily  Irradiated  Individuals,"  J.  Nucl. 
Med. ,  Vol .  1,  1960,  pp.  221-238. 

Borison,  H. ,  and  S.  C.  Wang,  "Physiology  and  Pharmacology  of 
Vomiting,"  Pharmacol.  Rev.,  Vol.  5,  1953,  pp.  193-230. 

Brown,  W.  M.  Court,  "Symptomatic  Disturbance  after  Single- 
Therapeutic  Dose  of  X-Rays,"  Br .  Med.  J.,  11  April  1953, 
pp.  802-805. 


Brown,  W.  M.  Court,  and  J.  D.  Abbatt,  "The  Effect  of  a  Single 
Dose  of  X-Rays  on  the  Peripheral  Blood  Count  of  Man," 

Br.  J.  Haematol.,  Vol.  1,  1955,  pp.  75-85. 

Brown,  W.  M.  Court,  and  R.  Doll,  Leukemia  and  Aplastic  Anemia 
in  Patients  Irradiated  for  Ankylosing  Spondylitis,  British 
Medical  Research  Council,  special  report  series  1-50,  Her 
Majesty's  Stationery  Office,  London,  1957. 

- ,  "Mortality  from  Cancer  and  Other  Causes  after  Radio¬ 
therapy  for  Ankylosing  Spondylitis,"  Br.  Med.  J.,  Vol.  2, 
1965,  pp.  1327-1332. 

Brucer,  M.  B.  (comp.),  The  Acute  Radiation  Syndrome:  A 
Medical  Report  on  the  Y-12  Accident,  June  16,  1958,  U.S. 
Atomic  Energy  Commission,  Report  ORINS-25,  April  1959 . 

Cairnie,  A.  B.,  and  H.  A.  Robitaille,  Arguments  for  the 
Greater  Importance  of  the  Prodromal  Syndrome  Than 
Incapacitation  (Involving  Early  Transient  Incapacitation) 
the  Consideration  of  Radiation.Ef fects  in  Irradiated 
Military  Personnel,  Together  with  a  Proposal  to  Simulate 
the  Prodromal  Effects  Using  Lithium  Carbonate,  [Canadian] 
Defense  Research  Establishment,  Ottawa,  Report  836, 

December  19BQ. 


Carpenter,  D.  0.,  Early  Transient  Incapacitation:  A 
Review  with  Consideration  of  Underlying  Mechanisms, 

Armed  Forces  Radiobiological  Research  Institute, 
Bethesda,  Maryland,  Scientific  report  SR-79-1,  1979. 

Chin  H.,  and  S.  C.  Wang,  "Locus  of  Emetic  Action  Following 
Irradiation,"  Proc.  Soc.  Exp.  Biol.  Med..  Vol.  85,  1954, 
pp.  472-474. 

Conard,  R.  A.,  "Acute  Myologenous  Leukemias  Following 
Fallout  Radiation  Exposure,"  J.  Am.  Med.  Assn.,  Vol. 

232,  1975,  pp.  1356-1357. 

- ,  A  Twenty-Year  Review  of  Medical  Findings  in  a 

Marshallese  Population  Accidentally  Exposed  to  Radio¬ 
active  Fallout,  Brookhaven  National  Laboratory,  Long 
Island,  New  York,  Report  BNL  50424  (TLD-4500) ,  1975. 

Conard,  R.  A.,  et  al.,  Medical  Survey  of  Rongelap  People 
Eight  Years  after  Exposure  to  Fallout,  Brookhaven 
National  Laboratory,  Long  Island,  New  York,  Report  BNL 
780  (T-296) ,  1963. 

Conard,  R.  A.,  et  al..  Review  of  Medical  Findings  in  a 
Marshallese  Population  Twenty-six  Years  after  Accidental 
Exposure  to  Radioactive  Fallout,  Brookhaven  National 
Laboratory,  Long  Island,  New  York,  Report  BNL  51261 
(TLD-4500) ,  1980, 

Cronkite,  E.  P.,  et  al.,  "Response  of  Human  Beings  Acci¬ 
dentally  Exposed  to  Significant  Fallout  Radiation,"  J. 
Am.  Med.  Assn..  Vol.  159,  1955,  pp.  430-434. 

Cronkite,  E.  P.,  V,  P.  Bond,  and  C.  L.  Dunhan,  Some 
Effects  of  Ionizing  Radiation  on  Human  Beings,  Atomic 
Energy  Coraiission  Report  TID  5358,  1956. 

Cronkite,  E.  P.,  and  V.  P.  Bond,  "Acute  Radiation  Syndrome 
in  Man,"  U.S.  Armed  Forces  Med.  J.,  Vol.  9,  1958, 
pp.  313-324. 

Cronkite,  E.  P.,  et  al.,  "Diagnosis  and  Therapy  of  Acute 
Radiation  Injury,"  Ch.  2  of  Atomic  Medicine,  3d  ed.. 

The  Williams  &  Williams  Co.,  Baltimore,  Maryland,  1959. 

Cronkite,  E.  P.,  and  V.  P.  Bond,  Radiation  Injury  in  Man, 
Charles  C.  Thomas  Co.,  Springfield,  Illinois,  1960. 


175 


34.  Cronkite,  E.  P., 


m 


"The  Effects  of  Dose,  Dose  Rate  and  Depth 
Dose  upon  Radiation  Mortality,"  Proceedings  of  a  Sym¬ 
posium  on  the  Control  of  Exposures  of  the  Public  to 
Ionizing  Radiation  in  the  Event  of  Accident  or  Attack, 
National  Council  on  Radiation  Protection  and  Measure¬ 
ment,  Bethesda,  Maryland,  May  1982,  pp.  21-27. 


35.  Dienstbier,  Z.,  M.  Arient,  and  J.  Pospisil,  "Hamatologische 

Veranderung  bei  der  Strahlenkrankheit-IV;  Hamokoagulations- 
veranderungen,"  Atompr axis ,  Vol.  9,  .1963,  pp.  189-194. 

36.  Edsall,  D.  L. ,  and  R.  Pemberton,  "The  Nature  of  the  General 

Toxic  Reaction  Following  Exposure  to  X-Rays,"  Am.  J.  Med. 
Sc. ,  Vol.  133,  1970,  pp.  426-431. 

37.  Ellinger,  F.,  et  al. ,  "A  Clinical  Study  of  Radiation  Sick¬ 

ness,"  Am.  J.  Roentgenol.,  Radium  Ther.  and  Nucl.  Med., 
Vol.  68,  1952,  pp.  275-280. 


$ 
to 

D 

:»5| 

r 


Fanger,  H.,  and  C.  C.  Lushbaugh,  "Radiation  Death  from 
Cardiovascular  Shock  Following  a  Criticality  Accident," 
Arch.  Path.,  Vol.  83,  May  1967,  pp.  446-460. 

Finch,  S.  C.,  "Recognition  of  Radiation-Induced  Late  Bone 
Marrow  Changes,"  Ann.  N.  Y.  Acad.  Sci.,  Vol.  145,  1967, 
pp.  748-754. 

Fliedner,  T.  M. ,  and  D.  van  Beckum,  private  communication, 
July  1983. 

Gerstner,  H.  B.,  "Acute  Clinical  Effects  of  Penetrating 
Nuclear  Radiation,"  J.  Am.  Med.  Assn.,  Vol.  168, 

27  September  1958,  pp.  381-388. 


- ,  "Acute  Radiation  Syndrome  in  Man,"  U.S.  Armed 

Forces  Med.  J.,  Vol.  9,  1958,  p.  313. 


- ,  "Reaction  to  Short  Term  Radiation  in  Man,"  Annual 

Rev.  Med.,  Vol.  11,  1960,  pp.  289-302. 


- ,  "Practical  Implication  of  the  Initial  Reaction  to 

Penetrating  Ionizing  Radiation,"  unpublished  manuscript, 
U.S.  Air  Force  School  of  Aerospace  Medicine,  Brooks 
Air  Force  Base,  Texas,  1970. 


L 


•J® 


Gilberti,  M.  V.,  "The  1967  Radiation  Accident  near 
Pittsburgh,  Pennsylvania,  and  a  Follow-Up  Report," 
in  K.  F.  Hubner  and  S.  A.  Fry  (eds.),  The  Medical 
Basis  for  Radiation  Accident  Preparedness,  Elsevier 
North  Holland,  Inc.,  New  York,  1980. 


176 


Glasstone,  S.,  and  P.  J.  iolan,  The  Effects  of  Nuclear 
Weapons,  3d  ed.,  U.S.  Departments  of  Defense  and 
Energy,  1977. 

Grant,  G.  A.,  et  al.,  A  Predictive  Study  of  the  Incidence 
of  Vomiting  in  Irradiated  Military  Personnel.  [Canadian] 


Defense  Research  Establishment,  Ottawa,  Report  817,  19 

Hall,  E.  J.,  Radiobiology  for  the  Radiobiologist,  Harper 
&  Row,  Hagerstown,  Maryland,  1978. 

Hempleraann,  L.  H.,  C.  C.  Lushbaugh,  and  G.  L.  Voetz, 
"What  Happened  to  the  Survivors  of  the  Early  Los 
Alamos  Nuclear  Accidents?"  in  K.  F.  Hubner  and 
S.  A.  Fry  (eds.).  The  Medical  Basis  for  Radiation 
Accident  Preparedness,  Elsevier  North  Holland,  Inc., 
New  York,  1980. 

Howland,  J.  W. ,  "Injury  and  Recovery  from  Ionizing  Radia 
tion  Exposure,"  Annual  Rev.  Med.,  Vol.  7,  1956,  pp. 
225-244. 


Howland,  J.  W. ,  M.  Ingram,  H.  Mermagen,  and  C.  L. 

Hansen,  Jr.,  "The  Lockport  Incident:  Accidental  Par¬ 
tial  Body  Exposure  of  Humans  to  Large  Doses  of  X- 
Irradiation,"  Diagnosis  and  Treatment  of  Acute  Radia¬ 
tion  Injury,  International  Atomic  Energy  Agency  and 
World  Health  Organization,  International  Documents  Ser 
vice.  New  York,  1961,  pp.  11-26. 

Hubner,  K.  F.,  and  S.  A.  Fry  (eds.).  The  Medical  Basis 
for  Radiation  Accident  Preparedness,  Elsevier  North 
Holland,  Inc.,  New  York,  1980 . 

Ingram,  M.,  J.  W.  Howland,  and  C.  L.  Hansen,  "Sequential 
Manifestation  of  Acute  Radiation  Injury  vs.  'Acute 
Radiation  Syndrome'  Stereotype,"  Ann.  N  Y  Acad.  Sci., 
Vol.  114,  1964,  pp.  356-367. 

Ishida,  M.,  and  I.  Matsubayashi ,  An  Analysis  of  Early 
Mortality  Rates  Following  the  Atomic  Bomb  in  Hiroshima, 


Atomic  Bomb  Casualty  Committee,  Hiroshima,  Technical 
report  20-61,  1948. 

Ishimaru,  T.,  et  al.,  "Leukemia  in  Atomic  Bomb  Survivors 
Hiroshima  and  Nagasaki,  1  October,  1950-30  September, 
1966,"  Radiat.  Res.,  Vol.  45,  1971,  pp.  216-233. 


M 

h 

i 

f 


m 

i*Ci 

o 


m 

i 


I 

p 

1 

| 

n 

M" 

i*Jl 

»s 

•ij 

ishg 

a 

1 

1 


Vi 

I 


C 


56.  Jammet,  H.  P.,  "Treatment  of  Victims  of  the  Zero-Energy 

Reactor  Accident  at  Vinca,"  Diagnosis  and  Treatment 
of  Acute  Radiation  Injury,  International  Atomic 
Energy  Agency  and  World  Health  Organization,  Inter¬ 
national  Documents  Service,  New  York,  1961. 

57.  Jammet,  H.  P.,  et  al.,  "Etude  de  Six  Cas  d' Irradiation 

Totale  Aigue  Accidentale , "  Rev.  Franc.  Etud.  Clin. 

Biol. ,  Vol.  4,  1959,  pp.  210-225. 

58.  Jammet,  H.  P.,  R.  Gongora,  R.  Le  Go,  and  M.  J.  Doley, 

"Clinical  and  Biological  Comparison  of  Two  Acute 
Accidental  Irradiations:  Mol  (1965)  and  Brescia 
(1975),"  in  K.  F.  Hubner  and  S.  A.  Fry  (eds.),  The 
Medical  Basis  for  Radiation  Accident  Preparedness, 
Elsevier  North  Holland,  Inc.,  New  York,  1980. 

59.  Karas,  J.  S.  and  J.  B.  Stanbury,  "Fatal  Radiation  Syn¬ 

dromes  from  an  Accidental  Nuclear  Excursion,"  New  Eng. 

J.  Med.,  Vol.  272,  No.  15,  1965,  pp.  755-761. 

60.  Kumatori,  T. ,  "Hematological  Effects  on  Heavily  Irradiated 

Japanese  Fishermen,"  in  T.  Sugahara  and  O.  Hug  (eds.), 
Biological  Aspects  of  Radiation  Protection,  Igaku 
Shoin,  Tokyo,  1971. 

61.  Langham,  W.  H.  (ed.),  Padiobiological  Factors  in  Manned 

Space  Flight,  National  Academy  of  Sciences,  National 
Research  Council  Publication  1487,  1967. 

62.  Laumets,  E.,  Time  History  of  Biological  Response  to 

Ionizing  Radiation,  U.S.  Naval  Radiobiological  Defense 
Laboratory,  San  Francisco,  California,  Report  USNRDL- 
TR-905 ,  November  1965. 

63.  Levin,  S.,  et  al.,  Early  Biological  Effects  from  Initial 

Nuclear  Radiation  in  Hiroshima  and  Nagasaki,  Armed 
Forces  Radiobiological  Research  Institute,  Bethesda, 
Maryland,  1983  (forthcoming) . 

64.  Liebow,  A.  A.,  S.  Warren,  and  E.  De  Coursey,  "Pathology 

of  Atomic  Bomb  Casualties,"  Am.  J.  Path.,  Vol.  25, 

1949,  p.  853. 

65.  Lushbaugh,  C.  C.,  "What  Can  We  Expect  to  Happen?"  Rocky 

Mount.  Med.  J.,  January  1962. 


178 


- ,  Recent  Progress  in  Assessment  of  Human  Resistance 

to  Total-Body  Irradiation,  National  Academy  of  Sciences, 
National  Research  Council  conference  paper  671135, 

April  1968. 

- ,  "Reflections  on  Some  Recent  Progress  in  Human  Radio- 

biology,"  Advances  in  Radiation  Biology,  Vol.  3, 

Academic  Press,  New  York,  1969,  pp.  277-315. 

- ,  "Theoretical  and  Practical  Aspects  of  Models 

Explaining  ' Gastroinest inal  Death'  and  other  Lethal 
Radiation  Syndromes,"  in  V.  P.  Bond  and  T.  Sugahara 
( eds . ) ,  Comparative  Cellular  and  Species  Radiosen¬ 
sitivity,  Igaku  Shoin,  Tokyo,  1969,  pp.  288-297. 

- ,  "Human  Radiation  Tolerance,"  Ch.  10  in  J.  Parker,  Jr 

and  V.  R.  West  (eds.),  Bioastronautics  Data  Book, 

National  Aeronautics  and  Space  Administration,  report 
NASA-S-30G0S,  1973. 

- ,  "The  Impact  of  Estimates  of  Human  Radiation  Toler¬ 
ance  upon  Radiation  Emergency  Management,"  Proceedings 
of  a  Symposium  on  the  Control  of  Exposure  of  the  Public 
to  Ionizing  Radiation  in  the  Event  of  Accident  or 
Attack,  National  Council  on  Radiation  Protection and 
Measurement,  Bethesda,  Maryland,  May  1982,  pp.  46-57. 

Lushbaugh,  C.  C.,  et  al.,  "Clinical  Studies  of  Radiation 
Effects  in  Man,"  Radiat.  Res,  Suppl.,  Vol.  1,  1967, 
pp.  398-412. 

Martin,  E.  J.,  and  R.  H.  Rowland,  Castle  Series.  1954, 
Defense  Nuclear  Agency,  Report  DNA  6035F,  April  1982. 

Martinez,  R.  G.,  et  al.,  "Observations  on  the  Accidental 
Exposure  of  a  Family  to  a  Source  of  Cobalt-60,"  Rev. 

Med.  Inst.  Mex.,  Seguro  Social,  Vol.  3,  1964,  pp. 

14-68  (translated  by  F.  V.  Comas) . 

McCandless,  J.  B.,  "Accidental  Acute  Whole-Body  Gamma 
Irradiation  of  Seven  Clinically  Well  Persons,"  J.  Am. 

Med.  Assn..  Vol.  192,  1965,  pp.  85-88. 

McFarland,  W. ,  and  H.  A.  Pearson,  "Hematological  Events 
as  Dosimeters  in  Human  Total-Body  Irradiation," 

Radiol. ,  Vol.  80,  1963,  pp.  850-855. 

McLean,  A.  S.,  "Early  Adverse  Effects  of  Radiation," 

Br.  Med.  Bull.,  Vol.  29,  1973,  pp.  69-73. 


Messerschmidt ,  0.,  Medical  Procedures  in  a  Nuclear 
Disaster ,  Verlag  Karl  Thieming,  Munich,  1979. 


Miller,  L.  S.,  G.  H.  Fletcher,  and  H.  B.  Gerstner, 
"Radiobiologic  Observations  on  Cancer  Patients 
Treated  with  Whole-Body  X-Radiation,"  Radiat.  Res., 

Vol .  8,  1958,  pp.  150-165. 

NATO  Handbook  on  the  Medical  Aspects  of  NBC  Defensive 
Operations,  U.S.  Departments  of  the  Army,  Navy,  and 
Air  Force,  Report  AMED  P-6,  August  1973. 

Nuclear  Weapons  Employment  Doctrine  and  Procedures, 

U.S.  Army,  FM  101-31-1,  1977. 

Ohkita,  T.  II,  "A  Review  of  Thirty  Years  Study  of  Hiroshima 
and  Nagasaki  ATomic  Bomb  Survivors,"  Jpn.  J.  Radiat. 

Res. ,  supplement  16,  1975,  pp.  49-66. 

Oughterson,  A.  W. ,  and  S.  Warren,  Medical  Effects  of  the 
Atomic  Bomb  in  Japan,  McGraw-Hill  Book  Company,  New 
York,  1956. 

Pendic,  B.,  "The  Zero-Energy  Reactor  Accident  at  Vinca," 
Diagnosis  and  Treatment  of  Acute  Radiation  Injury, 
International  Atomic  Energy  Agency  and  World  Health 
Organization,  International  Documents  Service,  New 
York,  1961. 

Porvaznik,  M.,  "Tight  Junction  Disruption  and  Recovery 
after  Sublethal  Gamma  Irradiation,"  Radiat.  Res., 

Vol.  78,  1979,  pp.  233-259. 

Prasad,  K.  N.,  Human  Radi ation  Biology,  Harper  &  Row, 
Hagerstown,  Maryland,  1974. 

Radiological  Factors  Affecting  Decision-Making  in  a 
Nuclear  Attack,  National  Council  on  Radiation  Pro¬ 
tection  and  Measurement,  B  thesda,  Maryland,  Report 
42,  November  1974. 

R  &  D  Associates,  "Collateral  Damage  Implications  of 
Low  Radiation  Dose  Criteria  for  Battlefield  Nuclear 
Operations,"  unpublished  manuscript,  January  1980. 


Rider,  W.  D.,  and  R.  Hasselback,  "The  Symptomatic  and 
Hematological  Disturbance  Following  Total  Body  Radia¬ 
tion  at  300-rad  Gamma-Ray  Irradiation,"  lectures  pre¬ 
sented  at  McGill  University,  Montreal,  August  1967, 
in  Guidelines  to  Radiological  Health,  U.S.  Public 
Health  Service,  Washington,  D.C.,  1968,  pp.  139-144. 

Rubin,  P.,  and  G.  W.  Casarett,  Clinical  Radiation  Path¬ 
ology  ,  W.  B.  Saunders  Company,  Philadelphia,  1968. 

Saenger,  E.  L.  (ed.).  Medical  Aspects  of  Radiation 
Accidents,  U.S.  Atomic  Energy  Commission,  1963. 

Saenger,  E.  L.,  et  al..  Metabolic  Changes  in  Humans 
Following  Total  Body  Irradiation,  Defense  Atomic 
Support  Agency,  Washington,  D.C.,  Report  1633,  1964. 

Saenger,  E.  L.,  et  al..  Radiation  Effects  in  Man:  Mani¬ 
festations  and  Therapeutic  Efforts,  Defense  Nuclear 
Agency,  Report  2751,  October  1971. 

Salazar,  0.  M. ,  et  al.,  "Systemic  (Half-Body)  Radiation 
Therapy:  Response  and  Toxicity,"  J.  Radiat.  Oncol. 

Biol.  Phys.,  Vol.  4,  1978,  pp.  937-950. 

Shipman,  T.  L. ,  "Acute  Radiation  Death  Resulting  from 
an  Accidental  Nuclear  Critical  Excursion,"  J.  Occup. 
Med. ,  Vol.  3,  No.  3,  1961,  pp.  145-192. 

Storb,  R. ,  "Total-Body  Irradiation  and  Marrow  Trans¬ 
plantation,"  Transplant.  P roc.,  Vol.  9,  1977, 
pp.  1113-1119. 

Thoma,  George  E.,  Jr.,  and  N.  Wald,  "The  Diagnosis  and 
Management  of  Accidental  Radiation  Injury,"  J.  Occup. 
Med. ,  Vol.  1,  August  1959,  pp.  421-447. 

- ,  "The  Acute  Radiation  Syndrome  in  Man,"  in  Epi¬ 
demiology  of  Radiation  Injury,  postgraduate  course 
syllabus,  St.  Louis  University,  School  of  Medicine, 
August  1961. 

Thomas,  E.  D.,  et  al.,  "Allogeneic  Marrow  Grafting  for 
Hematologic  Malignancy  Using  HL-A  Matched  Donor- 
Recipient  Sibling  Pairs,"  Blood ,  Vol.  38,  No.  3, 
September  1971,  pp.  267-287. 

Thomas,  E.  D.  K.  Dicke,  and  G.  Santos,  private  communi¬ 
cation,  July  1983. 


18] 


100.  Thomson,  J.  F.,  et  al. ,  "Studies  on  the  Effects  of  Con¬ 

tinuous  Exposure  of  Animals  to  Gamma  Radiation  from 
Cobalt-60  Plane  Sources,"  Am.  J.  Roentqenol.,  Vol. 

69,  1953,  pp.  830-838. 

101.  Upton,  A.  C.,  "Effects  of  Radiation  on  Man,"  Annual  Rev. 

Nucl.  Sci..  Vol.  18,  1968,  pp.  495-528. 

102.  - ,  Radiation  Injury:  Effects,  Principles,  and 

Per spect ives ,  The  University  of  Chicago  Press, 

Chicago,  1969. 

103.  Vodopick,  H.,  and  G.  A.  Andrews,  The  Clinical  Effects 

of  an  Accidental  Radiation  Exposure,  Oak  Ridge 
Associated  Universities  Publications,  Oak  Ridge, 
Tennessee,  July  1973. 

104.  - ,  "The  University  of  Tennessee  Comparative  Animal 

Research  Laboratory  Accident  in  1971,"  in  K.  F.  Hubner 
and  S.  A.  Fry  (eds.).  The  Medical  Basis  for  Radiation 
Accident  Preparedness,  Elsevier  North  Holland,  Inc., 

New  York,  1980. 

105.  Wald,  N.,  and  G.  E.  Thoma,  Jr.,  Radiation  Accidents; 

Medical  Aspects  of  Neutron  and  Gamma-Rav  Exposures, 

Oak  Ridge  National  Laboratory,  Report  OnRL-2748,  Part 
B,  March  1961. 

106.  Wald,  N.,  G.  E.  Thoma,  and  G.  Broun,  "Hematologic  Manifes¬ 

tations  of  Radiation  Exposure  in  Man,"  Proq.  Hematol., 
Vol.  3,  1962,  pp.  1-5. 

107.  Warren,  S.,  "The  Early  Changes  Caused  by  Radiation,  J.  Mt . 

Sinai  Hospital,  Vol.  19,  1952,  pp.  443-455. 

108.  - ,  "You,  Your  Patient,  and  Radioactive  Fallout,"  New 

Eng.  J.  Med..  Vol.  266,  1962,  pp.  1123-1125. 

109.  Yochmowitz,  M.  G. ,  and  G.  C.  Brown,  "Performance  in  a 

12-Hour,  300-Rad  Profile,"  Aviat.  Space  Environ.  Med., 
1977,  pp.  241-247. 

110.  Zellmer,  R.  W.,  "Human  Ability  to  Perform  after  Acute 

Sublethal  Radiation,"  Mil.  Med.,  Vol.  126,  September 
1961,  pp.  681-687. 

111.  Zellmer,  R.  W.,  and  J.  E.  Pickering,  Biological  Effects 

of  Nuclear  Radiation  in  Primates,  U.S.  Air  Force 
School  of  Aviation  Medicine,  Brooks  Air  Force  Base, 
Texas,  Technical  Report  60-77,  1960. 


182 


SECTION  2 

GLOBAL  CLIMATE  SIMULATIONS 


18$ 


A  COMPARISON  OF 

EULERIAN  AND  LAGRANGIAN  METHODS 
FOR  TRACER  TRANSPORT  IN  A  GCM 


Bob  Malone,  Gary  Glatzmaier, 
and  David  Langley 


Earth  and  Space  Sciences  Division 
Los  Alamos  National  Laboratory 


A  Comparison  of  Eulerian  and  Lagrangian  Methods 
for  Tracer  Transport  in  a  GCM 

Robert  C.  Malone  and  Gary  A.  Glatzmaier 
Los  Alamos  National  Laboratory 


Recent  studies  of  "nuclear  winter"  with  three-dimensional  models  have 
highlighted  the  need  for  numerical  algorithms  for  the  interactive 
transport  of  smoke  that  have  high  spatial  resolution  but  do  not  suffer  from 
excessive  numerical  diffusion.  Two  issues  of  importance  for  nuclear 
winter  are  (1)  the  residence  time  of  smoke  in  a  high-altitude,  stably 
stratified  layer  heated  by  sunlight  and  (2)  the  degree  of  nonuniformity 
("patchiness")  of  areal  coverage  by  smoke.  Numerical  diffusion,  which 
particularly  afflicts  coarse  resolution  Eulerian  transport  methods,  may 
lead  to  significant  errors  in  simulation  studies  of  these  effects. 

Although  no  observations  exist  to  compare  with  nuclear  winter 
simulations  in  which  the  atmospheric  circulation  and  structure  are 
modified  by  solar  heating  of  smoke,  analogues  in  the  normal  atmosphere  do 
exist  that  can  be  used  as  benchmarks  for  large-scale  transport 
simulations  with  a  GCM.  Two  for  which  good  data  are  available  are  (1 )  the 
residence  time  in  the  stratosphere  of  radioactive  particulates  and  gases 
from  high-yield  atmospheric  nuclear  weapon  tests  and  (2)  the  dispersal  in 
the  troposphere  of  deuterated-methane  tracer.  Simulations  of  both 
phenomena  with  Eulerian  and  Lagrangian  transport  models  will  be 
presented  and  compared  with  observational  data  to  illustrate  the  relative 
merits  of  the  methods. 

Simulations  of  interactive  smoke  transport  with  solar  heating  of 
smoke  will  also  be  presented  using  both  methods. 


s 

!« 
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What’s  New  in  the  Model 


Absorption  of  infrared  radiation  by  smoke 

Stability-dependent  vertical  mixing  in  the  planetary 
boundary  layer  and  free  atmosphere 

Subsurface  storage  and  vertical  diffusion  of  heat  and 
water  in  soil 

Annual  cycle  of  solar  declination 

Diurnal  cycle  of  solar  zenith  angle 

High-resolution,  non-diffusive  Lagrangian  tracer 
transport 
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EULER1AN  TRANSPORT  MODEL 


Calculates  tracer  concentration  at  discrete  grid  points 
Has  same  resolution  as  underlying  GCM 
Limitation:  number  of  grid  points  affordable  in  GCM 


LAGRANGIAN  TRACER  MODEL 


Tracer  is  represented  by  a  large  number  of  discrete 
particles 

Particle  coordinates  not  restricted  to  GCM  grid  points 

Limitation:  large  number  of  particles  needed  to 
adequately  represent  local  concentrations  on  grid 


HEAVY  METHANE  EXPERIMENTS 

•  Three  releases  of  1  kg  CD4  in  midtroposphere  near  55S, 
165E  (January,  June,  October  1984) 

•  Precisely  characterized  source:  time,  place,  altitude, 
amount.  Inert  tracer. 

•  Samples  collected 

-  on  airplane  flights  between  Antarctic  stations,  and 
to/from  New  Zealand  (1  hour  averages) 

-  at  surface  stations  (3  day  averages) 

-  for  60  days  after  release 

•  Detection  limit  - 10“17  g  CD4  /  g  air 

•  Samples  above  limit  found  out  to  20  days 
(longest  data  series  for  a  tropospheric  tracer?) 

•  Principal  shortcoming:  limited  sampling  precludes 
complete  picture  of  large-scale  distribution. 
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MASS  MIXING  RATIO 


COMPARISON  OF  OBSERVATIONS  WITH  LAGRANGIAN 
AND  EULERIAN  TRANSPORT  MODELS 
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STRATOSPHERIC  RESIDENCE  TIME 

Some  disagreement  among  observations,  but  recent 
data  indicates  1 /e-times  of  1  year  in  lower  stratosphere 

Observations  show  comparable  residence  times  in 
stratosphere  for  gases  (HTO),  radionuclides,  and 
volcanic  sulfates 

Radionuclide  data  suggests  no  removal  for  a  long  period 
(3-6  months)  following  lower  stratospheric  injections 

Eulerian  transport  model  gives  immediate  fallout  with 
1/e-times  in  lower  stratosphere  of  2-3  months 

Lagrangian  transport  model  gives  negligible  removal  for 
2-3  months  and  slow  removal  thereafter  (longer  run  with 
annual  cycle  needed  to  determine  residence  time) 
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Changes  in  Smoke  Transport 


Smoke  lofting  faster  in  new  model  than  in  old  model 

-  buoyancy  greater  in  hot  spots  (stability-dependent 
vertical  mixing  of  heat  less  efficient  than  DCA) 

-  more  inhomogeneous  smoke  distribution 
(Lagrangian  model)  drives  smaller-scale  circulations 
than  does  old  (Eulerian)  model 

-  diurnal  cycle:  greater  heating  in  daytime,  followed  by 
dilution  due  to  dispersion;  then  subsidence  at  night? 

Less  perturbation  of  tropical  troposphere  (less 
numerical  diffusion  horizontally  and  vertically) 

Longer  residence  times  (=  1  yr)  than  before  (=  1/2  yr), 
evolving  plausibly  with  changing  solar  declination 


Smoke  Mass  Fraction  Remainin 

180  Tg  injections,  NAS  profile 

Lagrangian  transport  of  smoke 
Annual  cycle 
Precipitation  scavenging 
Gravitational  sedimentation  (1  \±) 

No  coagulation 
No  reactions  with  ozone 


Smoke  mass  fraction  after  20  days 


Injected  in 

Laaranaian  model 

Eulerian  model 

late  June 

47% 

42% 

mid-April 

42% 

mid-August 

42% 

mid-October 

32% 

late  December 

27% 

25% 

Smoke  mass  fraction  after  200  davs 

Injected  in 

Laaranaian  model 

-  , 

{ 

i 

late  June 

I 

31%  (annual  cycle  to  January) 

Chemical  Reaction  of  Smoke  with  Ozone 

( Conjectures  about  likely  interactions  with  transport) 

If  reaction  time  is  short  (<  few  days) 

-  thermal  perturbations  by  smoke  reduced 

-  smoke  lofting  and  spreading  reduced 

-  ozone  depletion  concentrated  in  northern 
midlatitudes 

If  reaction  time  is  long  (>  1  month) 

-  thermal  perturbation  and  smoke  distribution  similar 
to  present  model 

-  ozone  depletion  will  occur  in  both  hemispheres 

-  solar  heating  of  smoke  further  reduces  ozone 
through  temperature-dependence  of  ozone 
reactions  (Vupputuri) 


Improvements  in  the 
Los  Alamos  Global  Climate  Model 
and  their  Effects  on 
Nuclear  Winter  Simulations 


Gary  Glatzmaier,  Bob  Malone, 
and  David  Langley 


Earth  and  Space  Sciences  Division 
Los  Alamos  National  Laboratory 
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IMPROVEMENTS  IN  THE  LOS  ALAMOS  GLOBAL  CLIMATE  MODEL 
AND  THEIR  EFFECTS  ON  NUCLEAR  WINTER  SIMULATIONS 


Gary  A.  Glatzmaier  and  Robert  C.  Malone 
Los  Alamos  National  Laboratory 

We  have  improved  the  Los  Alamos  global  climate  model  in  several  ways 
in  order  to  obtain  more  realistic  nuclear  winter  simulations.  Our  new 
model  has  absorption  of  infrared  radiation  by  smoke,  stability-dependent 
vertical  diffusion  in  the  atmosphere,  heat  and  moisture  storage  in  the  soil, 
Lagrangian  transport  of  smoke,  and  diurnal  and  annual  cycles. 

The  absorption  of  infrared  radiation  by  smoke  moderates  surface 
cooling  by  about  15-25%,  in  agreement  with  previous  investigations. 

The  fixed,  linear,  stability-independent  vertical  diffusion  coefficients 
for  momentum,  heat,  and  water  vapor  and  the  infinitely  efficient 
"dry-convective-adjustment"  algorithm  of  the  old  model  have  been 
replaced  with  time-varying  nonlinear  vertical  diffusion  coefficients  that 
depend  on  the  local  atmospheric  stability  and  wind  shear.  For  the 
simulated  normal  atmosphere,  this  nonlinear  diffusion  is  significant  only 
in  the  lower  one  to  two  kilometers,  the  planetary  boundary  layer,  where 
the  atmosphere  can  become  unstable  to  small-scale  turbulence.  However, 
for  a  simulated  nuclear  winter,  solar-heating  of  lofted  smoke  produces  a 
stable  inversion  in  the  troposphere  which  severely  reduces  the  vertical 
diffusion  coefficients  during  the  first  week. 

For  the  calculation  of  surface  temperatures,  we  replaced  the  diagnostic 
"heat-balance"  algorithm  of  the  old  model  with  prognostic  equations 
describing  heat  storage  and  thermal  diffusion  in  six  soil  layers  down  to 
five  meters  below  the  surface.  In  addition,  the  time-independent 
saturated  soil  moisture  condition  of  the  old  model  was  replaced  with 
prognostic  equations  describing  surface  and  subsurface  storage  of 
moisture. 

With  the  nonlinear,  time  and  spacially  dependent,  vertical  diffusion  in 
the  atmosphere  and  the  heat  storage  in  the  soil,  we  are  now  able  to  run 
with  a  diurnal  cycle. 

The  Eulerian  transport  of  smoke  in  our  old  model  has  been  replaced  by  a 
Lagrangian  transport  scheme  (described  by  R.  Malone)  which  allows  us  to 
make  longer  nuclear  winter  simulations  with  an  annual  cycle. 

The  influences  of  all  these  model  improvements  on  surface  cooling  and 
smoke  lofting  and  residence  time  will  be  described. 


Improvements  in  the 
Los  Alamos  Global  Climate  Model 


•  Subsurface  storage  and  vertical  diffusion  of 

moisture  in  a  2-level  soil  model  (10cm,  50cm). 

•  Subsurface  storage  and  vertical  diffusion  of 

heat  in  a  6-level  soil  model  (5cm  -  5m). 

•  Stability-dependent  vertical  diffusion  in  the 

20-level  atmospheric  model  (30m  -  30km). 

•  Absorption  of  infrared  radiation  by  smoke. 

•  High-resolution,  non-diffusive  Lagrangian 

transport  of  smoke. 

•  Annual  cycle  of  solar  declination. 

•  Diurnal  cycle  of  solar  zenith  angle. 
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Richardson  number-dependent  vertical  diffusion 
coefficient  computed  locally 

Replaces  "dry  convective  adjustment"  and 
neutral-stability  linear  vertical  diffusion 

In  normal  atmosphere,  active  only  in  PBL 

During  nuclear  winter,  active  in  PBL  arid  in  free 
atmosphere  above  smoke  clouds 

Vertical  mixing  in  PBL  ceases  over  land  during  first 
week  of  nuclear  winter  due  to  stable  inversion. 

Because  new  method  mixes  less  efficiently  than  DCA, 
hot  spots  in  the  atmosphere  are  more  buoyant. 

Smoke  lofting  is  fester  than  in  old  model  with  DCA. 


(Ton  24  -  Tvi 


Subsurface  Heat  and  Water  Storage  and  Diffusion 

•  For  temperature,  6  nonuniform  levels  to  5m 

•  For  water,  2-layer  "bucket  model" 

•  Replaces  simple  energy  balance  assumption  (no  soil 
heat  capacity)  and  fixed  wetness 

•  Thermal  inertia  needed  for  diurnal  cycle 

•  Upward  diffusion  supplies  heat  to  surface  when 
radiation  deficit  occurs  (polar  night  and  nuclear  winter) 


•  When  atmosphere  is  stable,  upward  diffusion  of  heat 
approximately  compensates  for  absence  of  erroneous 
heat  flux  supplied  in  old  model  by  "neutral-stability 
linear  vertical  diffusion". 
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Modification 

of  predicted  surface  cooling 
by  model  improvements 
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AT  =  change  in  surface  air  temperature 
relative  to  an  unperturbed  July. 
averaged  over  the  USA. 
averaged  over  days  5-15. 


AT(°C) 

5(AT) 

"Old  Model" 

-15 

heat  and  moisture 
storage  in  soil 

-10 

+5 

stability-dependent 
vertical  diffusion 

-11 

-1 

absorption  of  IR 
by  smoke 

-8 

+3 

Lagrangian  transport 
of  smoke 

-8 

0 

diurnal  cycle 

-8 

0 
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Seasonal  dependence 
of  predicted  surface  cooling 


AT  =  change  in  surface  air  temperature 
relative  to  an  unperturbed  season. 
averaged  over  the  USA. 
averaged  over  davs  5-15. 


AT(°C) 

July 

-8 

October 

-6 

April 

-5 

January 

-4 
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200  day  evolution 
of  predicted  surface  cooling 
(Jun  29  -  Jan  15) 

AT  =  change  in  surface  air  temperature 
relative  to  an  unperturbed  season, 
averaged  over  the  USA, 
averaged  over  a  ten-dav  period. 

x  =  optical  depth, 

averaged  over  the  USA. 
averaged  over  the  ten-dav  period. 


Days 

AT(°C) 

X 

5-15  (Jul  4  -  Jul  14) 

-8 

2.0 

30-40  (Jul  29 -Aug  8) 

-1 

0.6 

190-200  (Jan 5 -Jan  15) 

-4 

0.4 

1 
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A  Theory  of  Radiative-Dynamical  Instability* 


STEVEN  J.  GHAN 

Lawrence  Livermore  National  Laboratory 
Livermore,  CA  94550 


ABSTRACT 


Numerical  simulations  of  the  global  scale  response  to  summertime  injections  of  absorp¬ 
tive  smoke  into  the  atmosphere  exhibit  a  self-lofting  mechanism,  in  which  radiative  heating 
and  the  dynamical  circulation  interact  to  loft  the  smoke  to  ever  higher  altitudes.  After 
making  several  plausible  approximations,  I  shall  demonstrate  that  for  small  perturbations 
this  mechanism  is  an  unstable  interaction. 

Linearized  equations  governing  the  conservation  of  vorticity,  enthalpy,  and  aerosol 
concentration  are  solved  analytically  under  certain  assumptions,  producing  an  expression 
for  the  exponential  growth  rate  of  internal  waves  propagating  in  a  stratified  fluid  with 
a  strong  vertical  gradient  in  the  distribution  of  solax^absorbing  aerosols.  Comparison  of 
the  analytical  growth  rate  with  smoke  lofting  rates  simulated  by  GCMs  yields  reasonable 
agreement,  thus  providing  theoretical  support  for  the  numerical  results. 


*This  work  was  performed  under  the  auspices  of  the  U.S.  Department  of  Energy  by  the 
Lawrence  Livermore  National  Laboratory  under  Contract  W-7405-Eng-48. 
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Unstable  Radiative-Dynamical  Interactions 

Steven  J.  Ghan 

Lawrence  Livermore  National  Laboratory 

An  outstanding  feature  of  the  terrestrial  and  Martian  atmospheres  is  their  near¬ 
transparency  with  respect  to  solar  radiation.  Although  terrestrial  water  clouds  scatter 
a  significant  fraction  of  the  incoming  solar  radiation,  and  ozone  absorbs  much  of  the  ul¬ 
traviolet  radiation,  most  absorption  of  visible  radiation  normally  occurs  at  the  surface. 
Radiative  heating  rates  throughout  most  of  the  troposphere  are  dominated  by  infrared 
radiative  cooling. 

In  the  Martian  atmosphere  this  situation  is  ocassionally  disrupted  by  global-scale  dust 
storms,  which  increase  the  solar  opacity  of  the  atmosphere,  leading  to  substantial  tropo¬ 
spheric  warming.  Although  such  global  storms  fortunately  do  not  develop  in  the  terrestrial 
atmosphere,  recent  numerical  simulations  (Malone  et  al.,  1986)  involving  the  terrestrial  at¬ 
mospheric  response  to  massive  injections  of  absorptive  smoke  produced  by  hypothetical 
post-nuclear  war  fires  have  exhibited  a  similar  phenomenon.  Substantial  lofting  of  the 
smoke  is  found  to  occur  in  these  simulations,  suggesting  that  the  feedback  between  short¬ 
wave  radiative  heating  and  the  dynamical  response  to  the  heating  can  be  important. 

Although  the  interaction  between  the  aerosol  distribution,  the  shortwave  radiative 
heating,  and  the  dynamical  circulation  in  Martian  dust  storms  and  in  hypothetical  terres¬ 
trial  smoke  storms  is  appreciated  in  the  literature,  the  unstable  nature  of  the  interaction 
has  not  been  rigorously  demonstrated.  Gierasch  et  al.  (1973)  discuss  a  radiative  instability 
mechanism  involving  clouds,  but  it  is  driven  by  longwave  radiation.  Gierasch  and  Goody 


(1973)  define  a  growth  phase  for  Martian  dust  storms,  but  treat  the  dynamical  circulation 
diagnostically  rather  than  prognostically.  Leovy  et  al.  (1973)  treat  the  dynamical  circula¬ 
tion  prognostically,  but  simply  prescribe  the  aerosol  distribution.  Aside  from  these  early 
analytical  attempts  to  explain  extraterrestial  cloud  bands  and  Martian  dust  storms,  most 
subsequent  work  on  the  radiative-dynamical  interaction  has  been  numerical  (Haberle  et 
al.,  1982;  Haberle  et  al.,  1985;  Malone  et  al.,  1986).  Thus,  a  rigorous  theoretical  demon¬ 
stration  of  the  shortwave  radiative-dynamical  instability,  and  an  analysis  of  the  conditions 
under  which  it  is  important,  does  not  presently  exist.  Given  the  current  interest  in  Martian 
dust  storms,  hypothetical  terrestrial  smoke  storms,  and  radiative-dynamical  interactions 
in  general,  such  a  theory  is  desirable.  Here  I  shall  summarize  the  progress  I  have  made  in 
developing  just  such  a  theory. 

The  theory  is  based  on  a  variety  of  simplifying  approximations.  Necessary  assumptions 
are  that  the  atmosphere  is  hydrostatic,  that  perturbations  to  an  assumed  basic  state  are 
small,  that  short-wave  scattering  and  infrared  emission  by  the  absorber  are  unimportant, 
and  that  the  basic  state  is  horizontally  uniform  with  no  vertical  variations  in  the  basic 
state  flow,  the  Brunt-Vaisala  frequency  N,  or  in  the  radiative-dynamical  feedback  rate 

_  RSqclT  dq 
“  ~  ~  cpN2H  di 

where  R  is  the  gas  constant,  So  is  the  solar  constant,  a  is  the  absorption  cross-section  of  the 
absorber,  T  is  the  basic  state  short-wave  transmissivity,  cp  is  the  specific  heat  at  constant 
pressure,  H  is  the  density  scale  height,  z  is  altitude,  and  q  is  the  basic  state  absorber  mixing 
ratio.  Additional  simplifying  assumptions,  which  can  be  relaxed  for  a  more  general,  but 
complicated  analysis,  are  that  solar  absorption  is  “grey”,  that  damping  of  perturbations  is 
equally  rapid  for  momentum,  heat,  and  absorber  mixing  ratio,  that  density  is  constant  (the 
Boussinesq  approximation),  and  that  motions  are  quasi-geostrophic  and  localized  enough 
to  neglect  the  sphericity  of  the  earth  and  latitudinal  variations  in  the  Coriolis  parameter 


The  general  solution  procedure  is  to  first  express  the  solar  heating  in  terms  of  pertur¬ 
bations  in  the  distribution  of  absorber  q.  The  equations  governing  conservation  of  mass, 
momentum,  heat  and  absorber  mixing  ratio  are  then  linearized  about  the  basic  state.  The 
linear  system  of  equations  is  then  reduced  to  a  single  partial  differential  equation  (PDE). 
By  expressing  solutions  in  terms  of  orthogonal  basis  functions  with  a  time  dependence 
given  by  exp(-tot),  the  PDE  reduces  to  an  algebraic  equation  for  a.  In  the  absence  of 
dissipation,  solutions  are  unstable  if  Im(o)  >  0. 

In  the  absence  of  radiative-dynamical  feedback  (a=0),  solutions  correspond  to  the 
Rossby  mode,  the  advective  mode,  and  (for  the  primitive  equations)  two  inertia-gravity 
modes.  The  Rossby  and  inertia-gravity  modes  propagate  with  respect  to  the  basic  state 
flow,  while  the  advective  mode  does  not.  All  of  these  modes  are  neutrally  stable. 

In  the  presence  of  radiative-dynamical  feedback,  the  advective  mode  becomes  unstable 
for  positive  a  (ie.,  basic  state  absorber  decreasing  with  altitude).  For  a  given  feedback 
rate,  the  growth  rate  of  disturbances  is  found  to  be  greatest  for  short,  deep  perturbations, 
but  never  exceeds  the  feedback  rate.  For  feedback  rates  less  than  the  Coriolis  parameter  /, 
the  quasi-geostrophic  theory  is  accurate,  while  for  feedback  rates  greater  than  the  Rossby 
wave  frequency  the  assumption  of  constant  /  is  adequate.  For  feedback  rates  greater  than 
/  but  less  than  the  Rossby  wave  frequency,  the  growth  rate  is  proportional  to  the  feedback 
rate. 

For  negative  a  (absorber  concentration  decreasing  with  altitude),  the  advective  mode 
is  damped,  but  the  Rossby  and  inertia-gravity  modes  are  unstable.  However,  the  net 
vertical  transport  of  absorber  mixing  ratio  for  these  circulation  modes  is  downward  rather 
than  upward,  so  they  cannot  be  responsible  for  the  smoke  lofting. 

The  feedback  rate,  because  it  characterizes  the  growth  rate  of  perturbations,  is  the 
most  important  parameter  of  the  problem.  If  one  assumes  that  the  absorber  mixing  ratio 
decreases  exponentially  with  altitude,  the  feedback  rate  is  found  to  be  largest  when  the 
absorption  optical  depth  equals  the  cosine  of  the  solar  zenith  angle.  Above  this  level, 


the  vertical  gradient  of  absorber  mixing  ratio  becomes  small,  while  below  this  level,  the 
transmission  of  solar  radiation  becomes  small.  For  typical  terrestrial  values  for  R,  S0,  a, 
cp,  N,  and  H ,  the  maximum  feedback  rate  (ie.,  when  the  zenith  angle  is  zero  and  the 
absorption  optical  depth  is  unity)  is  found  to  be  3.4  x  10-6  s-1  for  an  absorber  scale 
height  of  10  km,  and  1.7  x  10-4  s-1  for  an  absorber  scale  height  of  1  km,  corresponding 
to  time  scales  of  3.4  days  and  1.6  hours,  respectively.  These  feedback  rates  are  consistent 
with  the  time  scales  of  smoke  lofting  exhibited  in  the  GCM  simulations  of  Malone  et  al. 
(1986). 
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CBKONIC  EFFECTS  07  LAMS  ATMOSHBXC  SHOO  BUECTIOW:  DfflH- 

actions  win  ns  ocean  mixed  um,  sia  in,  and  giqund 

HYD10L0GY 


Steven  J.  Ghan,  Michael  C.  MacCracken,  Jobs  J.  Waltoa  (all  at: 
Lawrence  Liveraore  National  Laboratory,  L-282,  P.0.  Boa  BOS, 
Liveraore,  CA  94550). 

A  sot  of  extended  integrations  is  porforaod  using  a  gosoral 
circulation  sodol  of  ths  troposphere  asd  ocean  nixed  layer 
initialised  aith  a  largo  suasertiae  injection  of  absorptive 
saoke  over  North  Aaoriea  aad  Eurasia.  The  saeks  is  ouhject  to 
coagulation  aad  reaoval  as  sell  as  advective  aad  convective 
transport. 

Several  aspects  of  the  eliaatic  response  support  ths  con¬ 
clusions  of  lobock  using  as  energy  balance  eliaate  aodel.  Ths 
ocean  aixed  layer  in  northern  aidlatitudss  cools  5-10  K  rela¬ 
tive  to  the  control  eliaate  vitkia  three  neat  he.  During  the 
northers  sister  following  ths  injection,  this  cooling  decrea¬ 
ses  as  a  result  of  the  redaction  is  sasko  aad  ths  saallsr 
vintertiao  insolation.  During  ths  following  northern  manor, 
enhanced  sea  ice  coverage  insulates  the  ocean  sized  layer  iron 
the  ataoophere,  hence  redaciag  senasrtiae  earning  at  the 
ooraal  ice  aargia.  During  the  suhoeqpant  viator,  ths  eahssrsd 
sea  ice  coverage  redness  the  surface  boot  rapacity,  hence 
increasing  vintertiao  cooling  of  ths  ocean  surfaoo,  but 
reducing  the  cooliag  of  ths  ocean  aixed  layer. 

Other  aspects  of  ths  response  are  quite  different,  however. 
Due  to  a  saoke-iaduced  redaction  in  land  precipitation,  sub¬ 
stantial  surface  warning  aad  drying  occurs  ia  ths  tropica 
during  the  first  aoaths  following  ths  injection.  Ia  addition, 
the  ocean  aixed  layer  deepens  considerably  an  a  result  of 
•ntraiaaoat  induced  by  both  greater  surface  cooliag  aad 
enhanced  aschsaicsl  stirring. 


Ibis  work  was  perfomod  under  ths  auspices  of  ths  0.8.  Depart* 
««IF  ky  ths  L serenes  Liveraore  Rational  Lahoratery  under  Goatn 
405-8ac-48. 


INFLUENCE  OF  SOLAR  HEATING  ON  SCAVENGING 
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Fig.  4.  Time  evolution  of  the  calibrated  (a)  net  surface  heat  flux  Q,  and  (6)  surface  wind  stress  t  at  JOHN,  150*E  used  in 
the  control  experiment  (solid  lines)  and  the  nuclear  winter  experiment  (dashed). 


By  examining  Table  2  and  Figure*  2  and  3,  it  is  possible  to 
determine  for  each  geographical  location  whether  the  model- 
simulated  SST  decrease  during  the  nuclear  winter  experiment 
was  due  to  enhanced  air-sea  cooling  or  entrainment.  For  ex¬ 
ample.  at  the  first  location  (S0°N,  150°W),  the  net  heat  flux 
into  the  ocean  was  33  W  ra~a  less  during  the  nudear  winter 
experiment  than  during  the  control  experiment  (Table  2).  If 
the  average  MLD  at  this  location  during  nudear  winter  is 
taken  to  be  50  m  (Figure  3),  this  decrease  in  surface  heat  flux 
cannot  account  for  more  than  a  1°C  difference  in  SST  at  the 
end  of  the  month.  This  implies  that  well  over  half  of  the  2.8'C 
decrease  in  SST  in  the  nuclear  winter  experiment  (Figure  2) 
was  caused  by  entrainment  This  interpretation  is  also  consis¬ 
tent  with  the  much  greater  wind  stress,  an  important  (actor  in 
produdng  entrainment,  experienced  during  nuclear  winter  at 
this  location  (Table  2).  In  the  same  way,  we  find  that  more 
than  half  of  the  decrease  in  SST  at  the  second  geographical 
location  (50°N,  150°E)  is  due  to  a  decrease  in  the  net  surface 
heat  flux,  with  the  rest  due  to  enhanced  entrainment  (see 
Figure  5b,  and  additional  discussion  below).  At  the  third  lo¬ 
cation  (50*N,  50°  W),  where  there  was  little  change  in  the  wind 
stress  and  MLD  due  to  nuclear  winter,  all  of  the  SST  decrease 
is  attributed  to  the  rather  large  decrease  (126  W  m"1)  in  the 
net  surface  heat  flux.  The  only  other  location  where  the  SST 
response  was  significant  is  the  East  China  Sea  point  (30°N. 
I23°E),  and  there  most  of  the  SST  decrease  is  accounted  for 
simply  by  the  very  large  change  in  the  net  surface  beat  flux 
caused  by  the  extraordinary  sensible  and  latent  heat  fluxes 
noted  above. 

We  now  examine  the  response  of  the  upper  ocean  to  nuclear 
winter  in  greater  detail  by  describing  the  time  evolution  of  the 
vertical  thermal  structure  and  its  relation  to  the  atmospheric 
forcing  at  one  of  the  experimental  ocean  locations.  The  point 
at  50*N,  130*E  is  chosen  for  this  purpose  because  the  atmo¬ 
spheric  forcing  during  the  control  and  nuclear  winter  experi¬ 
ments  at  this  location  is  fairly  representative  of  the  other  lo¬ 
cations  as  well.  This  point,  located  approximately  300  km  west 
of  the  southern  tip  of  the  Kamchatka  peninsula  in  the  sea  of 
Okhotsk,  experiences  a  doubling  of  the  average  wind  stress 
and  a  139  W  m'*  decrease  in  the  net  surface  heat  flux  due  to 
the  smoke  injection  (Table  2).  The  results  for  this  point  are 
shown  in  Figures  4  a:  d  5. 


Looking  first  at  the  atmospheric  forcing  (Figure  4),  we  see 
that  the  net  heat  flux  is  affected  by  the  smoke  as  early  as  July 
3,  and  that  an  exceptionally  strong  wind  event  lasting  several 
days  occurs  shortly  thereafter.  For  16  consecutive  days  during 
the  nuclear  winter  experiment,  from  July  II  to  26,  the  solar 
radiation  is  reduced  essentially  to  zero  by  the  thick  smoke 
cloud  originating  over  nearby  Siberia,  and  this  produces  a  net 
heat  loss  (Q.  <  0)  which  lasts  almost  to  the  end  of  the  month. 
By  the  end  of  the  month,  the  net  heat  flux  has  partially  re¬ 
turned  to  normaL  Except  for  the  enormous  wind  storm  during 
the  first  week,  the  wind  stress  is  quite  similar  in  the  two  exper¬ 
iments.  This  characteristic  of  a  strong  wind  event  early  in  the 
month,  followed  by  severe  reductions  in  solar  radiation  and 
subsequent  partial  recovery,  is  typical  of  the  atmospheric  forc¬ 
ing  during  nudear  winter  at  other  locations  as  welL 

The  response  of  the  upper  ocean  thermal  structure  (Figure 
3)  is  easily  understood  in  terms  of  the  above  forcing.  During 
the  control  experiment,  the  upper  ocean  responds  to  a  cycle  of 
typical  summer  wind  events  on  July  6, 13, 16, 21,  and  23.  The 
ocean’s  response  to  these  wind  events  is  seen  in  the  deepening 
and  coale  scans*  of  the  isotherms  representing  an  enhanced 
vertical  temperature  gradient  at  the  base  of  the  mixed  layer 
caused  by  the  wind-generated  downward  mixing  of  warm  sur¬ 
face  layer  water.  In  between  these  characteristically  gentle 
wind  events,  the  isotherms  tend  to  spread  out  in  the  vertical  as 
the  mixed  layer  shallows  and  warms  owing  to  the  net  surface 
heating.  The  evolution  during  the  nuclear  winter  experiment  is 
very  different  The  wind  storm  on  July  5  deepens  and  cools  the 
mixed  layer,  and  the  subsequent  negative  (upward)  heat  flux 
and  relatively  normal  winds  produce  a  gradual  cooling  and 
slow  deepening  of  the  mixed  layer  during  the  remainder  of  the 
month.  In  spite  of  the  drastic  changes  in  the  atmospheric  forc¬ 
ing  during  the  nuclear  winter  experiment  the  response  of  the 
ocean  is  confined  to  the  upper  60  m. 

Discussion 

In  this  study  we  have  examined  the  short-term  response  of 
the  upper  ocean  to  the  simulated  forcing  of  a  hypothetical 
atmospheric  nuclear  winter.  The  atmospheric  forcing  was 
derived  from  fields  produced  by  a  two-level  atmospheric  gen¬ 
eral  circulation  model  in  both  a  control  and  a  nuclear  winter 
simulation.  The  ocean  was  represented  by  a  one-dimensional 


231 


MrmACH  et  «l:  Brief  Retort 


fate  4,s 

OC«an  mixed  lay f  response 
cm  por+om  t? 


Uhect  Arc  jde  rclc s  &f  so*  /c* 

-&ri*?at/on  **/  hydro Icycal 

cyc.lt  /Vi  *^c  if/ron/c  res  pa  ^  re  7 


i 


234 


OSU/LLNL  TfcO POSPHEZIC  GCM 

LAGZA TZAce  SPECIES 
TtAMSPoZT 

POLLAZb  2-LAV£A  OCSAti 

•  hvxeb  layez 

•  THERMIC  CL/W6 

•  0 YMAtflCAL 

•  PAZArOETEUWE  0  EWTRAlbifflENT 

. THem o  b YtJA m/c  sea  ice 


ExDenmeot 


SURFACE  STRESS  DYNE/CM2 


MO  9 


z 

z 

z 

z 

z 

CO 

CO 

CO 

CO 

CO 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O) 

r*. 

IO 

CO 

CO 

to 

r*. 

O) 

259 


g> V 

•§  s 

.5  o 

6 


w 

p 

o 

S 

w 

H 

◄ 

2 

HH 

*4 

u 

Sr 

g 

o 

2 

w 

n3 

Z 

n9 

nJ 


Cti 

P 

O 

H 

£ 

<J 

ftj 

O 

m4 

►-3 

hJ 


V 
53 

o 

E 

V 
** 

o 

s 

CO 


ho 

c 


« 

09 

hJD 

C 

• 

hO 

C 

V 

§ 

O 

09 

** 

s 

o 

J2 

3 

ho 

c« 

o 

CJ 


u 

o 

a 

09 

c 

rt 

u. 

v 

j< 

o 

£ 

CO 

CO 

V 


J3 

£ 

O 

o 

cs 


£ 

rX 

o'  o 
o  o 
o  ^ 

^  X 

CO  £ 
oj 

K  © 

ccS  O 

a:* 


V 

CJ 

u. 

*» 

a 

_c 

rt 

« 


3 

cj 

V 

JO 

> 

3 

h 

o 

cfl 

u 

a 

V 

CCS 

> 

a 

XI 

g 

V 

o 

CJ 

o 

•m 

53 

S 

s 

V 

°  o 

.2 

■4A 

ccS 

Q) 

*w 

s  § 

-S  !S 

t§ 

O 

a 

t- 

V 

x— N  53 

V 

3 

CJ 

09 

1.8 
•  f. 

CD 

V 

g 

> 

s 
•  »— 

CQ 

V 

hO 

v 

o 

53 

S 

«3 


V 

£ 

3 

S 

O 

c 


CQ 

V 
•  ^ 

o 

V 

a 

CO 


09 

3 

E 


I* 

0) 


co 


V 

-X 

o 

£ 

CO 


JU 

3 

Cfl 


0)  a 

P  <3 


ns 

s 

4 

« 


a 

E  2 

<4  I* 

*  *3 

a  * 

.12  4 
h0  3  *-> 
0*  rt 

^  w  > 

§>•  • 

1-4 


a  §  o 
o  .2  s 

ho  r 

4  §1 

f!  i_.  W 
^  H  •_, 

rt  4 

8  J  t 

to" 

CSV 

•Go  .2  § 

g  O  53 

V  QJ 

5  >  I 

o  55  1 

CO  < 

•  • 


"  *V  G 
>»  *C  4 
.tS  ho  ^ 

5: 

<8  3  • - 

Sr«  55  73 

ho  ^ 

o  cr 

s  4 

-  2 

a  £ 

u  > — ' 


C- 

ho 

o 

LO 


09  h0 

QJ  C 

1:3 

-  E 

h  _ 
V  « 

ho  y 

c  *- 

o  4 

—  a 

i  u 

V 


x  “ 

o  .5 
^  XS 

g  O 

O  o 

g  E 

§  * 
53  ccS 

ho  .2 

3  « 

-  s 

S  c 


cj 

eS 

V 

V 

•  «N 

•P 

«S 

•  0«4 

•s 

U. 

V 

ho 

4 


cfl 

> 

o 

E 

V 

t-l 

■w 

S 

v 

cj 

09  CC 

a>  «2 

N  V 

*53  ® 

—  •— 

v  .. 
y  6 

’5  §. 

Ch  ^ 

a  \y 

o 

* 

H 


26 1 


mmmmmmmmm- 


>  1/im:  efficient  removal,  smaller  radiative  effect 
Scavenging 

Convective:  fractional  time  and  grid  area 
Large-scale:  fractional  time 


LLNL  SMOKE/CLIMATE  MODEL 


ULf4u  SauJinc  Com  ^T®  -8.Z  &T- -‘ZZ 


oe  325 


o 


•Vo 


o 


275 


t*iVt 


■'**&*) 


.  i V \\ v* :4*ax*%“*v»*#V>V 4*»%%V.  !iv$V»v*SV»V»V%V*,iV iVvV*VtVA**V  V* V»V<VtN' *tV< \$JV, 


p  I  o  i|  $  a  1 1|  j,  *o  |  a  g  i  u  s  a  j  s  j 

278 


S8S8a^^ 


282 


ssssssssss^^ 


sssssssssss^^ 


p  I  o  ij  s  a  j  ij  x  m  o  |  a  g  j  u  a  o  j  a  j 

284 


'W 


T'TTtTtvTT' 


W  f.*  =**» 

A'.VV 


;*w*:;v-/*  •  w-v** 


ssssssss^^ 


ssssssssssss^^ 


p  I  o  q  s  a  j  q  j.  mo  |  a  g  \  u  a  o  j  a  g 

285 


o> 


sssssssssss^^ 


Saaaaaaa^^ 


>,♦1 


j 


tin 


>  C^l 

SSSSSSSSSSSSS^^ 

J  ’’"“  o 

!«!l 


p  |  o  ij  s  i  j  i|  i  m  o  |  a  g  i  u  a  o  j  a  d 


293 


296 


Reprinted  from  Nature,  Vol.  325.  No.  6106.  pp  ?01-*03,  W  February  1V8" 
Macmillan  Journals  Ltd..  1987 
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Climate  perturbation  caused  by  mddea  injections  of  aerosols  into 
the  atmosphere  have  received  increased  atteatioa  with  the  realiza¬ 
tion  that  fires  started  by  a  nuclear  war  might  release  sufficient 
quantities  of  smoke  to  block  sunlight  from  the  Earth’s  surface 
and  dramatically  reduce  temperatures  over  land1-1.  Three- 
dimensional  atmospheric  general  circulation  models  which  include 
aerosol  transport  indicate  that  much  of  the  smoke  could  survive 
in  the  atmosphere  for  months  or  more4-*,  implying  significant 
climatic  effects  on  such  timescales,  but  these  models  generally 
assume  fixed  sea  surface  temperature  and  sen  ice,  and  therefore 
cannot  properly  simulate  cooling  for  mote  than  a  few  weeks.  This 
report  describes  a  general  circulation  model  run  in  which  the 
foregoing  assumptions  ate  relaxed  and  the  climate  impact  of 
massive  smoke  loading  is  smessfd  for  four  months  of  simulated 
time.  Two  competing  long-term  effects  are  apparent:  positive  feed¬ 
back  on  surface  cooling  due  to  enhanced  sea  ice  formation,  and 
surface  warming  at  some  latitudes  due  in  part  to  heat  transport 
from  the  upper  troposphere,  where  smoke  aerosols  absorb  sunlight 

It  is  instructive  to  compare  results  from  the  general  circulation 
model  (GCM)  to  those  obtained  by  Robock7  using  a  one¬ 
dimensional  (in  latitude)  energy  balance  model  (EBM)  in  multi¬ 
year  simulations  of  nuclear  war  climate  effects.  In  the  present 
study,  total  column  absorption  optical  depth  of  smoke  r  is 
prescribed  as  a  function  of  time  to  be  twice  the  value  normally 
assumed  by  Robock  except  that  we  restrict  r  *  3,  that  is,  under¬ 
neath  the  smoke  r  =  3  for  the  first  30  days  and  thereafter  r 
decreases  with  a  half  life  of  30  days.  This  prescription  ensures 
that  GCM  results  of  interest  are  obviously  significant  compared 
to  model-generated  variability  (that  is,  weather  fluctuations)  and 
also  that  during  the  first  20  days  after  aerosol  injection  the  smoke 
concentrations  assumed  for  this  study  are  identical  to  those  used 
in  our  earlier  GCM  experiments  with  prescribed  sea  surface 
temperature  (SST)*'9.  In  all  studies  discussed  here,  including 
Robock’s,  smoke  is  assumed  to  exist  in  a  uniform  band  of 
latitudes  between  30°  and  70°  N  and  is  injected  into  the  atmos¬ 
phere  at  approximately  the  Northern  Hemisphere  spring 
equinox.  In  the  GCM  runs  smoke  is  assumed  to  be  distributed 
between  1  and  10  km  altitude  at  constant  mixing  ratio. 

The  assumed  initial  condition  r  >  3  in  the  Northern  Hemi¬ 
sphere  mid-latitudes  is  equivalent  to  34  x  10>:  g  of  amorphous 
elemental  carbon  with  absorption  cross  section  I0mJg-1  and 
corresponds  to  ‘baseline’  estimates  of  nuclear  war  smoke  pro¬ 
duction  published  by  the  US  National  Academy  of  Sciences3 
and  the  International  Council  of  Scientific  Unions1.  However, 
longer  term  optical  depths,  r»2.4  at  40  days  and  0.5  at  110 
days,  greatly  exceed  values  implied  by  GCMs  which  compute 
smoke  removal1’*'10,  and  or  course  the  assumption  that  smoke 
is  fixed  in  space  is  unrealistic.  Thus,  the  present  study  is  an 
intercomparison  of  climate  model  responses  to  hypothetical 
perturbations,  with  application  to  nuclear  war  climate  effects; 
it  is  not  an  attempt  to  directly  simulate  'nuclear  winter’.  The 
use  of  extreme  smoke  amounts  and  the  assumption  that  smoke 
docs  not  spread  beyond  the  Northern  Hemisphere  mid-latitudes 
mean  that  long-term  temperature  effects  are  probably  overesti¬ 
mated  by  the  present  study.  Indeed,  in  reporting  a  long-term 
‘nuclear  winter’  simulation  including  smoke  spreading  by  a 
two-level  GCM,  Stenchikov"  does  not  elucidate  the  surface 
temperature  feedback  processes  discussed  below,  perhaps 
because  these  processes  are  less  obvious  with  more  realistic 
(smaller)  smoke  concentrations. 

The  GCM  used  in  this  study  is  the  NCAR  (National  Center 


Time  Irom  smoke  iniecnon  (days) 


Fig.  I  Zonaily  averaged  land  surface  temperature  difference  in 
°C  between  prescribed  smoke  experiments  and  controls.  The  upper 
panel  shows  results  obtained  in  the  present  study;  the  lower  panel 
is  redrawn  from  ref.  7  and  shows  results  from  a  one-dimensional 
energy  balance  climate  model.  In  both  cases  solid  line  contours 
are  drawn  at  10  “C  intervals,  and  shaded  areas  indicate  tem¬ 
peratures  more  than  10  °C  cooler  than  control  runs  while  hatched 
areas  indicate  temperatures  more  than  10  ”C  warmer  than  control 
runs.  GCM  results  are  shown  for  130  days  and  the  thin  dashed 
lines  indicate  the  corresponding  time  period  in  Robock's  simula¬ 
tion.  Results  from  the  GCM  are  not  plotted  at  latitudes  which  are 
more  than  90%  ocean. 

for  Atmospheric  Research)  Community  Climate  Model 
(CCM)12.  Results  presented  here  were  obtained  with  a  research 
version  of  the  CCM  developed  by  Washington  and  Meehl13,14 
for  studies  of  the  C02  greenhouse  effect.  In  this  model  SSTs 
are  computed  assuming  a  static  mixed  layer  of  constant  heat 
capacity,  and  sea  ice  is  computed  with  a  simple  thermodynamic 
model.  This  GCM,  hereafter  referred  to  as  the  CCM/mixed 
layer  model,  was  modified  as  described  by  Covey  el  al *'9  to 
include  absorption  of  sunlight  by  smoke.  Similar  representations 
of  mixed  layer  heat  capacity  and  sea  ice  are  used  in  Robock's 
EBM.  (The  earlier  studies  of  Covey  et  al'-9  made  use  of  a  version 
of  the  CCM  in  which  SST  and  sea  ice  are  prescribed  at  seasonally 
varying  climatological  values11;  this  model  will  be  referred  to 
as  the  prescribed  SST  version.) 

Figure  1  shows  the  smoke-induced  land  surface  cooling  (that 
is,  the  temperature  difference  in  °C  between  prescribed  smoke 
experiments  and  controls)  averaged  over  longitude  and  plotted 
as  a  function  of  latitude  and  time.  The  upper  panel  shows  results 
from  the  CCM/mixed  layer  model  and  the  lower  panel  shows 
results  from  Robock's  EBM  (‘Run  4‘  of  ref.  7).  During  the  first 
half  of  the  120  day  simulation  with  the  CCM/mixed  layer  model, 
land  surface  cooling  underneath  the  assumed  smoke  cloud  is 
about  twice  that  exhibited  by  Robock's  model.  Maximum  cool¬ 
ing  in  the  GCM  is  47  °C  at  day  37  and  60*  N  latitude  (poleward 
of  the  sea  ice  margin,  slightly  to  the  north  of  maximum  cooling 
in  Robock's  model).  During  this  time  period  smoke  optical 
depths  assumed  in  both  the  CCM/mixed  layer  model  and 
Robock's  simulation  are  >1,  so  the  discrepancy  in  the  amount 
of-civ,!"g  between  the  two  models  is  difficult  to  explain  by  the 
dilfcsctaJ in  smoke  loading.  The  CCM/mixed  layer  model  may 
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Fl«.  2  Soil  moisture  content  averaged  over  longitude  and  over 
days  1 10  to  120  after  smoke  injection  (that  is,  the  final  10  days  of 
the  GCM  simulation)  for  the  smoke-perturbed  case  (solid  line) 
and  the  control  case  (dashed  line).  Results  are  not  plotted  at 
latitudes  which  are  more  than  90%  ocean. 

exaggerate  the  amount  of  cooling  because  it  overpredicts  the 
unperturbed  amount  of  sea  ice,  possibly  due  to  neglect  of 
poleward  heat  transport  by  the  oceans13,14;  the  maximum  tem¬ 
perature  drop  exhibited  by  the  CCM/mixed  layer  model  is 
almost  twice  that  of  the  prescribed  SST  version  during  the  time 
that  smoke  loadings  in  the  two  GCMs  are  identical.  On  the 
other  hand,  Robock's  model  may  well  underpredict  the 
maximum  degree  of  land  surface  cooling  because  its  limited 
resolution —  10°  latitude,  no  longitude  resolutions  within  land 
and  ocean  sectors,  no  height  resolution — requires  the  cooling 
to  be  spread  out  over  a  large  volume  of  space.  The  latter 
consideration  may  also  explain  in  part  the  rapid  cooling 
exhibited  by  the  CCM  during  the  first  few  days  after  aerosol 
injection';  the  EBM,  in  contrast,  takes  several  weeks  to  obtain 
land  surface  cooling  >10°C.  Another  factor  contributing  to 
rapid  initial  cooling  in  the  GCM  is  atmospheric  dynamics,  not 
included  in  the  EBM.  An  examination  of  the  lower  tropospheric 
wind  and  temperature  fields  (not  shown)  indicates  that  ocean-to- 
land  temperature  advection  takes  about  two  weeks  to  increase 
to  <ts  maximum  value  as  the  circulation  in  the  perturbed  case 
evolves  to  a  more  zonally  symmetric  state.  Two  weeks  is  approxi¬ 
mately  the  thermal  relaxation  time  of  the  lower  troposphere, 
that  is,  the  quantity  CtTg~'\P/(<rT*)  where  C,  is  the  heat 
capacity  of  air  at  constant  pressure,  T  is  the  air  temperature,  g 
is  the  acceleration  of  gravity,  \P  is  the  pressure  difference 
between  the  surface  and  the  85  kPa  pressure  level,  and  a  is  the 
Stefan-  Boltzmann  constant. 

Unlike  the  situation  early  on,  from  about  100  to  120  days 
after  smoke  injection  the  amount  of  cooling  in  the  CCM/mixed 
layer  model  is  about  as  small  or  smaller  than  that  in  Robock's 
EBM,  despite  doubled  smoke  loading.  Furthermore,  throughout 
the  second  half  of  the  CCM  simulation  land  surface  tem¬ 


peratures  immediately  south  of  the  assumed  smoke  cloud  are 
significantly  warmer  than  control  values  (see  the  hatched  areas 
in  Fig.  1).  The  EBM's  surface  temperatures  are  everywhere 
colder  when  smoke  is  injected  into  the  model,  for  the  following 
reason.  Since  the  EBM  has  no  height  resolution  the  method  of 
assuming  smoke  is  to  simply  decrease  the  solar  constant  by  the 
amount  corresponding  to  absorption  by  the  appropriate  optical 
depth.  This  situation  is  quite  different  from  that  expected  in  the 
real  atmosphere,  in  which  smoke  warmed  by  absorption  of 
sunlight  should  overlie  a  cool  surface  (see  Fig.  la  in  ref.  8).  It 
is  possible  that  in  the  CCM  simulation  some  of  the  heat  from 
the  upper  troposphere  is  eventually  transported  (presumably 
via  infrared  radiation16)  down  to  the  surface,  reversing  long-term 
surface  cooling  in  those  areas  which  are  not  covered  by  smoke. 
This  effect  would  no  doubt  be  smaller  in  a  more  realistic  treat¬ 
ment  of  protracted  nuclear  war-generated  smoke,  which  would 
be  more  widespread,  thinner,  less  black,  and  higher  than  the 
idealized  fixed  smoke  considered  here;  nevertheless  it  is  possible 
that  the  effect  could  moderate  long-term  surface  cooling. 

Partial  support  for  the  foregoing  hypothesis  comes  from  an 
analysis  of  the  GCM’s  land  surface  energy  budget  averaged  over 
the  final  10  days  of  the  model  run.  At  20°  N  latitude  infrared 
flux  from  the  atmosphere  to  land  surface  is  92  W  m“2  greater  in 
the  smoke-perturbed  case  than  in  the  control  case.  A  smaller 
warming  effect  (28  W  m~2)  is  produced  by  loss  of  upward  latent 
heat  flux  due  to  a  drying  out  of  soil  moisture  in  the  15-30°  N 
latitude  zone  (Fig.  2).  The  latter  type  of  wanning  mechanism 
has  been  observed  in  a  previous  GCM  simulation  of  nuclear 
war-generated  smoke  effects17.  However,  underneath  the  smoke 
(47°  N)  the  dominant  change  in  surface  energy  balance  is 
between  decreased  absorption  of  solar  radiation  and  decreased 
upward  infrared  and  sensible  heat  fluxes;  within  this  latitude 
band,  the  infrared  flux  received  by  the  surface  from  the  atmos¬ 
phere  is  not  significantly  greater  in  the  smoke-perturbed  case 
than  in  the  control  case. 

The  differences  in  soil  moisture  between  the  control  and 
perturbed  cases  (Fig.  2)  are  for  the  most  part  as  expected  from 
previous  experiments  with  the  CCM9.  In  general,  the  CCM 
responds  to  fixed  upper  tropospheric  smoke  loadings  by  decreas¬ 
ing  cloudiness  where  the  smoke  is  heated  by  sunlight  while 
increasing  cloudiness  at  lower  levels,  resulting  in  only  small 
changes  in  precipitation  underneath  the  smoke.  (This  result  is 
not  consistent  with  results  from  the.  Oregan  State  University 
two-layer  model'7  and  may  be  due  to  the  CCM's  assumptions 
of  100%  relative  humidity  immediately  above  the  surface 
together  with  boundary  layer  condensation  at  only  80%  relative 
humidity.)  Accordingly,  soil  moisture  is  greater  underneath  the 
smoke  in  the  perturbed  case,  where  surface  cooling  inhibits 
evaporation.  As  pointed  out  above,  soil  moisture  essentially 
disappears  from  latitudes  1 5-30°  N,  where  surface  warming  from 
enhanced  downward  infrared  radiation  increases  evaporation 


Fig.  3  Surface  temperature 
difference  (prescribed 
smoke  case  minus  control 
case)  averaged  over  days 
tIOto  1 20  after  smoke  injec¬ 
tion.  On  the  left  side  the 
zonal  average  is  shown  as  a 
function  of  latitude.  On  the 
right  side  the  geographical 
distribution  is  shown,  with 
shaded  areas  indicating 
greater  than  10  °C  cooling 
and  hatched  areas  indicat¬ 
ing  greater  than  10X  warm¬ 
ing;  also  the  control  case 
(thick  dashed  line)  and  per¬ 
turbed  case  (thick  solid  line) 
sea  ice  margins  are  indi¬ 
cated  (the  nro  coincide  in 
the  Southern  Hemisphere). 


during  the  first  half  of  the  smoke-perturbed  model  integration. 
For  latitudes  15-30°  S.  the  expected  location  of  the  downward 
branch  of  an  enhanced  cross-equatorial  Hadley  circulation  in 
the  perturbed  case  (see  Fig.  4ft  in  ref.  8),  a  decrease  in  precipita¬ 
tion  leads  to  the  drying  out  of  soil  moisture  apparent  in  Fig.  2. 
Finally,  soil  moisture  more  than  doubles  in  the  perturbed  case 
in  a  narrow  latitude  band  centered  on  8°  N  due  to  a  burst  of 
precipitation  during  the  last  10  days  of  the  model  integration, 
which  occurs  for  unknown  reasons. 

The  geographical  distribution  of  surface  cooling  and  sea  ice 
averaged  over  the  final  10  days  of  the  CCM/ mixed  layer  model 
simulation  is  depicted  in  Fig.  3.  As  expected  from  Robock's 
EBM  results1  the  extent  of  sea  ice  is  significantly  greater  in  the 
smoke-perturbed  case.  The  GCM  obtains  greatest  cooling  in  the 
western  parts  of  the  North  American  and  Eurasian  land  masses, 
downwind  of  the  sea  ice,  in  contrast  to  early  times  when 
maximum  cooling  occurs  in  the  eastern  parts8.  Cooling  of 
Europe  due  to  enhanced  North  Atlantic  sea  ice  has  a  palaeocli- 
matic  analogue:  the  Younger  Dry  as  period,  10,000-11,000  years 
ago,  exhibited  this  phenomenon  according  to  both  observa¬ 
tional18  and  model19  studies. 
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In  summary,  the  three-dimensional  GCM  simulation,  while 
qualitatively  confirming  a  positive  feedback  due  to  sea  ice  that 
enhances  long-term  cooling  due  to  nuclear  war-generated 
smoke,  reveals  additional  effects  that  can  moderate  or  perhaps 
even  reverse  long-term  cooling.  The  latter  appear  to  dominate 
by  the  end  of  the  GCM  run,  when  the  extent  of  land  area  that 
is  cooled  in  excess  of  10°  C  is  much  smaller  than  it  is  in  the 
analogous  EBM  experiment  (‘Run  4‘  in  Fig.  1  of  ref.  7)  despite 
the  fact  that  the  GCM  experiment  assumed  more  smoke. 
However,  both  classes  of  effects  are  probably  exaggerated  by 
the  assumption  of  very  large  quantities  of  smoke  that  is  fixed 
in  space,  and  by  other  model  simplifications  such  as  lack  of 
ocean  heat  transport  (leading  to  overprediction  of  sea  ice).  It 
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ments,  particularly  including  model-computed  smoke  transport 
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ance  of  the  positive  and  negative  feedbacks  discussed  here. 

I  thank  Warren  Washington  and  Gerald  Meehl  for  providing 
both  the  model  and  valuable  advice  on  its  use,  and  Lynda 
Verplank  for  expert  assistance  with  programming.  This  work 
was  supported  by  the  NSF  (grant  ATM-8504807). 


10.  Thompson.  S.  L,  A  Giorgi,  F.  /  geophyx  Rex  (in  the  press). 

11.  Stenchikov,  G.  L  in  The  Night  After  . . .  Climatic  and  Biological  Consequences  of  Nuciear 

WV  <ed.  Velikov.  Y.)  53-82  (Mir.  Moscow.  1985). 

12.  Pitcher.  E.  J.  et  al  J  atmox  Sci  4B,  580-604  (1983). 

13.  Washington.  W.  M.  A  Meehl,  G.  A.  J  geophyx  Rex  B4.  9475-9503  (1984) 

14.  Meehl.  G.  A.  A  Washington,  W.  M  J.  phyx  Oceanogr.  15,  92-104  ( 1985). 

15.  Chervin.  R.  M.  /  atmos.  Sci  43,  233-251  (1986). 

16.  Cess.  R.  D.,  Potter,  G.  L.  Ghan.  S. A  Gates.  W  L.  J.  geophyx  Rex  9B.  1 2937- 1 2950  ( 1985 ) 

17.  Ghsn,  S.  J..  MacCracken,  M.  C.  A  Walton.  J.  J.  Lawrence  Uoermore  nam.  Laboratory  tech 

Not  UCRL-92324  ( 1985). 

18.  Broecker.  W.  S.,  Pweet.  D.  M.  A  Rind.  D.  Nature  JIS.  21-26  (1985) 

19.  Schneider.  S.  H.,  Peteet,  D.  M.,  A  North.  G.  R.  in  Abrupt  Climate  Changes  feds  Berger. 

W.  H.  A  Lsbeyne,  L.)  (Reidcl,  Dordrecht,  in  the  press). 


i 


30J 


Mod*l  -for  S4adjj  at  prolonged 
( 10  monHt;  Clinrilc  impact  of 
odbmospkerlc  Smoke.  injections 


Climate  model  (  WaskiVg+ori  V  Meek!  WM) 

3  O  Qjtmogphe>'/c  ft  CM  ( MCAR  CC.M) 

ShotocoueA,  /  Soil  moisfiune  (  Ma/noJbe) 
*tWmodjjr\<wn‘it  Sea.  ice  model  (Semtrver) 
mixed.  (ajcA,  keaJb  (joupaciX^  E  “oc ea*,  model* 

Assumptions  dbouJt  Smoke  (Covey  eJb  oJL.  19?^) 

Puea.  in  Spate:  1-i.Okrvi,  30°W  -  ?0°M 

VaAAAfole  in  -time 

CoAs+QM3t  for  first  30  dUijs 
TKen.  kaaJ-Sty*  crj  3odLo^s 
Ve/iy  Ia/i0e  Smoke.  OLmouA^S 

Absorption,  optical  deptk  =  3  fin.  30  dc uj 
"^e/EaffeA.,  Smoke.  ~  3.X  TTAPS  baseline 


302 


jj*A\ 


ft 


iSi 


v*l 


m 


m 


: 


MS 


j 


V 


•.('I 


ss 


m 


m 


8 


isa 


m 


Mi!* 


cm 


» M 


,va 


£S1 


t.' 


'♦.R 


rJ 


‘‘.*1 


w 


C.HAM&E  IN  UNO  SURFACE  TEMPERATURE 


(Fe/turbtA  Case)  -  (Cow+rol  Case) 


EQUATOR 


EQUATOR 


TIME  FROM  SMOKE  INJECTION  (DAYS) 


u 


«20l 


3* 


wm 


spy  w 


'«  *ofMv>r>  *  a  to  i .  ^ 


3D  GCM 


COKJCLULS/OKiS 


1 •  Long  -  -ie/tm  feejibacks 

CL.  Sea.  102.  -feceLbadc  6nJi ounces  Coo/ iruj  * 

O-  LLppeA.  iroposphwk,  Leodt  diminishes  Coo/ii^ 

6o4K  oJ-f&dts  ojul,  Lkaej^tAoltd.  by  iluu  id&alijid 
OSSujmpiions  CLboujt  Smoke. 

X.  Ctaunges  in  kycUdlopy  budg&f 


Cj.  Roback, 


306 


Effects  of  Vertical  Injection  Profile  and  Infrared  Opacity 
in  GCM  Simulations  of  Nuclear  War  Smoke 


Stephen  H.  Schneider 

National  Center  for  Atmospheric  Research 
P.O.  Box  S000,  Boulder,  Colorado  80S07 


A  global  general  circulation  model  (GCM)  that  allows  for  the  transport,  removal 
and  radiative  effects  of  smoke  aerosols  has  been  used  to  examine  the  effects  of  initial 
smoke  vertical  injection  profile  and  smoke  thermal  infrared  opacity  on  the  subsequent 
evolution  of  land  surface  temperatures.  Two  vertical  injection  profiles  are  considered: 
constant  mass  mixing  ratio  (CMR)  from  0-7  km  altitude,  and  constant  density  (CD) 
concentration  from  0-10  km  altitude.  Both  profiles  prescribe  the  same  column  loading.  A 
variety  of  simulations  have  been  done  with  both  profiles.  In  general,  the  CD  cases  produce 
midlatitude  land  temperature  decreases  roughly  3°C  larger  than  the  CMR  cases  in  the  first 
two  weeks  of  the  simulations.  The  CD  cases  also  show  a  substantially  slower  overall  rate 
of  smoke  removal  from  the  atmosphere  as  compared  to  removal  rates  in  the  CMR  cases. 
This  has  a  much  greater  relative  impact  on  subtropical  acute  temperature  responses  than 
on  mid  to  high  latitude  temperature  declines.  1-D  radiative/convective  modeling  and  3-D 
GCM  modeling  show  the  ameliorative  effect  of  smoke  infrared  opacity  generally  becomes 
larger  with  increasing  height  of  smoke  injection  for  smoke  injection  cloud  tops  in  the  low 
to  middle  troposphere.  Variation  of  smoke  injection  cloud  tops  in  the  middle  to  upper 
troposphere  produces  little  variation  in  the  surface  temperature  effect  of  smoke  infrared 
opacity.  Increasing  the  solar  “blackness”  of  smoke  by  increasing  the  fractional  carbon 
content  while  holding  the  total  carbon  mass  fixed  is  shown  to  reduce  the  ameliorative 
effect  of  smoke  infrared  opacity. 


Analysis  of  Surface  Temperature  Variability  and  Extremes 
in  GCM  Simulations  of  Large  Smoke  Injections 


Star  ley  L.  Thompson 

National  Center  for  Atmospheric  Research 
P.O.  Box  3000,  Boulder,  Colorado  80801 


Simulations  of  massive  nuclear  war  smoke  injections  performed  with  a  global  general 
circulation  model  (GCM)  have  been  analyzed  to  determine  some  statistics  of  surface  tem¬ 
perature  variability  and  extremes.  In  the  context  of  environmental  effects,  such  analyses 
may  be  more  meaningful  than  the  traditional  reliance  on  temporal  or  spatial  means  because 
biological  processes  are  particularly  susceptible  to  environmental  variability.  Global  maps 
of  coldest  temperature  reached  at  any  time  during  an  integration  imply  a  more  severe,  but 
more  realistic,  picture  of  potential  environmental  effects  than  either  daily  “snapshots”  or 
time  averages.  Probability  distributions  of  temperatures  for  selected  geographic  regions 
highlight  regional  differences  as  well  as  providing  probability  estimates  for  temperatures 
below  arbitrary  threshold  values.  The  enhanced  information  content  in  the  probability 
distributions,  as  compared  to  simple  means,  could  make  a  substantial  difference  in  envi¬ 
ronmental  inferences  drawn  from  the  data.  For  example,  some  regions  whose  time  mean 
temperatures  are  well  above  freezing  nonetheless  show  a  substantial  probability  of  near  or 
sub-freezing  events.  While  the  particular  analyses  currently  available  cannot  be  considered 
quantitatively  reliable  due  to  current  model  limitations,  the  analysis  methodology  appears 
promising  and  should  be  applied  to  future  GCM  simulations. 
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10  km  CD  180  Tg  Medium  -14.0  -15.3  -10.6 

10  km  CD  50  Tg  Sooty  -17.8  -17.0  -10.1 
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SECTION  3 


Emission  factor  measurements  for  Lodi  test  fires 
using  surface  based  emissions  sampling  equipment 

by 


Darold  E.  Ward  and  Colin  C.  Hardy 
Fire  and  Air  Resource  Management  Project 
USDA — Forest  Service 
4043  Roosevelt  Way  NE 
Seattle,  WA  98105 


ABSTRACT 

Surface  systems  of  towers  and  cables  were  used  for  supporting  sample  packages 
over  small  fires  of  about  0.5  hectares  In  size  preliminary  to  the  large  LODI 
fire  of  December  1986.  The  purpose  of  the  experiment  was  to  measure  emissions 
of  particulate  matter  and  gases  and  compare  emission  factors  for  these 
materials  to  the  measurements  acquired  through  using  airborne  sampling 
techniques.  The  subplots  were  located  to  represent  a  diversity  of  fuel 
conditions  of  the  larger  LODI  area:  chamlse,  standing  chaparral,  and  felled 
chaparral.  Burning  of  the  test  fires  was  accomplished  by  lighting  the  fire 
along  the  top  of  the  subunit  and  allowing  the  fire  to  back  downhill  for 
approximately  30  meters  and  then  lighting  the  fire  along  the  bottom  of  the  plot 
to  get  maximum  fire  Intensity  as  the  fire  headed  uphill  toward  the  location  of 
the  sampling  equipment.  Samples  were  collected  for  different  periods  of  the 
flaming  and  smoldering  combustion  processes  to  obtain  samples  representing  as 
early  as  possible  emissions  produced  by  each  of  these  processes.  In  addition, 
real  time  measurements  of  temperature  at  the  location  of  the  sample  packages 
and  the  vertlcle  velocity  were  measured  In  real  time.  The  concentration  of  CO 
and  C02  were  measured  In  real  time. 

Results  from  these  tests  are  compared  with  tests  of  emissions  from  logging 
slash  fires  In  the  Pacific  Northwest  states  of  Oregon  and  Washington.  Some  of 
the  results  for  the  LODI  subplot  emissions  study  are  summarized  In  the  table 
below: 


Emission  Factors 
PM2.5 
CO 

CHjCl 


8.3  g/kg 
61.5  g/kg 
.036  g/kg 


Percent  composition  of  PM2.5 
Cl 
Pb 
K 

Organic  carbon 
Elemental  carbon 


3.630 * 

0.339* 

8.419* 

51.4* 

8.9*  maximum  of  16.9* 
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Absorption  Measurements  from  Preliminary 
Site  Burns  in  Lodi  Canyon 


E.  M.  Patterson 

Georgia  Institute  of  Technology 
Atlanta,  Georgia 


D.  A.  Ward 

USDA  Forest  Service 
Seattle,  Washington 


We  have  measured  the  absorption  properties  of  smoke  samples 
collected  during  preliminary  burns  in  the  Lodi  Canyon.  These 
absorption  properties  are  expressed  in  terms  of  the  specific 
absorption  of  the  material.  These  measurements  show  a  range  of 
values  of  specific  absorption. 

We  will  present  these  data  and  will  discuss  these  results  in 
terms  of  the  range  of  variation  of  the  data,  the  relation  of  the 
absorption  data  to  the  elemental  carbon  concentrations,  and  the 
comparison  of  these  results  with  earlier,  similar  measurements. 
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Vertical  Plume  Velocities  in  a  Controlled  Burn. 

John  Hallett 

Desert  Research  Institute 
Reno,  Nevada  89506 

Time  lapse  photography  from  a  fixed  point  on  the  ground  were  made 
for  the  smoke  plumes  rising  from  the  Lodi  Canyon  controlled  burn  and  a 
controlled  JP-4  burn  in  Sandia  Albuquerque,  New  Mexico.  The  plume 
ascent  rate  was  measured  at  the  visual  periphery  of  a  rising  plume 
top,  and  of  individual  elements  following  puffing  oscillatory  motion 
from  the  source  (1/2  -  1  Hz).  In  the  case  of  the  Lodi  burn,  the  slant 
angle  was  3  km/15  km  ~  10°,  in  the  oil  fire  ~  15°,  so  error  of  order 
Cos  15°,  ~  5%  are  inherent  in  these  measurements.  If  we  assume  a 
simple  plume  or  thermal  model,  then  the  central  velocities  are 
approximately  twice  the  cap  velocity. 

Results 

Lodi  Canyon:  A  sequence  of  photographs  shows  a  plume  rising  with  top 
velocity  3.0  m  s-l  below  the  main  inversion  at  -  7000  ft.  As  the 
plume  passes  through  the  inversion  level  the  rise  velocity  falls  to 
1  m  s"l.  Fig.  1  decelerates  until  the  maximum  rise  is  -  plume  width 
(1  km)  and  then  sinks  back  with  a  velocity  of  order  of  the  rise  velo¬ 
city  to  a  level  somewhat  above  the  visual  (aerosol)  inversion.  The 
smoke  subsequently  moves  laterally  and  is  on  occasion  diluted  by 
mixing  in  a  Kel vin-Helmholtz  instability  wave  system  in  the  prevailing 
shear. 
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Sandia  Oil:  The  plume  rose  to  a  height  of  ~  1500  ft.  above  the  fire 
site,  and  was  carried  away  in  the  shear  towards  the  west.  The  cap  in 
velocity  (Fig.  4)  showed  a  maximum  about  halfway  up  the  plume 
(~  7  m  s-1)  being  ~  4  m  s"*  near  the  base  decelerating  to  4  m  s~*  near 
1000  ft.  and  falling  to  zero  at  1500  ft.  Lower  level  winds  are  a  few 
m  s*1  increasing  to  10-12  m  s“*  above  ~  1000  ft.  (Fig.  3). 

Visually,  the  strong  turbulence  visible  in  the  plume  above  the 
fire  and  in  the  cap  as  it  rose  into  the  environment  died  out  quickly 
as  the  plume  bent  over  and  was  carried  downstream  in  the  higher  velo¬ 
city  aloft.  The  "statistical"  plume  had  a  half  angle  expansion  of 
20  +  2°  during  its  initial  rise  stage,  and  gave  a  smoke  layer  carried 
downstream,  ~  400  ft.  deep  with  top  at  1500  ft.  Above  the  burn  site, 
it  was  estimated  that  overall  dilution  was: 

3 

vertical  volume  flux  at  200  m  =  20  x  20  x  5  m  _  50 

horizontal  volume  flux  at  500  m  100  x  100  x  10  m 

The  turned  over  cloud  was  everywhere  optically  thin  as  the  sun's  disc 

could  always  be  seen;  it  appeared  white,  indicating  a  large  particle 

size  for  the  soot,  *  urn. 
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Fig.  1:  Lodi  Canyon  burn  -  plume  cap  shape 


Fig.  3:  Sandia  Oil  pool  fire  velocity  variation  with  height. 

Note  the  necking  in  of  the  plume  just  above  the  fire. 


Paper  Submitted  for  Presentation  at  The  "Global  Effect?  Program" 
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Anomolous  Production  of  Sulfur  arid  Nitrogen  Species  in  the  Lodi  Fire 


Dean  A  Hegg,  Lawrence  F.  Radke,  Peter  V.  Hobbs  and  Charles  A.  Brock 
Department  of  Atmospheric  Sciences,  University  of  Washington 

Seattle  WA,  98195 
and 

Philip  J.  Riggan 

U.  S.  Forest  Service,  Riverside,  CA  92507 


ABSTRACT 

The  Lodi  burn  near  Los  Angeles  on  12  December  1986  produced 
unexpectedly  large  emissions  of  N0X,  SO2,  and  NO3”.  N0X  to  CO2  molar 

ratios  of  0.6  were  observed,  compared  to  the  value  of  2  X  10"^  adopted  by 
SCOPE  and  measured  by  us  in  prescribed  burns  of  forest  products  in  rural 
areas  of  Washington  State. 

It  is  postulated  that  the  high  emissions  of  N0X,  SO2  and  NC^-  in  the 

Lodi  burn  were  due  to  the  resuspension  of  previously  deposited  air 
pollutants  onto  the  vegetation.  As  such,  the  emissions  measured  in  the 
Lodi  burn  may  be  more  relevant  to  the  burning  of  biomass  near  urban  areas 
that  have  been  used  in  modeling  the  nuclear  winter  scenario. 


PLUME  CHARACTERIZATION  STUDIES  ON  THE  LODI  CANYON 
PRESCRIBED  BURN*  -  W.  Einfeld.  B.  D.  Zak,  B.  V.  Mokler  and 
D.  J.  Morrison,  Sandia  National  Laboratories **PO  BOX  5800, 
Albuquerque,  New  Mexico,  87185 


The  prescribed  chaparral  forest  bum  conducted  at  Lodi 
Canyon,  CA  in  December,  1986  provided  an  opportunity  for 
several  agencies  to  collaborate  in  a  major  plume 
characterization  effort.  The  Sandia  Twin  Otter  research 
aircraft  was  used  to  make  an  extensive  set  of  measurements 
on  plume  characteristics  both  at  the  head  of  the  fire  zone 
and  5-10  kilometers  downwind.  Instruments  flown  on  the 
aircraft  included  a  mobility  analyzer  with  a  measurement 
range  of  0.01  •  0.8  ftm  and  optical  particle  counters  with 
measurement  ranges  of  0.1  -  16  jim.  Plume  aerosol  samples 
were  collected  using  a  1  cubic  meter  grab  bag  sampler  and 
associated  filter  collection,  gas  sampling  and  pumping 
system.  Other  instruments  on  the  aircraft  for  this  fire 
included  a  NASA  sun  photometer  and  a  Desert  Research 
Institute  cloud  condensation  nuclei  counter. 

Results  to  be  reported  will  include  particle  size 
distributions,  aerosol  emission  factors  as  determined  by 
the  carbon  balance  technique,  and  smoke  optical  properties. 


★This  work  performed  at  Sandia  National  Laboratories,  Albuquerque,  New 
Mexico,  supported  by  the  United  States  Department  of  Energy  under  Contract 
DE-AC04-76DP00789 

**A  United  States  Department  of  Energy  Facility 
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SANDIA  TWIN  OTTER  AIRCRAFT  -  LODI  FIRE  SAMPLING  CONFIGURATION 

Particle  Size/Count 

PMS  ASASP-100-X . 0.1  -  3  uM 

PMS  FSSP-100-X . 0.5  -  16  uM 

TSI  DMPS . 0.01  -  0.8  uM 

TSI  Condensation  Nuclei  Counter 
DRI  Cloud  Condensation  Nuclei  Counter 
DRI  Formvar  replicator 

Aerosol  Chemisty  -  1  cubic  meter  grab  sampling  system 
Aerosol  mass,  aerosol  carbon,  gaseous  carbon 

Optical  Measurements 

MRI  NEPHELOMETERS  (bag  and  probe) 

NASA  SUN  TRACKING  PHOTOMETER 
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ANALYTICAL  METHODS 


Aerosol  carbon  -  Sunset  Labs 


furnace  combustion  /  reduction  /  FID  detection 


Gaseous  carbon  -  Sandia  Labs 

2  liter  grab  sample  /  GC  analysis  /  sonic  detector 

Aerosol  elemental  composition  -  Lawrence  Berkely  Lab 
X-ray  fluorescence  on  teflon  filters 
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LODI  FIRE  -  MEASUREMENT  SUMMARY 
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Emission  factor 


0.012  /  0.005 


EC/MASS  ratio 


0.089  /  0.032 


Sp  Mass  Scatt  Coeff  3.6  /  1.1 
(m2/g) 


Single  scatt  albedo  0.79  /  0.04 
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Measurements  of  Cloud  Active  Aerosol  in  Large  Scale  Fires. 

James  6.  Hudson 
Desert  Research  Institute 
Reno,  Nevada  89506 

CN  and  CCN  concentration  vs.  time  plots  are  shown  in  the  following 
10  figures.  These  measurements  were  made  in  the  Lodi  Canyon  brush 
fire  of  12  December  1986,  and  the  Sandia  Oil  fire  of  13  March  1987. 

It  is  noted  that  the  changes  in  the  CCN  concentrations  lag  the 
similar  changes  in  the  CN  concentrations  by  nearly  a  minute.  This  is 
approximately  the  time  required  for  the  sample  to  pass  through  the 
cloud  chamber  where  droplet  growth  occurs  before  entry  to  the  optical 
counter  where  recording  takes  place.  When  this  explanation  is  taken 
into  account  it  is  apparent  that  the  two  instruments  are  observing  the 
same  features  of  the  plumes.  Greater  resolution  is  observed  with  the 
CN  counter  simply  because  the  data  is  recorded  at  1  second  intervals 
instead  of  8  second  intervals.  Nevertheless,  the  time  widths  of  the 
plumes  as  observed  by  the  two  instruments  are  similar.  They  are 
generally  within  the  difference  in  time  resolution  (7  seconds). 
Therefore  the  CCN  instrument  does  not  appear  to  be  smearing  the  con¬ 
centration  in  time.  The  particles  are  being  recorded  with  the  same 
time  separation  which  they  had  in  the  sample  tube  when  they  went  into 
the  instrument.  Hence,  the  CCN  spectrometer  is  capable  of  better  time 
resolution  if  the  software  can  be  changed  to  allow  shorter  recording 
intervals.  This  is  now  currently  under  development.  This  would  be 
quite  valuable  because  this  is  the  only  instrument  with  multiple  chan¬ 
nel  data  which  can  be  resolved  on  such  a  time  scale. 
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Figures  1-7  for  the  Lodi  brush  fire  show  rather  high  CN  and  CCN 
concentrations.  The  most  surprising  aspect  of  these  measurements  is 
the  very  high  ratio  of  CCN  to  CN.  The  majority  of  the  plume  penetra¬ 
tions  show  a  ratio  of  nearly  one.  That  is  nearly  all  of  the  particles 
appear  to  be  CCN  active  at  1%  supersaturation.  This  ratio  is  quite 
variable  as  some  of  the  penetrations  display  a  much  lower  ratio. 
Differences  among  the  plume  penetrations  will  be  examined  to  see  if 
there  is  a  reason  for  these  variations  in  the  CCN/CN  ratio.  Neverthe¬ 
less,  it  appears  that  the  smoke  from  brushfires  is  perhaps  more  hygro¬ 
scopic  than  previously  thought. 

The  CCN  concentrations  which  are  presented  here  are  unprecedented. 
These  are  by  far  the  highest  CCN  concentrations  ever  measured.  CCN 
instruments  are  usually  limited  to  a  few  thousand  per  cubic  centimeter 
because  of  depletion  of  water  vapor  in  the  cloud  chamber.  The  main 
reason  this  chamber  can  measure  such  high  concentrations  is  that  it  is 
possible  to  reduce  the  sample  flow  rate  to  very  low  values 
(<5  cm3/min.)  and  still  measure  and  record  this  rate.  Apparently, 
this  technique  works  correctly  because  the  CN  measurements  indicate 
that  the  CCN  concentrations  could  not  possibly  be  any  higher. 
Therefore,  vapor  depletion  and  other  problems  associated  with  high 
count  rates  seem  to  have  been  overcome.  There  are  a  number  of  errors 
such  as  these  which  could  reduce  the  measured  CCN  concentrations  but 
the  only  known  error  which  would  raise  the  measured  concentration  is  a 
leak. 
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The  last  possibility  was  checked  during  one  of  the  plume  penetra¬ 
tions.  The  fifth  penetration  of  Figure  1  and  the  first  penetration  of 
Figure  4  display  a  zero  CCN  concentration.  This  is  as  it  should  be 
because  a  filter  was  placed  on  the  CCN  sample  inlet  during  this  pene¬ 
tration.  Therefore,  no  particles  were  getting  into  the  instrument 
except  those  which  were  measured  to  enter  through  the  sample  tube. 

The  last  penetrations  of  Figures  5  and  7  seem  to  display  CCN  con¬ 
centrations  higher  than  the  CN  concentrations.  However,  the  CCN 
sample  flow  rates  in  these  cases  slipped  to  such  low  values 
(<1  cm3/min.)  that  they  were  probably  not  accurate. 

Figures  8-10  show  CN  and  CCN  data  from  the  oil  pool  fire.  The 
most  obvious  difference  from  the  brush  fires  is  the  much  lower  par¬ 
ticle  concentrations.  The  next  important  contrast  with  the  brush 
fires  is  the  fact  that  the  CCN  concentrations  in  the  oil  fire  show 
very  small  increases  over  the  background  concentrations.  The  ratio  of 
CCN  to  CN  is  much  smaller  for  the  oil  fire.  There  the  oil  smoke  is 
much  less  hygroscopic  than  the  brushfire  smoke.  Those  results  are 
consistent  with  laboratory  work  presented  in  another  paper  in  these 
proceedings. 


Expansion  of  the  third  and  forth  penetrations  shown  in  Figure 


penetrations 


Fig.  5  for  later  penetrations.  The  last  penetration  shown  here  has  an 
pecially  strong  updraft  which  caused  some  minor  cabin  damage  and  interrup- 
on  of  the  instruments. 


Fig.  1  except  that  this  is  for  the  oil  pool  fire  of  13  March.  Shown  here 
background  concentrations  before  the  fire,  a  filter  on  the  sample  inlet 
monitor  the  instrument  noise,  and  three  penetrations  of  the  plume. 
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OPTICAL  DEPTH  MEASUREMENTS  OF  FOREST  FIRE  SMOKES:  TWO  CASE 
STUDIES. 

R. F. Pueschel,  P. B. Russell,  T. P. Ackerman,  D. A. Allen  and 
D. S. Colburn.  NASA  Ames  Research  Center,  Moffett  Field,  CA 
94035. 

B. D. Zak  and  W. Einfeld.  SANDIA  National  Laboratories, 
Albuquerque,  NM  87185. 


Since  installing  the  Ames  Tracking  Airborne  Sunphotometer  on 
the  Sandia  Laboratories'  Twin  Otter  research  aircraft,  we 
had  opportunities  to  measure  optical  depths  between  380  and 
1020  nm  in  smokes  from  a  fire  of  opportunity  near  Fresno,  CA 
on  11  Sep  1986,  and  from  the  Lodi  Canyon  controlled  burn  on 
12  December  1986. 

The  optical  depths  of  both  smokes  are  similar  in  magnitude 
and  wavelength  dependence.  At  short  wavelengths,  the  optical 
depths  of  the  smokes  are  increased  by  up  to  one  order  of 
magnitude  over  those  of  background  air.  The  increase  at  1020 
nm  wavelengths  i3  up  to  2  orders  of  magnitude  for  the  Lodi 
Canyon  fire,  but  less  than  one  order  in  the  smoke  of  the 
Fresno  fire  over  the  respective  background  optical  depths. 

The  optical  depth  in  background  air  in  the  vicinity  of 
Los  Angeles  differs  markedly  from  that  of  the  Central 
Valley.  A  strong  wavelength  dependence  suggests  the 
prevalence  in  the  Los  Angeles  atmosphere  of  very  small 
particles  and  fewer  large  aerosols  than  are  found  in  central 
California. 

Aging  of  the  Lodi  Canyon  controlled  burn  smoke  indicates  a 
decrease  of  optical  depths  due  to  dilution,  and  a  wavelength 
dependence  that  suggests  an  alteration  of  the  smoke  particle 
size  distribution;  the  aged  smoke  aerosol  contains  fewer 
small  particles  and  has  a  larger  modal  particle  size  than 
does  the  fresh  smoke. 
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POLYTECHNIC  OF  NEW  YORK 
ATTN:  BJ  BULKIN 
ATTN:  GTESORO 

PRINCETON  UNIVERSITY 

ATTN:  J  MAHLMAN 

QUADRI  CORP 

ATTN:  H  BURNSWORTH 

R  &  D  ASSOCIATES 

ATTN:  AKUHL 
ATTN:  F  GILMORE 
ATTN:  G JONES 
ATTN:  J SANBORN 
ATTN:  RTURCO 

R  &  D  ASSOCIATES 

ATTN:  B  YOON 
ATTN:  C  KNOWLES 

R J  EDWARDS  INC 

ATTN:  R  SEITZ 

RADIATION  RESEARCH  ASSOCIATES,  INC 
ATTN:  B  CAMPBELL 
ATTN:  M  WELLS 

RAND  CORP 

ATTN:  P  ROMERO 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  J  KELLEY 

SCUBED 

ATTN:  B  FREEMAN 
ATTN:  KDPYATT,  JR 
ATTN:  R  LAFRENZ 

SCIENCE  APPLICATIONS  INC 
ATTN:  REDELMAN 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  CHILL 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  B  MORTON 


ATTN:  B SCOTT 
ATTN:  D SACHS 
ATTN:  GT  PHILLIPS 
ATTN:  JBENGSTOM 
ATTN:  D  HAMLIN 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  D BACON 
ATTN:  DRLGOURE 
ATTN:  F  GIESSLER 
ATTN:  J COCKAYNE 
ATTN:  JGAHAN 
ATTN:  J SHANNON 
ATTN:  J STUART 
ATTN:  MSHARFF 
ATTN:  WLAYSON 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  J  SONTOWSKI 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  T  HARRIS 

SCIENTIFIC  RESEARCH  ASSOC,  INC 
ATTN:  B  WEINBERG 

SRI  INTERNATIONAL 

ATTN:  CWITHAM 
ATTN:  D GOLDEN 
ATTN:  D  MACDONALD 
ATTN:  D ROBERTS 
ATTN:  EUTHE 
ATTN:  GABRAHAMSON 
ATTN:  J  BACKOVSKY 
ATTN:  W  CHESNUT 

STAN  MARTIN  AND  ASSOCIATES 
ATTN:  SB  MARTIN 

STANTON  CONSULTING 

ATTN:  M STANTON 

SWETL,  INC 

ATTN:  TY  PALMER 

SYSTEM  PLANNING  CORP 
ATTN:  B  GARRETT 
ATTN:  C  FELDBAUM 
ATTN:  J  SCOURAS 
ATTN:  MBIENVENU 
ATTN:  R  SCHEERBAUM 

SYSTEMS  AND  APPLIED  SCIENCES  CORP 
ATTN:  M  KAPLAN 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  AORTELL 
ATTN:  F  LEOPARD 
ATTN:  J  FORD 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  DGUICE 

TEXAS  ENGR  EXPERIMENT  STATION 
ATTN:  WH  MARLOW 
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TOYON  RESEARCH  CORP 
ATTN:  CTRUAX 
ATTN:  J  GARBARINO 
ATYN:  J  ISE 

TRW  INC 

ATTN:  HBURNSWORTH 
ATTN:  JBELING 

TRW  INC 

ATTN:  FFENDELL 
ATTN:  G  KIRCHNER 
ATTN:  H  CROWDER 
ATTN:  JFEDELE 
ATTN:  M  BRONSTEIN 
ATTN:  R  BACHARACH 
ATTN:  S  FINK 
ATTN:  T  NGUYEN 

TRW  SPACE  &  DEFENSE,  DEFENSE  SYSTEMS 
ATTN:  M  HAAS 

VISIDYNE,  INC 

ATTN:  H  SMITH 
ATTN:  J CARPENTER 

WASHINGTON,  UNIVERSITY  OF 
ATTN:  J  I  KATZ 

FOREIGN 

AERE  ENVIRONMENTAL  AND  MEDICAL  SC 
ATTN:  SPENKETT 

ATOMIC  WEAPONS  RSCH  ESTABLISHMENT 
ATTN:  PFA  RICHARDS 

ATOMIC  WEAPONS  RSCH  ESTABLISHMENT 
ATTN:  DL  JONES 
ATTN:  DM  MOODY 

AUSTRALIA  EMBASSY 
ATTN:  ATEBBS 

ATTN:  AIR  VICE  MARSHAL  BGRATION 
ATTN:  DR  A  J  BEDFORD 

BRITISH  DEFENCE  STAFF 
ATTN:  J  EDMONDS 

CANADIAN  FORESTRY  SERVICE 
ATTN:  B STOCKS 
ATTN:  TLYNHAM 

CSIRO 

ATTN:  IGALBALLY 

CSIRO:  ATMOSPHERIC  RESEARCH 
ATTN:  APITTOCK 

EMBASSY  OF  BELGIUM 

ATTN:  LARNCULD 

ISRAEL  EMBASSY 

ATTN:  NBELKIND 


MAX-PLANCK  INSTITUTE  FOR  CHEMISTRY 
ATTN:  P  J  CRUTZEN 

MINISTRY  OF  DEFENCE 
ATTN:  R  RIDLEY 

NATIONAL  DEFENCE  HEADQUARTERS 
ATTN-  HAROBITALLE 

TRINITY  COLLEGE 

ATTN:  F  HARE 

DIRECTORY  OF  OTHER 

ATMOS.  SCIENCES 

ATTN:  GSISCOE 

BUCKNELL  UNIVERSITY 

ATTN:  0  ANDERSON 

CALIFORNIA,  UNIVERSITY 

ATTN:  R  WILLIAMSON 

CALIFORNIA,  UNIVERSITY  OF 
ATTN:  LBADASH 

COLORADO,  UNIVERSITY  LIBRARIES 
ATTN:  JBIRKS 
ATTN:  R  SCHNELL 

DREXEL  UNUVERSITY 

ATTN:  J  FRIEND 

GEORGE  MASON  UNIVERSITY 
ATTN:  PROFS  SINGER 
ATTN:  R  EHRLICH 

GEORGE  WASHINGTON  UNIVERSITY 
ATTN:  R  GOULARD 

GEORGIA  INST  OF  TECH 

ATTN:  E PATTERSON 

HARVARD  COLLEGE  LIBRARY 
ATTN:  W  PRESS 

HARVARD  UNIVERSITY 

ATTN:  G  CARRIER 

HARVARD  UNIVERSITY 

ATTN:  DEARDLEY 

IOWA,  UNIVERSITY  OF 
ATTN:  SPYNE 

MARYLAND  UNIVERSITY  OF 
ATTN:  AROBOCK 
ATTN:  AVOGELMANN 
ATTN:  RELLINGSON 

MIAMI  LIBRARY  UNIVERSITY  OF 
ATTN:  C  CONVEY 

MIAMI  UNIV  LIBRARY 

ATTN:  J  PROSPERO 
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NEW  YORK  STATE  UNIVERSITY  OF 
ATTN:  R  CESS 

OAK  RIDGE  ASSOCIATED  UNIVERSITIES 
ATTN:  C  WHITTLE 

SOUTH  DAKOTA  SCH  OF  MINES  &  TECH  LIB 
ATTN:  H  ORVILLE 

TENNESSEE,  UNIVERSITY  OF 
ATTN:  K  FOX 

UNIVERSITY  OF  SOUTH  FLORIDA 
ATTN:  SYING 

UNIVERSITY  OF  WASHINGTON 
ATTN:  CLEOVY 
ATTN:  LRAOKE 
ATTN:  P  HOBBS 

VIRGINIA  POLYTECHNIC  INST  LIB 
ATTN:  MNADLER 

WISCONSIN  UNIVERSITY  OF 
ATTN:  P  WANG 
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